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NOTICE TO AUTHORS 


1. Communications.—Papers must be communicated to the Society by a Fellow. They 
should be accompanied by a summary at the beginning of the paper conveying briefly the 
content of the paper, and drawing attention to important new information and to the main 
conclusions. The summary should be intelligible in itself, without reference to the paper, 
to a reader with some knowledge of the subject; it should not normally exceed 200 words 
in length. Authors are requested to submit MSS. in duplicate. These should be 

using double spacing and leaving a margin of not less than one inch on the 
left-hand side. Corrections to the MSS. should be made in the text and not in the 
margin. Unless a paper reaches the Secretaries more than seven days before a Council 
meeting it will not normally be considered at that meeting. By Council decision, MSS. of 
accepted papers are retained by the Society for one year after publication; unless their 
return is then requested by the author, they are destroyed. 


2. Presentation.—Authors are allowed considerable latitude, but they are requested to 
follow the general style and arrangement of Monthly Notices. References to literature 
should be given in the standard form, including a date, for printing either as footnotes or in 
a numbered list at the end of the paper. Each reference should give the name and 
initials of the author cited, irrespective of the occurrence of the name in the text (some 
latitude being permissible, however, in the case of an author referring to his own work). 
The following examples indicate the style of reference appropriate for a paper and a book, 
respectively :— 

A. Corlin, Zs. f. Astrophys., 15, 239, 1938. 

H. Jeffreys, Theory of Probability, 2nd edn., section 5.45, p. 258, Oxford, 1948. 


3. Notation.—For technical astronomical terms, authors should conform closely to the 
recommendations of Commission 3 of the International Astronomical Union (Trans. 
I.A.U., Vol. VI, p. 345, 1938). Council has decided to adopt the I.A.U. 3-letter abbrevi- 
ations for constellations where contraction is desirable (Vol. IV, p. 221, 1932). In general 
matters, authors should follow the recommendations in Symbols, Signs and Abbreviations 
(London : Royal Society, 1951) except where these conflict with 1.A.U. practice. 


4. Diagrams.—These should be designed to appear upright on the page, drawn 
about twice the size required in print and prepared for direct photographic 
reproduction except for the lettering, which should be inserted in pencil. 
Legends should be given in the manuscript indicating where in the text the 
figure should appear. Blocks are retained by the Society for 10 years; unless the author 
requires them before the end of this period they are then destroyed, 


5. Tables—These should be arranged so that they can be printed upright on 
the page. 


6. Proofs.—Costs of alterations exceeding 5 per cent of composition must be borne by 
the author. Fellows are warned that such costs have risen sharply in recent years, and it 
is in their own and the Society’s interests to seek the maximum conciseness and simplifi- 
cation of symbols and equations consistent with clarity. 


7. Revised Manuscripts—When papers are submitted in revised form it is especially 
requested that they be accompanied by the original MS. 


Reading of Papers at Meetings 


8. When submitting papers authors are requested to indicate whether they will be 
willing and able to read the paper at the next or some subsequent meeting, and approxi- 
mately how long they would like to be allotted for speaking. 


g. Postcards giving the programme of each meeting are issued some days before the 
meeting concerned. Fellows wishing to receive such cards whether for Ordinary 
Meetings or for the Geophysical Discussions or both should notify the Assistant Secretary. 
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MEETING OF 1953 JANUARY 9 
Professor H. Dingle, President, in the Chair 


The election by the Council of the following Fellows was duly confirmed :— 


*Alan Henry Batten, St. Regulus Hall, St. Andrews, Fife (proposed by 
E. Finlay Freundlich) ; 
*Joyce Gertrude Billings, B.Sc., 43 Rockland Road, Downend, Bristol, 
Somerset (proposed by W. H. McCrea); 
Eva Cassirer, Department of History and Philosophy of Science, University 
College, London, W.C.1 (proposed by H. Dingle); 
Leonard William Clarke, M.A., 115 Swakeley’s Road, Ickenham, Middlesex 
(proposed by 'T. G. Cowling); 
*Donald George Ewart, B.Sc., 15 Lilburn Gardens, South Gosforth, 
Newcastle-on-Tyne, 3 (proposed by W. M. Smart); 
Randolph D. Gifford, D.Sc., A.M.I.E.E., 15 Willington Road, Birmingham, 
20 (proposed by M. C. Johnson); 
Leonard John Heath, B.Sc., 20 Wren Gardens, Dagenham, Essex (proposed 
by C. A. Ronan); 
Samuel Herrick, B.A., Ph.D., Department of Astronomy, University of 
California, Los Angeles, California, U.S.A. (proposed by D. H. Sadler); 
Jan Hospers, Pembroke College, Cambridge (proposed by S. K. Runcorn); 
*Angela Margaret James, H.M. Nautical Almanac Office, Herstmonceux 
Castle, near Hailsham, Sussex (proposed by D. H. Sadler); 
*Leon Mestel, Department of Mathematics, The University, Leeds (proposed 
by T. G. Cowling); 
Jack Miller, 10 Bury Street, St. James’s, London, S.W.1 (proposed by 
M. Davidson); 
*Peter Coles Price, B.A., Trinity College, Cambridge (proposed by R. O. 
Redman); 
Michael John Seaton, Ph.D., 26 Kilmorie Road, London, S.E.23 (proposed 
by H. S. W. Massey); and 
*Robert Wilson, B.Sc., Royal Observatory, Edinburgh (proposed by W. M. H. 
Greaves). 
* Transferred from Junior Membership. 





Meeting of 1953 January 9 


The election by the Council of the following Junior Members was duly 
confirmed :— 


Gilbert Elliott Satterthwaite, 48 Upper Abbey Road, Belvedere, Kent 
(proposed by P. J. D. Gething); and 

Paul Anthony Young, 125 Saltram Crescent, London, W.g (proposed by 
R. W. B. Pearse). 


Ninety presents were announced as having been received since the last 
meeting. 


The President announced that the Council had awarded the Gold Medal 
of the Society to Professor Subrahmanyan Chandrasekhar, of Yerkes Observatory, 
Williams Bay, Wisconsin, U.S.A., for his contributions to mathematical 
astrophysics. 


The President announced that the Council had awarded the Eddington Medal 
to Canon Georges Lemaitre, of the University of Louvain, for his work on the 
expansion of the universe. 


The President announced that the Council had awarded a Jackson-Gwilt 
Medal and Gift to Mr John Philip Manning Prentice for his contributions to the 
study of meteors. 








AN INVESTIGATION OF THE HII REGIONS BY A RADIO METHOD 


P. A. G. Scheuer and M. Ryle 


(Received 1952 November 17*) 


Summary 


Observations have been made to determine the effect of the H II regions 
on the radio emission from the galaxy. Since both interstellar hydrogen 
and early-type stars show strong galactic concentration, regions of ionized 
hydrogen are likely to produce the greatest effect in directions near the 
galactic equator. At low frequencies, where the ‘“‘ brightness temperature ”’ 
of the general galactic radiation exceeds the kinetic temperature of the 
interstellar hydrogen, the absorption due to the hydrogen will exceed its 
emission, and a dark band is to be expected along the galactic equator. At 
frequencies greater than about 60 Mc/s, where the brightness temperature 
of the galactic radiation falls below the gas temperature, this band should 
appear in emission. 

An interferometric method of large resolving power has been used to 
determine the distribution of brightness near the galactic plane at 
frequencies of 81-5 Mc/s and 210 Mc/s. In each case a bright band about 2° 
wide has been observed superimposed on the general radiation. No accurate 
observations have yet been made at frequencies less than 60 Mc/s. 

A comparison of the observed brightness temperature of the band at the 
two frequencies makes it possible to derive a preliminary estimate of the 
kinetic temperature and effective density in the HII regions. The 
temperature derived in this way is considerably greater than the accepted 
value. 

Observations at different galactic longitudes have suggested that the 
distribution of the H II regions shows a marked concentration towards the 
centre of the galaxy. 





1. Introduction.—The earliest theories of the radio-frequency radiation from 
the galaxy were based on emission from the interstellar gas (1, 2). Later observa- 
tions at lower frequencies showed that it would be necessary to assume a mean 
interstellar temperature and density much greater than the values derived from 
the visual observations if the results obtained at frequencies between 30 and 
100 Mc/s were to be explained on this hypothesis (3). The discovery of radio 
stars suggested an alternative explanation of the galactic radiation (4, 5) and later 
work has been based on the assumption that both mechanisms may contribute to 
the total radiation. 

If the galactic radiation were due solely to radio stars, and the effects of the 
interstellar matter were negligible, there would be no reabsorption and the galaxy 
would show the characteristics of a region which was optically thin. The spectrum 
of the radiation would then be the same in all directions, unless there were popula- 
tions of radio stars having different spectra and different distributions within the 
galaxy. 

If the optical depth of the interstellar matter were appreciable, two effects 
would modify this result; firstly, the radiation from distant radio stars would 
suffer absorption, and secondly, the emission from the interstellar matter might 


* Received in original form 1952 July 29. 
1” 
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make an appreciable contribution to the total radiation. If it were assumed that 
there were only one population of radio stars, it would therefore be possible to 
detect any effects of interstellar matter by observing the spectrum of the galactic 
radiation in a number of different directions. 

Early observations by Moxon (6) indicated that between frequencies (v) of 
45 and 200 Mc/s, the brightness temperature of the radiation from near the galactic 
poles varied as »-*', whilst near the galactic centre the brightness temperature 
varied as v-*7, These results, though they were of great interest, were not con- 
clusive, owing to possible errors involved in making such observations with aerial 
systems of low resolving power. 

Inamore recent paper Piddington (7) has summarized the results of anumber of 
measurements which have been made to determine the distribution of intensity 
with direction at frequencies ranging from 18 to 3000 Mc/s. According to his 
results the contours of constant intensity near the galactic plane show a progressive 
broadening at the lower frequencies; if the derived contours at 64 and 18 Mc/s are 
correct it would be necessary to postulate that the optical depth of the interstellar 
gas was appreciable even in directions 15° away from the galactic equator.* 

Westerhout and Oort (9), on the other hand, have concluded that any interstellar 
emission and absorption processes must be associated with free-free transitions in 
regions of ionized hydrogen in the neighbourhood of early-type stars; in all other 
interstellar regions either the temperature or the density is too low for appreciable 
radio-frequency emission. These ionized hydrogen regions have been observed 
by their Balmer-line emission spectra (10) within a few kiloparsecs of the Sun, and 
the observations have shown that most of the regions are within about 200 parsecs 
of the galactic plane. If radio observations are made from the Earth, the total 
optical depth of the galaxy due to such regions should therefore decrease rapidly 
with increasing galactic latitude; on present evidence the optical depth might be 
expected to fall to half at about 1° from the galactic equator. Because of its 
limited extent Westerhout and Oort concluded that the interstellar matter would 
be unlikely to produce an effect sufficient to be detected with the comparatively 
poor resolving power which was used in the surveys they considered. 

The present paper describes an investigation which has been made to detect 
the presence of a distribution of regions of ionized hydrogen similar to that 
considered by Westerhout and Oort. An interferometric method of large 
resolving power has been used to obtain detailed measurements of the intensity 
distribution near the galactic plane, at frequencies of 81-5 and 210 Mc/s; _ the 
results of these measurements have then been compared with the predicted 
distribution of intensity, taking into account the general galactic radiation, and the 
emission and absorption associated with the regions of ionized hydrogen. 

2. Summary of theoretical analysis —The intensity of radio waves from a given 
direction depends on the emission from radio stars, the interstellar gas, and extra- 
galactic sources of radiation, and on the absorption of all components of the radia- 
tion in the interstellar gas. Westerhout and Oort have shown that the main part 
of the galactic radiation is consistent with emission from radio stars, distributed 
according to the model of Oort and van Woerkom (11). It is, however, possible 
that ‘‘radio stars”’ are in fact situated outside the galaxy, and that some other 


* Preliminary observations (8) at frequencies of 72 and 158 Mc/s, using the larger resolving 
power of the 220-ft. paraboloid at Jodrell Bank, have indicated that there is little difference 
between the contours at these two frequencies in a region between galactic longitudes 60° and 120°. 





No. 1, 1953 An investigation of the H II regions by a radio method 5 


mechanism is responsible for the radiation from the galaxy. In that case, the 
analysis of this section would still be applicable, provided that the new mechanism 
could be regarded as being “‘opticallythin”. From considerations such as those 
of Ryle (12) and Hutchinson(13) it seems likely that any interstellar mechanism 
for galactic radiation, other than that of free-free transitions, would be associated 
with particles of very great energy, and would therefore satisfy the requirement of 
small optical depth. For the present purposes it will be convenient to assume 
that radio stars are in fact responsible for the main component. 

In order to explain the intensity distribution observed by Bolton and Westfold 
(14) at 100 Mc/s, Westerhout and Oort further postulated the existence of an 
isotropic component of radiation, having an intensity about 10 per cent of that 
observed in the direction of the galactic centre. It is not yet certain whether this 
component represents a true extragalactic contribution, or whether the observations 
could be accounted for by a more spherical distribution of radio stars; some 
evidence for the former hypothesis has recently been obtained by Hanbury 
Brown and Hazard (15). 

The effect of the HII regions on the intensity distribution will now be 
considered. ‘The absorption coefficient for a frequency v in a region containing 
N, electrons, and N, ions of atomic number Z, per cm’, is given by Smerd and 
Westfold (16) : 

— 4e°A,(2)N,.N,Z* 
" 3(2n)'*(mkTYPCr*p,’ 


where y, is the refractive index (which may be taken as I in all the regions discussed 
here), 7, is the kinetic temperature in the region, and 


A,(2)=In {1 + (4kT,/eN,1°)?}. 


Following Westerhout and Oort, Z?N,/N, is taken to be 1-1. For T,=10‘deg. K 
and N,=8, A,(2) =32°6. If N, is changed by a factor of 10 or T, by a factor of 2, 
A,(2) is only changed by about 5 per cent, and A,(2) is therefore assumed to be 
32-6 throughout the following calculations. 

The absorption coefficient can now be written 





K,=0°177 T-*? v?N,?. 


It can be shown that, when the optical depth of a single H II cloud is small, it is 
permissible to treat a region containing many clouds as a continuous medium with 
an effective absorption coefficient given by substituting, in the above formula, 
the mean value of N,? over the clouds and the spaces between them. Westerhout 
and Oort estimate that the optical depth of a single cloud is 0-02 at 100 Mc/s, so 
that at 30 Mc/s it would be about 0-2; this approximation is therefore certainly 
valid for frequencies greater than 30 Mc/s. 

Now suppose that, in a particular direction, the brightness temperature due to 
the radio stars in the absence of absorption is 7,, and that the contribution to this 
temperature from radio stars situated at distances between s and s+ds is t,ds. 
Suppose, also, that the extragalactic component produces a brightness temperature 
T,'. If the kinetic temperature (7',) is assumed to be uniform throughout all the 
HII regions, the actual brightness temperature T is given by 


T=T, Ly —exp (- - «,ds)} + I. t,exp (- « ds) ds + T,’ exp (- | K, as) . 
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In the absence of absorption and emission by the HII regions, the brightness 
temperature would be simply T7*=7,+7,'. The presence of the HII regions 
thus increases the observed brightness temperature by an amount 


(T-T*)=T, {x —exp (- [- «,ds)} -|" tr —exp (- [eds,)} ds 


‘ Ty {1—exp(~|"syds) (1) 


The first and third terms of the right-hand member of equation (1) depend only 


on7,= | «, ds, the total optical depth of the H II regions in the direction considered, 
~ 0 


but the second term (the absorption of radio-star radiation) also depends to some 
extent on the distribution of the HII clouds and the radio stars along the line 
of sight. For calculating this term, it will be assumed initially that the HII 
regions are distributed uniformly over a disk of radius 12°5 kps, and that the radio 
stars are distributed according to a simplified version of Oort and van Woerkom’s 
model, due to de Kort (17).* j 

Equation (1) can now be expressed in the form 


(T—T*)=T, {t—exp(—7,)}—f(7,) T*, (2) 


where f (7,) is a known function of 7, for a given galactic longitude. 
Using the constants adopted by Strémgren (19), 
T.=10'deg. K, [N2ds=1-9 x 10% cm", 

the emission from the HII regions (first term of equation (2)) and the absorption 
of the radio-star and extragalactic components (second term of equation (2)) 
have been computed for a direction near the galactic centre. ‘These curves, 
together with their sum, which represents the overall effect of the HII regions, 
(T—T*), are shown in Fig. 1. It can be seen that with these constants, the HII 
regions should be visible in emission at frequencies greater than 60 Mc/s, and 
in absorption at lower frequencies. At other galactic longitudes the transition 
between emission and absorption should take place at a lower frequency. 

If measurements of the brightness temperature at a given galactic longitude 
were made at different galactic latitudes, a curve like Fig. 2 would be obtained, from 
which both 7* and (7— 7*) could be estimated. Using equation (2), a series of 
theoretical curves could be constructed relating the value of (T—7*) to T, and 
{N2ds. Fora different frequency or a different value of T* another set of curves 
would be obtained : examples of such curves (for 7* =6 500 deg. K at 81-5 Mc/s 
and 7'* =750 deg. K at 210 Mc/s) are shown in Fig. 3. From the point of inter- 
section of the curves representing the observed values of (7 — T*) at two different 
frequencies, it would then be possible to determine both 7, and {N,* ds for the 
direction of observation. 

If, now, such observations were made at a number of galactic longitudes, it 
would become possible to obtain information on the actual spatial distribution of 
the HII regions; it would then no longer be necessary to assume the model of a 
uniform disk. 

* In a recent paper Mills (18) concludes that there is a significant concentration of intense radio 
stars towards the galactic equator, and suggests that these may represent a new class of radio stars 


showing strong galactic concentration. The significance of these observations in relation to the 
present investigation will be discussed in Section 7. 
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Fic. 1.—Estimated changes in brightness temperature due to HII regions, expressed as 
fractions of T*, the brightness temperature which would be observed in the absence of H II regions. 


(a) Emission from H II regions. 
(b) Absorption of “‘ radio-star’’ and “ extragalactic ’’ components. 
(c) Net change due to H II regions. 








Fic, 2.— Sketch of the expected variation of brightness temperature with galactic latitude. 
Ordinates: Brightness temperature. 
Abscissae: Galactic latitude (b). 
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3. The experimental difficulties—Before describing the present method of 
observation, it may be useful to indicate why previous surveys of galactic radiation 
have failed to disclose the existence of a feature 2° wide, and probably more than 
go° long, having a brightness temperature perhaps twice as great as that of 
neighbouring areas of the sky. The two most complete surveys of galactic 
radiation (14, 20) were made with aerial systems having apertures of the order of 
four wave-lengths, giving reception patterns whose angular widths were about 15. 
The effect of the finite width of the reception pattern has been carefully corrected 
in these surveys, but it is possible to show that such methods cannot yield inform- 
ation on angular structure smaller than about 1/2” radians, where nd is the effective 
aperture of the aerial. Only two surveys appear to have been made with a 
sufficiently large resolving power: one of them, reported by Lovell and Clegg (8), 
made use of the 220-ft. paraboloid at Jodrell Bank, but was restricted to areas within 
15° of the zenith, and therefore gave no information on regions near the galactic 
centre; the other was made by Piddington and Minnett (21), at a frequency of 
1 200 Mc/s, in a restricted region near the galactic centre. 








0 10000 20 000 
Fic. 3.—Curves relating the value of (T—T*) to physical properties of the H II regions. 
Ordinates: T,~** [ Ng® ds x 10~*. 
Abscissae: Kinetic temperature of hydrogen clouds (Te). 
Curves for 81-5 Mc/s: —————— 
Curves for 210 Mc/s: — —— — 
Each curve is labelled with the corresponding value of (T—T*). 


It can be seen from Fig. 1 that an investigation of the HII regions is likely to 
involve accurate measurements of the intensity distribution at wave-lengths of 
several metres; at these wave-lengths there are serious practical difficulties in 
constructing aerials of sufficient aperture. 

The possibility of determining the distribution of brightness across a source by 
observations with an interferometer of variable spacing was first pointed out by 
Michelson and Pease (22), and the application of this principle to radio astronomy 
was suggested by McCready, Pawsey and Payne-Scott (23). ‘The method has 
been used by Stanier (24), Machin (25), and O’Brien (26) for determining the 
distribution of brightness across the undisturbed solar disk. 
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In the present problem a large resolving power is needed to determine the 
distribution of brightness temperature in galactic latitude, but the expected 
variation with galactic longitude is comparatively small; the observations could 
therefore be made with an aerial system having a relatively small resolving power 
in directions parallel to the galactic plane. 

The principle of the method which has been used will be discussed in the next 
section; the practical details are described in Section 5. 

4. The variable spacing interferometer.—Suppose that two aerials are mounted 
with a spacing of n wave-lengths, on an axis which passes through the galactic 
poles, and that they are directed at a point on the galactic equator (6 =0) whose 
galactic longitude is 4,. Let the reception pattern of each aerial be given by 
F(Z, 6). 

If the distribution of brightness is represented by the brightness temperature 
T (i, 6), the power intercepted by each aerial from an element of solid angle dw at 
(/, 6) may be represented by an effective aerial temperature 

T (1, 6) dw f (1, 5) 
bT,= , 
SJ (4 5) cos bdl db 
The denominator of this expression may be written as 47/G, where G is the gain of 
each aerial over an isotropic radiator. If the two aerials are now used in a phase- 
switching interferometer system (27), the deflection on the record due to the small 
element is given by 


CT ,=C g T (1, 6) f (1, 6) 8 cos (27n sin b), 





where C is a constant depending on the deflection sensitivity of the recording 


equipment. 
The deflection produced by the radiation from the whole sky may then be 
written 


d(n) -c= | T (1, 6) f (1, 6) cos (27m sin b) dw. (3) 
If the variation of 7 (/, b) with /can be regarded as linear within the reception pattern 


of the aerials, and if f(/,5) can be expressed in the form L(/) B(6), where L (J) is 
symmetrical about /=/,, equation (3) can be rewritten 


d(n=C © |" L(l)dl ¥ T (1g, 6) B(b) cos (27 sin b) cos b db 


-c2 | L(ladl | X (x) cos 27nx dx, 


x=sin 8, 
x)= {7D BO [x]<x 
or 0 |x |>r1. 


By using a number of different aerial spacings n, d(m) can be determined as a 
functionofm. Ifthe reception pattern of each aerial and the distribution in galactic 
latitude of the brightness temperature are both symmetrical, X(x) is an even 
function of x, and it follows from (4) that 


& ° 
— d(n) cos 27nx dn. 


“7 


CG L(l)dl°° 


7 ae 


T (4, 6) B(6) =X (x) = 
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Hence it is possible to derive the required distribution of brightness 
temperature. 

When the distribution of T is required only for values of 6 up to a few degrees, it 
is sufficiently accurate to put x =sinb =6 and B(b) =1. 

If it is not assumed that X (x) is an even function, it can still be expressed as the 
sum of an even and an odd function; the former will be given by (5), and the latter 
by the corresponding Fourier sine transform. ‘This sine transform could, in 
principle, be found by inserting an extra quarter wave-length of cable into one side 
of the interferometer; the cosines in (3), (4) and (5) then all become sines, and the 
right-hand side of (5) gives the antisymmetrical part of X (x). 

In practice, however, both the cosine and the sine transform of X (x) can usually 
be obtained from the same record in the following way. ‘The above calculations 
refer to the deflection d(m) at the instant when the interferometer axis passes 
through the galactic poles; owing to the rotation of the Earth the aerials will, after 
a short time, no longer be directed to the point (4, 0), but to (4,+4, 5,), and 
the deflection is now 


ce | T (1,6) f(l—1,, 6 —b,) cos (2mn sin (b—b,)) dw. (3a) 


Provided /, and 5, are small they can be neglected in all factors except the last ; 
and furthermore, for small 5, sin (b—5,) =b—5,. 
The deflection is therefore approximately 


c= f’ L(nat| [" T (h, b) Bb) cos2nnbdb | cos 27nb, 


+ | ‘* T (lb) B (b) sin 2nnb db | sin 2nnb,. (4a) 


/—2 


‘Thus except for very small values of the variation with time is approximately 
sinusoidal, and the amplitude and phase of the record enable both Fourier 
transforms of 7'(/,, 5) B(b) to be found. 

5. The experimental method. (a) The aerial system.—In the method des- 
cribed above, the resolving power in galactic latitude does not depend on the size 
of the aerials, but simply on the greatest spacing used. In order to obtain 
information on the distribution in galactic longitude, the aerials were made with 
the greatest practicable aperture in a direction parallel to the galactic plane. As 
the aerials must be moved many times during the experiments this dimension was 
limited to about gm. At 210 Mc/s the reception pattern in galactic longitude 
covered about 12°; at 81-5 Mc/s it was about 30°. The largest spacing between 
the aerials was about 60 wave-lengths in both experiments. Additional observations 
on 81-5 Mc/s were made with a recording system capable of measuring the total 
aerial power (28), so that the zero-order Fourier component could be determined. 
No measurements of total aerial power were made on 210 Mc/s, and the results 
therefore only give the variation of intensity over a few degrees on either side of the 
galactic plane. 

Owing to the large spacing employed (up to 200 metres), it is important that 
the line joining the aerials shall be horizontal. | Observations were therefore made 
at times when the galactic plane was vertical; that is, when one of the galactic poles 
was rising and the other setting. ‘This situation occurs twice daily at all places 
where the galactic poles are not circumpolar; at Cambridge, the axis joining the 
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galactic poles is 40° west or east of north at these times.* The former of these axes 
was used as the interferometer axis, since observations can then be made nearer 
to the galactic centre; unfortunately the galactic centre itself (at galactic longitude 
327°) cannot be observed at Cambridge by this method. 

By using the aerials at various angles of elevation, the average distribution in 
galactic latitude was found for regions centred on galactic longitudes between 
338° and 16°. A few observations were also made near the galactic anticentre. 
Observations at intermediate angles were prevented by the confusion due to the 
intense radio stars in Cygnus and Cassiopeia. 

In the present preliminary investigations no attempt was made to determine the 
phase of the recorded pattern and it was assumed that the distribution was 
symmetrical about some plane near the assumed galactic plane. More refined 
experiments should prove valuable not only in investigating any asymmetry in the 
distribution of brightness temperature, but also in providing a method of deter- 
mining the position of the galactic poles from observations of objects showing high 
galactic concentration and situated in the most distant parts of the galaxy. 
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Fic. 4.—Calibration of the recording system. 

















(b) Calibration of the recording system.—A full account of the phase-switching 
system, including the use of pre-amplifiers at the aerials, has been given elsewhere 
(27). The deflection sensitivity C was found by connecting the output of a noise 
signal generator to the pre-amplifiers through the two equal aerial cables; the 
output of the signal generator then corresponds to the aerial power produced by a 
source of small angular diameter situated in the central maximum of the 
interference pattern. 

At 81-5 Mc/s, the total aerial power is comparable with that generated in the 
pre-amplifiers, and if, as is usual with the phase-switching system, the gain of the 
receiver is controlled automatically by the total output power from the 
pre-amplifiers, the gain will not be the same during the calibration. ‘The small 
error due to this effect was measured and corrected in the manner explained by 
Ryle (27). 


* Assuming the galactic poles to be at declination +28°. 
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Fic. 5.—Amplitude-spacing curves at 81-5 Mc/s. 
(a) Galactic longitude 1=353°. 
(6) i=13°. 
Ordinates: Amplitudes of record (arbitrary units). 
Abscissae: Spacing between aerials, n (wave-lengths). 


6. The experimental results—The graphs of amplitude against spacing shown 
in Figs. 5 and 6 contain the immediate results of the experiment. 
The distributions of brightness temperature deduced from these curves are 


shown in Figs. 7 and 8. It is considered unlikely that the overall calibration of the 
recording system and the calculation of the aerial reception pattern could intro- 
duce systematic errors greater than 20 per cent. 
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Fic. 6.—Amplitude-spacing curves at 210 Me/s. 

(a) 1=338°. 

(b) 1=358°. 

(c) 1=16°. 

Ordinates: Amplitude of record (arbitrary units). 

Abscissae: Spacing between aerials, n (wave-lengths). 
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Fic. 7.—Distribution of brightness temperature with galactic latitude, at 81-5 M[cs. 
(a) Galactic longitude l= 353°. 
(b) l=13°. 
Ordinates: Brightness temperature (deg. K). 
Abscissae: Galactic latitude. 
Experimental curves: 
Bolton and Westfold survey: 
Westerhout and Oort: —+-—+—- 
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(a) (6) (c) 
Fic. 8.—Distribution of brightness temperature with galactic latitude at 210 Mc/s. 
(a) 1=338°. 
(b) 1=358°. 
(c) [=16°. 
Ordinates: Brightness temperature scale (deg. K). 
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For comparison with the present observations, the corresponding curves from 
the survey of Bolton and Westfold (14) and the theoretical curves derived by 
Westerhout and Oort (9) for the contribution of radio stars and the extragalactic 
component are also shown in Fig. 7. 

The values of 7* were obtained from Fig. 7 by drawing a smooth curve having 
the same form as that derived by Westerhout and Oort, to extrapolate the experi- 
mental curves outside b = 3° to the galactic equator. 

At 210 Mc/s no measurements of total power were made, and no observations 
were made with spacings less than four wave-lengths; the distributions deduced 
are therefore only reliable within a few degrees of the galactic equator and are only 
shown up to b=2°. In order to derive the magnitudes of (7 —7*) at 210 Mc/s 
from Fig. 8, the angular width of the HII regions as observed at 81-5 Mc/s was 
used. 

Some additional observations at 210 Mc/s were made towards the anticentre ; 
these showed that (7'— 7'*) in that direction is probably less than 100 deg. K. 

(a) The angular distribution.—A comparison of Figs. 7 and 8 with curves 
given by Westerhout and Oort shows that the distribution in latitude is in good 
agreement with that predicted theoretically ; the width of the region found experi- 
mentally is, however, slightly greater than that predicted from the distribution of 
O- and B-type stars. 

Although observations have not yet been made at a sufficient number of galactic 
longitudes for a detailed comparison to be made, it is interesting to compare the 
observed variation of (7 — 7*) with galactic longitude with that deduced from the 
two models suggested by Westerhout and Oort. The first of these models assumed 
that the HII regions were contained within a ring whose inner and outer radii 
were 5 and 11 kiloparsecs ; the second assumed that they were distributed uniformly 
over a disk of radius 12-5 kiloparsecs. 

The expected variation with galactic longitude of (7—7*) at 210 Mc/s is 
shown for the two models in Fig. 9 together with the observed values. ‘The crosses 
show the observed values and indicate the estimated limits of error and the 
range of longitude covered by the aerials. It can be seen that the observed 
distribution shows an even greater central concentration than the disk model. 

The absence of an appreciable component in directions near the anti-centre 
also provides some evidence in favour of a rather restricted radial extent. 

(b) The kinetic temperature and effective density of the HII regions.—If the 
results are now used in the analysis described in Section 2, it is found that they lead 
toa kinetic temperature 7", of at least 8000 deg. K. (This figure was derived from 
a high estimate (1 400 deg. K) of the 210 Mc/s temperature, and a low estimate 
(4200 deg.K) of the 81-5Mc/s temperature.) The corresponding value of 
J[N2ds is 3 x 10% cm-. 

These results may be compared with the values used by Westerhout and Oort 
in their theoretical analysis: 


T= 10000 deg. K, | NZ ds =1-9 x 10% cm-®. 


An examination of the form of the curves shown in Fig. 3 shows that for 
kinetic temperatures as great as 18 000 deg. K, the present two frequencies do not 
offer the possibility of a very accurate determination of T,; better accuracy could, 
in this case, be obtained by making observations at lower frequencies. Never- 
theless it appears impossible to account for the present observations unless the 
kinetic temperature is considerably greater than the value previously assumed. 
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7. Discussion.—The experiments described in this paper were made in an 
attempt to find the radio-frequency effects of the H II regions. Sucheffects should, 
according to numerous investigations at optical wave-lengths, be confined to a 
narrow belt close to the galactic plane, and a narrow bright belt near the galactic 
plane has indeed been found in experiments at two different frequencies ; the width 
of this belt is only slightly greater than that predicted from observations at optical 
wave-lengths. 

There is, however, a marked difference between the distribution found in 
galactic longitude and that predicted theoretically. Observations of extragalactic 
nebulae at optical wave-lengths (29) have suggested that the H II regions are largely 
confined to the spiral arms; they have not been observed in the nuclear regions. 
The present observations, on the other hand, suggest that there is a marked 
concentration towards the centre of the galaxy. 








T 7 
330 360 30 60 
Fic. 9.—Comparison of observed and predicted distribution of (T—T*) in galactic longitude. 


(a) Ring model. 
(6) Disk model. 


Ordinates: Brightness temperature (deg. K). 
Abscissae: Galactic longitude (1). 


Estimates of the average electron density of the HII regions do not differ 
significantly from those adopted by Westerhout and Oort, but the deduced kinetic 
temperature ( <18 000 deg. K) is considerably greater; Spitzer and Savedoff (30) 
estimate that the temperature in an HII region should lie in the range 6000- 
13 000 deg. K. 

In view of this discrepancy, it is important to examine the possibility that some 
other source of radiation may be responsible for the observed bright belt along 
the galactic equator. 

It is possible that one of the two intense radio stars in the constellations of 
Cassiopeia and Cygnus might produce by means of a subsidiary maximum of the 
reception pattern an appreciable signal in the aerial which would be interpreted as a 
source of restricted width; both of these sources are however situated at consider- 
ably greater declinations than the region under investigation, and the records would 
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therefore readily be distinguished by their longer period. All of the radio stars 
whose positions lie within the limits of the bright belt have intensities far too small 
to make a significant contribution, individually. 

Recently, however, Mills (18) has found evidence for a concentration of intense 
radio stars near the galactic equator, and he has suggested that they may represent 
a population of “Class I” radio stars having a real concentration towards the 
galactic plane; it is therefore important to see whether the integrated radiation 
from such a population of radio stars could account for the observed bright belt. 
It is first of all clear that the spectrum of the radiation produced in this way would 
be the same as that of an individual radio star. The majority of radio stars have a 
spectrum which would give rise to a brightness temperature proportional toA*? (31); 
the radio star in Taurus has been reported (21) to have a brightness temperature 
proportional to A*®, but there is, as yet, no evidence in favour of a general law of this 
type for those radio stars situated near the galactic equator. 

The brightness temperature of the observed belt has been found to vary, 
between 210 Mc/s and 81-5 Mc/s, with an exponent /ess than 2-0, as would be 
consistent with an origin in the HII regions. Furthermore, some recent pre- 
liminary observations at 38 Mc/s have shown that between 81-5 and 38 Mc/s the 
excess brightness temperature decreases with wave-length. __ It does not therefore 
appear possible to account for the present observations in terms of the emission 
from a new population of radio stars, unless these are further supposed to have a 
spectrum very different from the known radio stars. 

It is possible that the observed belt might be accounted for partly by the 
emission from H II regions having a temperature of the order of 10 000 deg. K, and 
partly by the emission from a new population of radio stars; such a model allows 
considerable freedom as regards spectrum, and might also allow the necessary 
concentration in galactic latitude. However, if Mills’s value of the population 
density of the “‘ Class I’’ radio stars is used to estimate the intensity of the integrated 
radiation, it is found that it could only account for a few per cent of the observed 
belt. Furthermore, if the estimate of density is increased, the integrated radiation 
is found to vary only as In1/d, where d is the average separation between stars, 
whilst at the same time the degree of galactic concentration must be increased if the 
width of the belt is to remain sufficiently small. 

Whilst a more definite conclusion should be possible when more complete 
information is available on the distribution and spectra of a large number of radio 
stars, it is clear that with the present data there are very serious difficulties in 
explaining the observed belt in terms of radiation from radio stars. 

The most powerful evidence for an origin in the HII regions would be provided 
by an observation of the belt in absorption, but if the kinetic temperature is greater 
than 18000deg.K it will be necessary to make observations at frequencies 
considerably lower than 60 Mc/s. Some preliminary measurements at 38 Mc/s 
have shown that the belt is still just visible in emission; a kinetic temperature of 
about 25000deg. K would account for these observations and would also be 
consistent with the measurements at 81-5 Mc/s and 210 Mc/s. 

Some additional supporting evidence for the contribution of the HII regions 
to galactic radio emission has been provided by other observations: the constants 
derived from the present experiments appear to account for the differences in the 
galactic contours obtained at 100 Mc/s (14) and 200 Mc/s (20), and would give 
intensities at 1200Mc/s and 3000Mc/s comparable with those found by 
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Piddington and Minnett (21). Furthermore, as pointed out by Piddington and 
Minnett (32), the extended source “‘ Cygnus X’’, which lies near the direction of a 
spiral arm, has a spectrum similar to that expected from ionized hydrogen. Since 
the HII regions are found to be concentrated into spiral arms (29), all the available 
evidence suggests that the radiation from “‘Cygnus X”’ may be the integrated 
radiation from all the H II regions in the spiral arm. 


We should like to express our thanks to our colleagues in the Laboratory 
for assistance in making the observations, and to Mrs Moore for help in the 
computation of the results. 
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CINEMATOGRAPHY OF PARTIAL SOLAR ECLIPSES 
I. GENERAL 


R. d’E. Atkinson 


(Communicated by the Astronomer Royal) 
(Received 1952 November 6) 


Summary 


When a central eclipse of the Sun is viewed from a station slightly 
outside the central track, the thin unobscured crescent changes its position 
angle rapidly near mid-eclipse; a detailed record of its behaviour should 
provide good corrections to the Moon’s relative place in both coordinates. 
An expedition was sent to Mombasa for the eclipse of 1948 November 1, 
and obtained a successful timed cinema record. The present paper is an 
introduction to this work, and deals with general principles, choice of 
apparatus, choice of site, and the suitability of different eclipses. If only the 
P.A. of the line of cusps is measured, small images are adequate, and the 
camera can be stationary; if the light is reflected into it from a plane mirror, 
it can usually also be horizontal. Formulae are developed for the azimuth 
(and, if necessary, tilts) required to make the image move parallel to the 
long edge of the frame. The optical features of the Mombasa apparatus 
are outlined. Convenient approximate formulae are developed for the rate 
of rotation of the crescent, the sharpness of the cusps, and their divergence 
from each other, so that the suitability of a site can be rapidly assessed in 
any eclipse ; the maximum admissible sharpness and divergence are discussed, 
and eclipse statistics are summarized from this point of view. Both total and 
annular eclipses are suitable. 





If an observer be stationed near the central line of an eclipse of the Sun, but a 
little outside the actual track of totality (or of annularity), he will naturally see a 
thin crescent of unobscured solar surface even at mid-eclipse. It is easily verified 
that the position angle taken up by this crescent will change very rapidly through- 
out a period of several minutes ; if the minimum apparent separation of centres 
were 57”:3, a displacement of the Moon of 1”, perpendicular to the line of centres, 
would change the position angle of this line (and therefore of the whole crescent) 
by 1°, and, since the mean motion of the Moon is over half a second of arc per second 
of time, the P.A. would change by half a degree per second, at its fastest, if the 
Moon’s motion were fully operative. The apparent relative motion of the Moon 
past the Sun is somewhat less than the actual lunar motion, because the observer 
himself is being carried eastwards by the Earth’s rotation; moreover, a separation 
of centres as small as 1’ may in some cases produce a crescent too short, or too long, 
for satisfactory direct measurement of its P.A., or even prevent the appearance of 
any crescent at all; however, even allowing for these considerations one can in 
most cases count on a swing at the rate of 1° in 2 to 6 seconds at mid-eclipse. 

‘The question of using this rapid run to obtain a precise correction to the Moon’s 
place (relative to that of the Sun) attracted the writer’s attention late in 1947; an 
important advantage is that if an accurately timed run in P.A. can be obtained, it 
will give corrections not merely to one coordinate but to both. It is qualitatively 
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clear that the smaller the minimum separation of centres the more rapid the swing, 
and vice versa, so that not only does the coordinate parallel to the motion come out 
of the mean difference between predicted and observed position angles, but also the 
coordinate at right angles to the motion comes out of the rate of change; both are 
thus obtained from P.A. measurements only; there is no need to determine any 
auxiliary quantity, such as the apparent separation of centres or the distance between 
the cusps, by measurements of absolute lengths or positions. Indeed, it has been 
found in practice that when corrections for the lunar limb detail are used, a value 
for the difference between the apparent semi-diameters of Sun and Moon is also 
obtained, still without any need for measurements of lengths. 

The accuracy to be expected is roughly the same in both coordinates. If we 
take rectangular coordinates, with the axis of x parallel to the apparent relative 
motion of the Moon past the Sun at mid-eclipse, then, if the Sun’s centre is at the 
origin (S, Fig. 1), the Moon’s centre, L, will move along a track which, near mid- 
eclipse, is given by y =const. =xy say. 
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If the apparent relative velocity is v, and if the P.A. of the line SL is M, and 
that of the y-axis is M,, then assuming v constant we have 
tan (M— M,)=vt/Ko, (1) 


where ¢ is measured from the time of mid-eclipse. With roughly the values 
already assumed, i.e. v=0":32/sec say, =0'-32/min, and «x)=57":3, we have 


Rea! os : 
tan (M— M,) a; if ¢ is in minutes; also, 


—3438 vt dk (2) 


vt\2 * K ~~ 
1+ (2) pai ae 
Ko 


Thus if we put dko/ky = 335 and keep v/ky = 4, we see that the difference between the 

initial and final values of M changes 4’-8 for each 0”-1 change of x, if the extreme 

values of ¢ are + 1-5 min, and 5’:5 for each o”*1 if they are +2 min; these are 

increased (though not quite in proportion) if v is increased, and v may in fact often 
2# 


dM= 
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be somewhat larger than 0’-32/min, especially towards the ends of the track. 
‘Taking the figures as they stand, they represent the effect, on the extreme values of 
M, of a correction to the Moon’s place, in the direction perpendicular to its 
apparent relative motion, of 0”-1, and they are well comparable with the figure 
(0°-1 =6’) by which all position angles near mid-eclipse are changed if the Moon’s 
coordinate parallel to its motion is changed by 0”-1. 

The advantage of getting both coordinates with approximately equal accuracy, 
from observations at a single station, is an important one; moreover, if something 
of the order of 2 x 30° of the Moon’s limb is traversed by each of the two cusps, the 
Moon so located will in fact be very well defined, at least provided that the limb 
detail at the libration in question is well known. It thus appeared very desirable 
to develop a method of accurately recording the rapidly changing position angle as 
a function of the time ; timed cinematography was naturally the method demanding 
first study. 

Since there would clearly be no difficulty in timing each of several thousand 
cinema-pictures to about ;4; sec, the question of the attainable accuracy reduced, 
in the first instance, to the question how accurately position angles could be 
determined from the photographic images; these must evidently be small, if they 
are to be obtained on 35-mm film at all, and it might be thought that the accuracy 
of measurement would prove to be insufficient. A little study sufficed to show 
that such fears were probably unfounded. _ It is true that one might need elaborate 
precautions if one attempted to determine the actual positions of the centres of 
Sun and Moon by measurement of the available arcs of their limbs on a small- 
scale picture, and that a P.A. inferred from those positions might be of low 
accuracy; there is, however, an alternative possibility. The line of cusps is 
exactly at right angles to the line of centres for two circular disks, whatever their 
radii, and to a first approximation there is thus no need to locate the centres them- 
selves at all, or to know or determine the radii precisely ; and if in fact one confines 
oneself entirely to a direct measurement of the P.A. of the line joining the twocusps, 
the accuracy to be expected is fairly high. With this method, incidentally, a 
number of disturbing factors such as (especially) irradiation and atmospheric 
unsteadiness tend to cancel out, by symmetry, on the average. 

It is well known that the centres of photographic star-images can be located with 
an internal scatter of the order of half a micron or so, largely because the human eye 
is very sensitive to symmetry; a single cusp of a crescent is obviously a much less 
favourable object than a star-image, but the differential location of a pair of cusps 
brings in questions of symmetry again. If the photographic image be rotated, in 
the measuring-machine, and moved parallel to itself, until the tips of the two thin 
cusps just project across, or just stop short of, a single fixed fiducial thread in the 
field of the eyepiece, and if the cusps are really well defined, most experienced 
measurers would probably agree that a relative change of position, in the direction 
perpendicular to the thread, of 30 u should be clearly detectable even at the fairly 
low magnification which must be used if the whole crescent is to be visible at one 
viewing. If the distance between the cusps is 10 mm, a relative displacement of 
30 represents a rotation of about 10’ in P.A., and it is thus reasonable to expect the 
random scatter to remain below this figure provided the two cusps are cleanly 
defined. It is true that differences in their shape may introduce errors, and that 
such differences are in fact sure to arise on account of the rough contour of the 
Moon’s limb; but on the average these differences also may be expected to favour 
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one cusp as much as the other and so to produce errors of measurement which will 
still be quasi-random. (The actual displacement of the cusp, which is really 
produced by a lunar feature, may be larger than this error in its determination, 
and must be allowed for directly, so far as the limb-contour is known.) A truly 
random scatter, even if somewhat larger than 10’, would be of little importance 
if two or three thousand pictures were obtained, and it thus appeared probable that 
the small size of the images would be no real drawback, and that the method would 
be worth trying out at an actual eclipse; it was decided to plan for the eclipse of 
1948 November I, visible in East Africa. 

This plan was carried through, and a successful film covering nearly 3 minutes 
of time, with a total run in P.A. of about 76°, was obtained. The results will be 
discussed in detail in later papers, but it may be stated here that the above pre- 
liminary estimate of the attainable accuracy in P.A. has proved to be very conser- 
vative ; in favourable cases, which include the majority of the 3 100 images secured, 
the actual probable error, for re-measurement of one frame by one observer, has 
turned out to be approximately 1’. There are some systematic effects in addition, 
but these are in general only of about the same magnitude, and it may even be 
doubted whether increasing the scale (if that were in fact practicable) would reduce 
them very much. Such accuracy, corresponding as it does to a differential 
displacement of 3 1 between the two cusps, is naturally only possible if the images 
are sharply defined, and it was an essential feature of the work that this aim was in 
fact attained. 

The use of partial eclipses for the purpose of determining the Moon’s place is 
not in itself new. On several occasions it has happened that the track of a total 
eclipse passed fairly near to an established observatory, and programmes of 
positional work were then arranged for the existing equipment, especially helio- 
meters and large refractors; in one of the earliest cases, in fact, the observatory 
in question (KG6nigsberg) actually lay in the track of totality, and heliometer 
measurements of the P.A. and separation of the cusps gave corrections to both 
coordinates of the Moon’s place and also to its radius, with internal probable errors 
of two or three tenths of a second (1). Measurements with ordinary equatorial 
telescopes (2) have also been made using micrometer eyepieces, and even on the 
projected image; all methods employed relatively few discrete determinations 
(photographs, drawings, or direct measurements, as the case might be). 
A number of other quantities besides the position angle of the line of cusps have 
been studied, but most of these alternatives suffer from some disadvantages; in 
particular, irradiation can falsify the distance apart of the cusps, the P.A. of a single 
cusp with respect to the centre, the “‘ thickness”’ of the crescent, and (probably) the 
radii of curvature of the concave and convex arcs which define it (3). In so far as 
position angles are used, either alone or in combination, the effect on them of the 
adjustment and flexure of an equatorial mounting enters directly into the results, so 
that the direction of the diurnal motion must be observationally determined at the 
same time; sufficient accuracy in this may not be easy to get with a small number 
of images, especially if the atmosphere is unsteady. Some of the results obtained 
have, none the less, shown remarkably good internal consistency; but the work 
appears in general to have been regarded as a “‘ faute de mieux”’ enterprise, and the 
method was not developed for actual expeditionary use, except for some instances 
(4) involving the partial phases as observed from inside the track of totality or 
annularity; the really rapid part of the run in P.A. is of course lost in that case. 
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Cinema equipment has also been used at eclipses, one of the earliest cases being 
the “ limiting”’ eclipse of 1912 April 17, when a battery of cine cameras was set up 
across the track in the hope of determining just where the centre of the track lay (5). 
The records from several cameras showed Baily’s beads more or less well distributed 
round the Moon at mid-eclipse, and the “‘centre”’ of the track was taken to lie 
near that camera for which the distribution was the most uniform. This method 
in fact (in addition to being only very rarely practicable) sacrifices much of the 
advantage which a cinema record is in principle well able to give, since only a small 
number of pictures from any one station can actually be used at all, and, in addition, 
only a very few points on the lunar limb, namely the bottoms of the deepest valleys, 
had any weight in the answer. Even for these points, corrections based on the 
actual lunar limb detail were not applied; the Moon was in effect defined by a circle 
through the centres of gravity of the beads in question, or some approximation 
to this, and it might be difficult to relate this “‘ Moon” to any other one. 

Much the most important application of cinematography to eclipse measure- 
ments was made by K. Kordylewski, in a programme due to T. Banachiewicz 
(6), at the eclipse of 1927 June 29. Cine cameras were sct up inside the path of 
totality and were exposed for fairly short runs near second and third contacts. It 
was deliberately intended to make use of the Baily’s beads which naturally appear 
in such a case (particularly as the eclipse was only just total), and the times of appear- 
ance and disappearance of beads, and of gaps in short residual stretches of crescent, 
were computed from Hayn’s lunar charts and compared with observation. ‘The 
internal probable error of the Moon’s relative place so obtained was exceptionally 
small. ‘The work was not known to the writer until after the Mombasa expedition 
had returned; the method evidently differs from the present one in principle, and 
the Moon was also still located mainly by isolated features of its limb profile and 
not by a continuous survey. Given an accurate knowledge of the limb at the 
libration in question, the fact that this knowledge is unevenly used need not of 
course weaken the result much; moreover, although Kordylewski was himself of 
the opinion that the accuracy would be less if the solar and lunar semi-diameters 
differed more than they did in this somewhat exceptional eclipse, it is quite possible 
that the loss of accuracy would in fact be slight. ‘The method appears to deserve 
further trial (7, 8). 

There does not seem to be any case of cinema work deliberately undertaken 
outside the track of totality; but it is here that the advantages of cinematography 
should probably be greatest. It is only outside the track that observations can be 
made at the time when the position angle is changing most rapidly; and the cusps 
are considerably blunter than in Kordylewski’s method. Also, a small image is 
more acceptable for measurements of position angle than for the detailed study of 
beads, though the degree of detail obtained by Kordylewski was undeniably good. 
The very large number of pictures obtained in the present method, coupled with 
the unexpectedly high precision with which they can be measured, suggests that 
very detailed information about the Moon’s limb can be used, if it is available. 

In all photographic methods, there is evidently a great advantage in having an 
entirely stationary camera, with no coelostat. Quite apart from the increased 
simplicity and the greater rigidity which result, it is valuable to have the direction 
of the diurnal motion directly recorded, whenever position angles are required. 
In the present case (and probably in most others) there is then no need to rely on 
any angles measured in the field; all the necessary data are contained in the 
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film itself, and can be evaluated in the laboratory at leisure and checked in case of 
doubt. However, if a stationary camera is to be used, the focal length and the 
dimensions of the frame must be suitably related to the intended duration of the run. 
With a 35-mm camera, the frame measures about 23 mm x 17:5 mm; if the diurnal 
motion of the image is arranged to run parallel to the long side, and if the solar 
image is 9 mm in diameter, then, since the diurnal motion is a little slower (even on 
the equator) than one solar radius per minute, it is possible to keep the camera 
stationary and still have the entire solar image “‘ on”’ fora period of rather more than 
three minutes. This will fill a standard 200-ft. roll of film if the speed is about 
174 frames/sec. The duration can in some cases be extended; if only the tips of 
the crescent are actually required to be ‘“‘on”’, and if the line of cusps is nearly 
parallel to the short frame-edge at the start or finish of the run, even the centre of 
the Sun may sometimes be outside the frame, and a total run of about 4 minutes is 
quite possible. Moreover, if either the north limb of the Sun or the south one is 
unobscured throughout most of this time, one can expect to get a very reliable 
direct determination of the angle (ideally zero) between the long frame-edge 
and the direction of diurnal motion, since this direction appears as the common 
tangent to all the limb images obtained. (There is a slight curvature, if the 
declination is not zero, but this can easily be allowed for.) If position angles are 
measured, on individual frames, with respect to the direction of this same frame- 
edge, they can then at once be converted to true directions on the celestial sphere. 
Although the azimuth of the camera thus does not appear directly in the results, 
it is none the less desirable to set it up in such a position that the direction of diurnal 
motion will in fact be closely parallel to the frame-edge (or to any other fiducial 
line which may be used instead). This is partly to facilitate adjusting the image 
to give a “‘central”’ run, but principally because the procedure for determining the 
direction of diurnal motion, from the photographs, will be to measure the distance 
of the limb from the fiducial line, on as many frames as seem necessary, and to 
divide the change in this distance by the distance moved parallel to the fiducial 
line; this latter distance may however be somewhat difficult to determine 
accurately, since the E.-W. coordinate of the north (or south) limb is of course 
poorly defined. The total distance of E.-W. travel is only about 14mm on the 
assumptions used above, and, if it cannot be obtained to better than (say) I per cent, 
the angle cannot either; it would in that case be desirable for the angle to be so small 
that an error of I per cent in its value was unimportant. It is true that the E.—W. 
travel can be determined by other measurements (e.g. of the east or west limb, 
if present, or of the focal length), so that if the camera were in fact wrongly set up the 
results could still be used; there are also alternative measuring-procedures, at 
least if point-like fiducial marks are available; but these procedures all seem to 
involve more measurement. It is therefore preferable to get the direction of 
travel correct within at any rate half a degree, and there is no real difficulty in this. 
This does not, of course, imply that the camera must be set up in a tilted position ; 
the obvious procedure is to have the camera itself horizontal, and to mount a plane 
mirror in front of the lens; the image of the crescent may then be brought into 
place just before the start of the run by adjusting this mirror empirically, and it 
will in general be possible to select such an azimuth that at the time of the eclipse 
the image will travel horizontally. (This use of a mirror is clearly also the most 
convenient arrangement from the point of view of focus-tests, rehearsals, etc., 
since it is easier to adjust a mirror, so as to bring the image on to the front of the 
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camera, than to swing the whole assembly about for this purpose; the diurnal 
motion will not normally be horizontal in such rehearsals, but this is immaterial.) 
The required azimuth may be determined as follows. 

In Fig. 2, S is the Sun, H-H the horizon, Z the zenith, N the point where the 
mirror-normal meets the celestial sphere, C the point (in the horizon) where the 
optical axis of the camera (in the sense lens-film) does so; z is the Sun’s zenith 
distance and A’ is the difference in azimuth, which it is required to obtain, between 
the Sun and the camera-axis. The figure is drawn for the case that the observer 
is looking ‘‘out’’ at the celestial sphere, that the hour angle and latitude are 
positive, and that the Sun would be south of the zenith at noon, and S’S is the 
direction of diurnal motion, differing from true westwards by the angle Y (always 
<5’) which arises from the Sun’s changing declination. (Evidently Y is given by 
, secd. Ad 
tan Y= ——— 

15 At 
At; if A& is in seconds of arc, At is in seconds of time.) C’C is the direction of 
motion of the Sun’s image, and it is required that this shall be horizontal, as shown. 
QO is the parallactic angle; it appears simpler to use Q, s and A’ than to work in 
terms of hour angle, declination, latitude, and the absolute azimuth of C. 


, where Aé is the change in (geocentric) declination in the time 
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From the laws of reflection, it is clear that SNC and S’NC’ are great circles, 
and that the two spherical triangles NCC’ and NSS’ are equal in all respects; 
NS =NC =1is the angle of incidence (and reflection) at the mirror. If W is the 
auxiliary angle ZSC, then from consideration of the angles at S we have 


NSS’+7=W+O4+ = =, 


while the angles at C give 
NCC’ == -U 
2 


if U is the angle ZCS. 
The congruence of the two triangles thus gives 


am—-U=W+O-Y 
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or, writing O-Y=V, 

sin U=sin(W+V). 
From the triangle ZCS we have directly 

sin U=sin W sin z, 
sin V 

sin z—cos V~ 
From the same triangle we also have directly 


so that tan W= 


tan W= —tan A’ sec zg, 
and eliminating W we thus obtain 
tan A’ = See. hn ‘ (3) 
cos V —sinz 
The sign of A’ in this is positive as drawn in the figure, i.e. when A’ is positive, 
C is to the left of S; it will be found in fact that the meridian always lies between 
ZS and ZC. 

If the line of fiducial marks (or the edge of the frame) is not horizontal when the 
base of the camera is, or if a slanting run is desired for any other reason, the formula 
for tan A’ may still be used by writing V = Q— Y—T, where T is the slant of the 
line in question, reckoned positive if the image is to rise with advancing ¢. 

For the special case of the 1948 eclipse in East Africa the circumstances 
were unfavourable; the values of cos z and cos V (or cos Q) were both fairly 
small, and the value of A’ would thus have approached 180°. Since 
cos 2/=sin z cos A’, 27 must also have approached 180°, so that an undesirably 
large angle of incidence on the mirror would have been involved. The altitude 
of the Sun at mid-eclipse was in fact only about 22° even on the coast, and the 
direction of diurnal motion was almost vertical, so that one could hardly expect 
to produce horizontal motion of the image, in a horizontal camera, by one 
reflection at an acceptable angle. In principle, the solution might have been to 
lay the camera on its side, so as to get the long edge of the frame vertical 
(i.e. to put T=90°); but the cameras available did not appear to have been 
designed for use in that position and it was considered unsafe. There seemed 
no objection, however, to tilting the camera forward; the reflected ray would 
then run up at a fairly steep slant, and it was qualitatively obvious that this 
would ease the situation greatly. On repeating the above analysis, but with 
the zenith distance of the camera axis taken as w, instead of 7/2, formula (3) is 
replaced by 





cos V —sin 2 cosec w 
F cos zsin V 


sin A’ = cos A’—cot w tan s. (4) 
This equation is evidently less convenient than the other; it can, however, be 
solved for A’ fairly readily by successive approximations, or it may be solved 
directly by writing it in the form cos (A’ +a) =cot w tan z . cos a, where 

cos V —sin z cosec w 


tan a= - : 
cos zsin V 





in either case, a value of w is first selected more or less arbitrarily, and the two 
factors (cos V —sin z cosec w)/cos z sin V and cot wtan z are computed for this. 
The formula for cos 22 is readily amended, and, if the final result does not give a 
small enough value of 7, a fresh value of w may be tried; in practice, however, 
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it seems unlikely that a case less favourable in this respect than the 1948 one 
could arise, and the value then selected for w, namely 45°, is probably small 
enough for all circumstances. It led to practical arrangements which were in 
most respects quite convenient. 

Choice of apparatus.—If the solar image is to be about 9 mm in diameter, 
the focal length must be roughly 96cm. Cinema lenses of this focal length are 
not common, and would not in fact be very suitable; but since the ends of the 
23-mm field will be only about 40’ from the centre, a good quality small-transit 
lens (i.e. a visual telescope lens corrected for coma) will give adequate cover 
even though it is only a simple achromat. Such a lens will normally have an 
aperture of 7 or 8cm; it may be stopped down somewhat in order to reduce 
the intensity (or to improve the geometrical definition), but it is in fact undesirable 
to reduce the aperture much since the resolving power would then suffer; it was 
assumed above that the final 30 » of the tip of a cusp would be clearly noticeable 
in the measuring machine, and it would seem advisable to keep the diffraction disk 
below this tigure, so that the true appearance of the cusp, as caused by the actual 
lunar contour from moment to moment, may be reflected in the image as far as 
possible, and also so that the image may be as sharply bounded as possible. 
On the scale we have been considering, the radius of the first dark ring, in the 
image of a point-source, will be 30, if the diameter of the stop is 2 cm; it would 
clearly be well to keep the aperture larger than this. The /f-ratio will thus be 
wider than f/48; indeed, the full aperture of f/12 appears preferable if the lens 
is good; and since a photoheliograph can work at an effective f/250, and with a 
very slow emulsion, and still give an image at = sec exposure, it is plain that 
there is a very considerable amount of light to spare. 

A reducing factor of nearly 20 can be obtained by using the mirror unsilvered ; 
this has the additional advantage that there is then no risk of burning any filter 
or of burning the film if the shutter happens to stop in the “open’’ position. 
(It is of course necessary to eliminate the risk of reflection from the rear surface 
of the mirror; what one requires, therefore, is a “ solar flat’’, or alternatively one 
face of a good spectrograph prism, of suitable size.) A further reduction in 
intensity will be obtained if a filter is used, and since fairly narrow-band green 
filters are commercially available one is thus in the fortunate position of being 
able to eliminate the secondary spectrum of a visual achromat very completely 
indeed. With a lens that is corrected for coma, as well as spherical aberration, 
approximately down to the Rayleigh limit, the secondary spectrum is much the 
most important remaining factor in photographic image-defects, especially with 
fully exposed images such as are desirable in the present case, and the improve- 
ment in photographic definition which a narrow green filter gives is very marked. 
‘Trial apparatus assembled on these lines, using Kodak ‘‘ Background X” 
(panchromatic) film, a Wratten 57 filter mounted close to the camera, a lens 
from one of the small transits made by Troughton and Simms for the transit 
of Venus expeditions of 1874, and a “‘ Vinten K”’ (obsolete Air Force pattern) 
cine camera (with its own lens removed), gave excellent test-images in the early 
summer of 1948; sunspots less than 10” in size were very readily seen under the 
magnifying glass, and could be checked against the photoheliograph pictures of 
the same day. 

In these tests the telescope tube was cut off some distance from the camera, 
so as to allow the solar image to be easily seen on the front of the camera and 
adjusted by moving the mirror (prism); it was found that with the green filter 
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in place, and using exposures of about ;4; sec, there was very little fogging 
by direct light from the sky or the surroundings, even when the camera was left 
completely uncovered after adjusting the image. This was in fact regarded as 
unfortunate, since it had been intended to use weak fogging to show up the edge 
of the frame, so that it could serve as the fiducial line in the measurements; it 
now appeared almost certain that when the general level of illumination was 
reduced by a factor of 40 or so by the eclipse, the scattered light would be 
decidedly too faint for the purpose. Actual fiducial marks were thus desirable, 
but no practicable system for impressing these on the film could be devised in 
the time available, and it was decided to rely on the perforations, which are 
punched with extraordinarily high accuracy (as indeed they must be if the 
projected image is not to flicker). It was realized that under these conditions a 
scatter could arise, both in the position angles as referred to the holes, and also 
in the positions themselves (which are used for the zero of position angles), 
if the film ran in any way variably in the camera; but it was thought that such 
effects would be random, and also that they must already have been shown to be 
small by the absence of appreciable flicker. Even if these expectations were 
disappointed (as, in fact, they proved to be), it would not interfere with the 
judgment of the method’s suitability otherwise; fiducial marks could be 
incorporated in any later work, though no satisfactory method of actually doing 
this was evident in 1948. 

Chotce of station.—-The site of the observation station will influence the shape 
of the crescent, the sharpness of its cusps, and the rate at which it will rotate in 
position angle. These quantities may be evaluated, with ample accuracy for 
the preliminary Purpose of choosing a station, as follows. 

As is usual in eclipse work, we pick a fundamental plane perpendicular at 
all times to the axis of the shadow; for present purposes, it is more convenient 
to think of the plane through the observer, satisfying this condition, than of the 
parallel plane through the centre of the Earth; the numerical values of m, M, n, 
and N are not affected by this change. In the standard notation (m, ) is then 
the vector giving the distance and direction from the observer to the point where 
the shadow-axis cuts the plane, and (n, NV) is the vector giving the relative velocity 
of that point with respect to the observer; directions are reckoned from the 
projection of the Earth’s axis on to this plane, in the sense north-east. Mid- 
eclipse is evidently (as always) defined by M= N+ 7/2, m being then a minimum 
equal to my, say, and we compute Mp), my, No, my for this instant in the standard 
manner. Since, now, the unit of length for m and n is the Earth’s equatorial 
radius, the apparent separation of centres of Sun and Moon, « say, is 

k=(7q—7Q)m, (5) 
where 7@ is the Sun’s (geocentric) parallax, and zq is the Moon’s geocentric 
parallax increased to allow for the fact that the observer’s distance to the Moon 
is appreciably less than the geocentric distance; with sufficient accuracy we may 
write (Q) 

mq =7qC (1+cosz.sinz¢@), (6) 
where zq is the geocentric value and 2 the zenith distance ; for present purposes 
we may also put 70 and 7¢ constant, so that « becomes simply proportional to m. 

A little reflection makes it clear that M, which is defined in effect as a 
‘position angle ’’ on the fundamental plane, is also the position angle, on the 
celestial sphere, of the line apex-Sun—Moon as seen by the observer (looking at 
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the apex); or, since the convergence of meridians is entirely negligible as between 
apex and Sun, M is the position angle of the Moon, reckoning from the Sun. 
It follows, if « is proportional only to m, that (n, N) is parallel in direction and 
proportional in magnitude to the vector rate of change of the vector x, i.e. to the 
apparent relative motion of the Moon past the Sun; in fact 


v=(7(—70)(n, N) (7) 
is this velocity, in the same angular units in which z is measured. ‘The relative 
velocity (n, N) of the shadow-centre on the fundamental plane is very nearly 
constant; the main portion of it is the shadow’s own motion, which changes 
only very slightly indeed in the few minutes involved, and even the secondary 
part (the projection on to the plane of the reversed diurnal motion of the observer) 
does not vary much. Accordingly, we also have v=const., as was assumed 
in Fig. 1; and since the axes of that figure were already chosen parallel to the 
directions given by Ny and M, respectively, we may confine ourselves to the 
magnitudes of the vectors again, and write 


+y =(7(-—70—) mM =K,y =const., (3) 
v =(7(—70)mMp =Const., 


x=vt, 


tan(M— M,)=x/y= + 2 (0) 


k =k, sec(M— M,), (10) 


where ¢ is in minutes of time; x, y and « are in minutes of arc if 7¢ and 76 are. 
There is no need to go through the somewhat lengthy calculations of M, m, x, 
etc., for any other times; these equations give all the accuracy necessary and with 
very much less labour, though of course the rigorous ones must be used in the 
actual analysis, and will indeed naturally be employed as a check as soon as the 
station has been definitely selected. ‘The double sign in (8) and (g) arises from 
the fact that m is always taken as positive by definition; y will be negative, in 
Fig. 1, and tan(M— M,) will decrease with time, if the station is north of the 
centre-line; and conversely. (The opposite holds in those rare eclipses where 
the track crosses the Arctic or Antarctic “‘ beyond’”’ the pole.) 

The total angle swung through, 2(M— M,), is immediately determined by 
(9) as soon as the limiting values of ¢ are given; ¢ is limited by the focal length, 
and cannot in fact greatly exceed + 2™, whatever the focal length and field, without 
making a more elaborate lens necessary. ‘The exact limiting times will depend 
also on the P.A. of the cusps at the times in question, and on the Sun’s declination, 
but in any one eclipse they cannot vary much once the lens is chosen. Since, 
noW, Mp is insensitive to the observer’s exact position, so long as the general 
neighbourhood is fixed, while m, is markedly sensitive, equation (9) shows that 
in effect the limit to the total swing that can be obtained is set by m); from this 
point of view the smaller my, can be made the better. 

It is however clear that if « is too small, an unsatisfactory crescent may result. 
If v and o are the apparent semi-diameters of Moon and Sun, then in an eclipse 
which is total on the centre-line (i.e. v>c), there will be no crescent at all unless 
the station is chosen so that « >v—a, and for values only slightly above this limit 
the crescent will still be short, and with strongly diverging and very sharp cusps. 
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The question of their sharpness must first be settled; it is not entirely clear 
a priori that a separation of centres as small as 1’ might not always involve such 
sharp cusps that the irregularities of the Moon’s outline would make them behave 
rather wildly, showing incipient or even complete Baily’s beads for a large fraction 
of the total time. ‘The two tangents at a cusp being at right angles to the radii 
there, it is evident that if «=1', and v ¥o 106’, the cusp-angle will not be more 
than about 3$°, and will indeed fall below even this figure if the crescent is 
markedly less (or more) than a semi-circle; and the slopes of mountains on the 
Moon being notoriously much steeper than 3}° in general, it might be thought 
that such thin cusps would be quite unmanageable. In fact, however, this is 
not the case; large-scale photographs of the Moon, visual observations, and the 
traces published by Hayn, all agree that although there are some fairly steep 
places on the limb, there really are relatively few of them, and the number of 
cases where a slope greater than three or four degrees persists for any appreciable 
distance is not in itself large enough to cause serious anxiety. The reason for this 
relative smoothness of the limb is that at the limb, if the region is a mountainous 
one, we are looking at a jumble of mountains one behind the other, and if one 
mountain falls steeply, it will often merely unmask the next in the line of sight, 
so that the limb is little affected. It would be increasingly unsatisfactory to 
sharpen the cusps still further, since even when Baily’s beads did not actually 
form the cusps would become very sensitive even to moderate slopes on the 
limb; but it seemed probable that provided the minimum cusp-angle was not 
less than 3° one might expect serious difficulty only over short isolated stretches. 
For the first trial, in 1948, it was decided to accept this, and the results did in 
fact confirm expectation in this respect. 

In Fig. 3, S and L are the centres of Sun and Moon as before, and C is a cusp; 
the cusp-angle is @, and the half-lengths of the two circular arcs which limit the 
crescent are defined by the angles % and x at the centres of Moon and Sun 
respectively. Then evidently 


v+K2—o" 
2KV 
2 


cos yf = 


: (11) 


2 


, (12) 


v—K?—o 


cos x = ‘nae 


er 
and cnnta Santen , (13) 


2vo 


From the last of these, we see at once that if « is varied (either by moving the 
station, or by the course of the eclipse), cos @ always moves in the opposite sense ; 
since @ is always positive (and small), it follows that the smaller « is, the smaller 
6 must be also. Once the station is fixed, @ will therefore have a minimurna at 
mid-eclipse, and in estimating whether a station will permit a suitable value of 
@ it is sufficient to consider mid-eclipse only. For this purpose it is convenient 
to use sin@=(«/v)sinx instead of (13); it will appear below that x (which is, 
of course, equal to 4+ or to (¥%+ x)/2+ 6/2) should probably not be allowed to 
fall as low as 70°, so that sin x will exceed 0-95 and will moreover vary only slowly 
as x changes; the value of @ can in that case be estimated at once for any given «x, 
by writing 
sin 6 = 0-g6x/v. (14) 
y and yx will also have minima at mid-eclipse, but in the case of y% it is not 
immediately clear that the minimum will be the smallest value since ¢ has also a 
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maximum, namely when «x? =v?—o”, i.e. when x=7/2. In practice, however, 
it will be found that unless v and o are equal within a few seconds, this maximum 
will occur further from mid-eclipse than the run will extend, and the minimum 
value will usually be the smallest also; if values for some other time than mid- 
eclipse are nevertheless required, they can always be readily obtained by evaluating 
the appropriate « from (10), and inserting it in (11). 





Fic. 3. 


What we shall in fact be most interested in is not % but (¢% + x)/2, since this is the 
real measure of the divergence of the cusps. Ifthese are not parallel to one another, 
they wll not be perpendicular to the cross-hair in the measuring machine, and 
although we can (and, indeed, must) accept some obliquity in this respect it 
would probably be undesirable to allow a cusp to make an angle of less than 
about 70° with the cross-hair. Indeed, the risk of systematic error due to residual 
astigmatism in the observer’s eye is considerably more serious than is sometimes 
realized, and, if the cusps diverge (as they would in that case) by more than 40° 
from each other, even the use of several different positions for every measurement 
may not altogether eliminate bias. 

The “‘direction’’ of a cusp may be taken as perpendicular to the bisector of 
the angle 0, and this bisector evidently makes an angle (¢ + x)/2 with the line SL. 
The (total) divergence, D say, is therefore D=x7—(+ x), and we require a 
convenient expression for the behaviour of D as « is varied. 

Since sin x =(v/c) sin, we have cos(+ x) =cos¢%cos x — (v/a) sin? yy, and, as 
sin* % is readily available from (11), we obtain after a little elementary reduction 


cos (p+ x)= —1+2 (=) [:- ebay) 


4vo 


or cos? o+x _ (=) [: ae ). 
2 K 4vo 


Now 8va* 2 000, and x? will seldom be as large as 20; we may thus write 
g+x v-—o (16) 


~ 
~ 


2 K 


7. 
sin — =cos 
2 


with high accuracy. 
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Accordingly, the minimum allowable value of « necessary to keep the (total) 
divergence D below any specified figure can be readily estimated in any given 
eclipse (i.e. for any given value of y—c) from the following table :— 

TasLe | 


D 45° 40° 35° 30° 20° 15° 


25 
«/(v—o) 2°61 2°92 3°26 3°86 4°62 5°76 7°66 


Before evaluating « from this table it is necessary to increase the geocentric 
value of v by the same factor, (I +coszsin7q°), as was applied to 7@°. For an 
average eclipse, this will increase v by about I per cent; in an extreme case 
(perigee eclipse in the zenith) the augmentation can amount to 1°8 per cent, 
i.e. 0'-30. It is clear that the increase in vy—o may occasionally be as large as 
v—o itself, or even larger; there is no objection to this, but it underlines the 
importance of the correction. 

A rough survey of the frequency with which different values of y—o occur 
in total eclipses indicates that in about half the cases v—o <o’-4 before applying 
the above augmentation to v, or say <o’-6 after applying it; even for the worst 
of these cases, Table I shows that if D is allowed to reach 40°, x can be brought 
down to about 1-7, which corresponds to a total change of M of over 50° if ¢ 
can reach + 2™ and if v is 0’-4/min. Such cases may perhaps all be classed as 
“‘favourable’’, in contrast to those in which larger values of x) must be accepted. 
Geocentric values of v—o as large as 0’-85 are still fairly common; even if near- 
zenith cases are avoided (by picking stations towards the ends of the track, when 
necessary), « will then have to be nearly 3’, and M will therefore change rather 
slowly; the value of such eclipses, or alternatively the dangers involved in 
accepting a larger value of D than 40°, would probably have to be ascertained by 
trial. It might be found, in practice, that x, could be still further increased, since 
possibly the larger values of 6 which would then occur would more than com- 
pensate for the small value of )/m , by diminishing the effect of lunar limb- 
irregularities. 

If we combine the two conditions we have been considering, namely that 
sin shall not be less than 0-060 say (@= 33°) and that D shall not exceed 40°, 
we obtain 

0-060 = sin 8 = 0°95x«/v =0°95 x 2:92(v—a)/v, 

and putting v 16’ we have roughly 

_ 0:060 x 16’ 
0°95 X 2°92 
if v—a is less than this, the lower limit to « is set by the sharpness of the cusps, 
and if v—c is greater it is set by their divergence. In the Mombasa case, v—o 
was a little over 0’-2, which is considerably better than an average value; the 
maximum value of D was 30°, and the minimum value of @ was actually 3}°; 
this, coupled with the high value of v which the low altitude involved, led to a 
rate of rotation (31'/s at mid-eclipse), which is about the maximum obtainable, 
with the above restrictions. However, it is possible that a slower rate of rotation, 
with blunter cusps, would have been at least as satisfactory, so that eclipses in 
which v—o was larger would be very little inferior to the Mombasa one. (For 
all eclipses in which the rate of rotation is limited by 6, and not by D, it is 
insensitive to v—a.) 

The considerations in the case of an annular eclipse are in general very similar 
to those for a total one, and may readily be derived from them. Fig. 4 differs 


v-—o = 035; 
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from Fig. 3 only in having v<o instead of v>o, and the only changes which 
result in the conclusions are the following: 

If the eclipse is annular on the centre-line, o >v, and if a crescent (and not an 
annulus) is to be visible, the station must be chosen so that « >o—v; the cusps 
will then always converge, instead of diverging, but otherwise what has already 
been said about them stands. Equations (11), (12) and (13) are not altered; 
cos x is always negative, and x is still always greater than %; x and ¢% both approach 
= as « approaches o—v, and they both have maxima at mid-eclipse (though 4 
still has a minimum there); x also has minima (=7/2+ 9), if « is increased until 
as=7/2. yf should probably be kept below 110°, and @ can be readily estimated 
trom 


sin 6 = * sin ®0-96x o. (17) 





Fic. 4. 


The convergence of the cusps is — D=s+x—7 and the same formula holds 
as before for sinD/2; Table I is unaltered numerically, but the numbers now 
refer to —D and «/(o—v) respectively. The augmentation of v now therefore 
allows smaller values of « than the geocentric v would admit; in limiting cases, 
as is well known, the augmentation may change the sign of o—v and so produce 
an eclipse which is total in the central part of its run and annular only at the ends. 

The statistics of annular eclipses appear at first sight somewhat unfavourable 
by comparison with those of total ones, since the geocentric values of o—v not 
infrequently exceed 1'-2 and occasionally exceed even 1':5; however, if we 
consider also the fact that annular eclipses happen oftener, the worst cases lose 
significance. If we again allow an average augmentation of 0’-2, we find that for 
the best of the total and annular cases, respectively, the probability of getting an 
augmented value of v —o not exceeding 0’-6 is about } per year, and the probability 
of getting a diminished value of c—v in the same range approaches } per year. 
Exact figures depend on the actual augmentation in each case, and this is partly 
a matter of the smallest zenith distance reached by the eclipse, and partly a matter 





No. I, 1953 Cinematography of partial solar eclipses 33 


of how far from the ends of the track the station is placed. (The augmentation is 
naturally negligible at the ends of any track, whether short or long.) Annular 
eclipses do of course also involve a smaller relative velocity of the Moon, on the 
average, and also a smaller value of 7¢ ; so that on balance there will be nearly 
equal numbers of favourable cases in either class. The total number of total and 
annular eclipses is roughly double the number of “‘ favourable”’ ones in the above 
sense, and, as we have seen, it may well be found that even the ‘‘ unfavourable” 
ones can give good results, owing to the relatively blunt cusps which occur if a 
slow rate of rotation is accepted. ‘Thus an eclipse programme devoted to this 
method, if it can in fact use most annular eclipses as well as all total ones (so far, 
that is, as questions of accessibility and weather permit), should have many more 
opportunities of being carried out than those programmes which are restricted to 
total eclipses. Whenever the track runs over, or close past, two regions which are 
far apart on the Earth’s surface, the method will have geodetic applications, and 
the fact that two coordinates, and not merely one, are obtained will obviously be 
in its favour in this respect also. ‘The fact that sites may be picked on either side 
of the centre-line will sometimes greatly widen the practical possibilities, especially 
when the actual track just fails to touch land over an important part of its length. 
The personnel requirements are modest, since in principle one man can do all 
that is required at any one station. 

A detailed account of the actual work at Mombasa, and an analysis of the 
results obtained, will be given in subsequent papers. 


Royal Observatory, 
Greenwich : 
1952 November 4. 
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Summary 


As Einstein has pointed out, general relativity does not account satis- 
factorily for the inertial properties of matter, so that an adequate theory of 
inertia is still lacking. This paper describes a theory of gravitation which 
ascribes inertia to an inductive effect of distant matter. In the rest-frame 
of any body the gravitational field of the universe as a whole cancels the 
gravitational field of local matter, so that in this frame the body is “‘ free ’’. 
Thus in this theory inertial effects arise from the gravitational field of a 
moving universe. For simplicity, gravitational effects are calculated in 
flat space-time by means of Maxwell-type field equations, although a 
complete theory of inertia requires more complicated equations. 

This theory differs from general relativity principally in the following 
respects : 

(i) It enables the amount of matter in the universe to be estimated 
from a knowledge of the gravitational constant. 

(ii) The principle of equivalence is a consequence of the theory, not 
an initial axiom. 

(iii) It implies that gravitation must be attractive. 

The present theory is intended only as a model. A more complete, but 
necessarily more complicated theory will be described in another paper. 





1. Introduction.—In this paper we construct a tentative theory to account 
for the inertial properties of matter. These properties imply that at each point 
of space there exists a set of reference frames in which Newton’s laws of motion 
hold good—the so-called “inertial frames”. If other frames are used, 
Newton’s laws will no longer hold unless one introduces “fictitious” (inertial) 
forces which depend on the motion of these frames relative to an inertial frame. 

The question then arises: what determines the inertial frames? Newton 
asserted that they were determined by absolute space. However, absolute 
space is not observable in any other way, and it has been suggested that it is 
more satisfactory to attempt to correlate the inertial frames with observable 
features of the universe. In particular, Berkeley (1) and Mach (2) maintained 
that inertial frames are those which are unaccelerated relative to the “fixed 
stars’’, that is, relative to a suitably defined mean of all the matter in the universe. 
This statement is usually known as Mach’s principle. As this principle will 
be used as a guide in constructing our theory, we shall first discuss its general 
implications. 

The view that the problem of motion can be completely discussed in terms 
of observables implies that a kinematical description of all the relative motions 
in the universe completely specifies the system, so that kinematically equivalent 
motions must be dynamically equivalent. For instance, the statement that the 
Earth is rotating and the universe is at rest should lead to the same dynamical 
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consequences as the statement that the universe is rotating and the Earth is at 
rest, whereas this is not true in a scheme based on absolute space. Using Mach’s 
principle we can predict that the angular velocity of the Earth, as deduced from 
a local dynamical experiment (such as the motion of a Foucault pendulum), will 
be the same as that deduced kinematically from the apparent motion of the fixed 
stars. This prediction cannot be made in Newton’s theory, because there is 
then no causal connection between the motion of the stars and the existence of 
inertial forces at the Earth; the two observations give the same result only 
because, as it happens, the stars are not rotating relative to absolute space. 

If the rest of the universe determines the inertial frames, it follows that 
inertia is not an intrinsic property of matter, but arises as a result of the 
interaction of matter with the rest of the matter in the universe. This immediately 
raises the problem of how Newton’s laws of motion can be so accurate despite 
their complete lack of reference to the physical properties of the universe, such 
as the amount of matter it contains. It was largely this problem which originally 
prevented the general acceptance of Mach’s ideas, and one of the requirements 
of a theory of inertia that is consistent with Mach’s principle is that it should 
account for the apparent irrelevance of the properties of the universe. 

The observed fact that a gravitational force is locally indistinguishable 
from an inertial force, in that each induces the same acceleration in all bodies, 
suggested to Einstein that it is the gravitational influence of the whole universe 
which gives rise to inertia. General relativity was devised to incorporate this 
idea, but, as emphasized by Einstein (3, 4, cf. 5, p. 97), it failed to do so. Einstein 
showed that his field equations imply that a test-particle in an otherwise empty 
universe has inertial properties. In view of this it seems to be worth while 
searching for theories of gravitation which imply that matter has inertia only 
in the presence of other matter. In this paper we describe what appears to be 
the simplest possible theory of gravitation that has this property, though this 
theory is incomplete in other respects. 

2. General formalism.—Our problem is to construct a formalism in which 
the motion of a body is influenced by the presence of other bodies, but in which 
the concepts of “‘inertia’’ and “inertial frames’’ do not have to be introduced 
a priori. We shall represent the influence of the bodies on each other by a set 
of quantities defined at all points of space and time. As we are ignoring 
electromagnetic effects in this paper, we say that these quantities describe the 
gravitational field. The field is determined by the bodies (the sources) by means 
of a set of differential equations, together with suitably chosen boundary 
conditions. ‘These equations show how the field can be determined from the 
motion (and other properties) of the sources. 

In addition we must know how the field affects the motion of the sources. 
For this purpose we introduce the following postulate: in the rest-frame of any 
body the total gravitational field at the body arising from all the other matter in the 
universe is zero. In Newtonian language we could say that the universe moves 
relative to any body in such a way that the body never experiences a force—the 
difference from the ordinary Newtonian theory being that the forces acting on 
the body are here derived entirely from the matter in the universe. 

We must now set up the equations relating the gravitational field to its sources. 
In the rest-frame of any particle we assume the field to be derivable from a 
potential in Minkowski space, that is, we do not describe it in terms of a curved 
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space. Kinematical considerations (Section 6) show that the potential should 
be a tensor of the second rank, but the use of such a potential leads to rather 
involved mathematics which tends to obscure the physical significance of the 
theory. It seems advisable to begin by working with the simplest mathematical 
scheme which contains the physically important aspects of the problem. The 
simplest type of potential we could use is a scalar, but as we shall see (Section 4 (ii)), 
this would not give rise to inertia. ‘The next simplest possibility is a vector 
potential, and with a theory based on such a potential we can reproduce the main 
properties of inertia. In this paper we shall confine our attention to such a 
potential. This simplification is useful for gaining insight into the problem, 
but naturally it has its limitations, some of which are mentioned in Section 6. 
The more elaborate equations that are needed for a tensor potential will be 
described in a subsequent paper (hereafter called II). 

In a theory based on a vector potential, the field is an antisymmetrical tensor 
—the curl of the potential. The only linear second-order differential tensor 
equations for a field of this type that imply the conservation of source are (6) 
Maxwell’s equations, which accordingly we shall adopt. We emphasize that 
although our equations have the same formal structure as Maxwell’s, they 
describe purely gravitational effects, electromagnetic phenomena being outside 
the scope of this paper. 

In order to apply the theory to even as simple a problem as the motion of a 
particle in the gravitational field of the Earth, we must know the distribution 
of matter in the universe. In practice we shall have to approximate to this 
distribution in some way. ‘The type of approximation that will be most useful 
depends on the relative importance of near and distant matter. Since the amount 
of matter at a given distance increases roughly with the square of the distance, 
it follows that if the influence of matter falls off more slowly than the inverse 
square of the distance, then very distant matter is of predominant importance. 
It is convenient to anticipate that this is indeed the case (Section 4(iii)). This 
means that a smoothed-out model of the universe should be a good first 
approximation, local irregularities having only a small effect which can easily 
be estimated. It also means that local phenomena are strongly coupled to the 
universe as a whole, not just to local conditions. ‘This in turn means that local 
experiments, if interpreted by means of this theory, can give us information 
about the structure of the universe as a whole (cf. 5, 7). ‘The correctness of this 
information can, in principle, be tested by independent and more direct 
considerations. 

We shall take as our smoothed-out model a homogeneous and isotropic 
distribution of matter of density p expanding (relative to any point as origin) 
according to the Hubble law v = r/r, where v is the velocity of matter at distance r, 
and 7 is a constant. This neglects certain relativistic difficulties such as the 
significance of velocities exceeding that of light, but for the tentative theory 
developed in this paper we shall not concern ourselves with these problems; 
a consistent treatment will be given in II. ‘This model is one of those in which 
there is a natural state of rest at each point, namely, that in which the observed 
distribution of the red-shifts of distant matter is isotropic. Thus we can speak 
of a body being at rest relative to the universe. 

3. Inertia-induction.—In order to show how inertia arises in this formalism, 
we consider the behaviour of a test-particle in the presence of a single body 
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superposed on the smoothed-out universe. The problem is to determine what 
motion of the system universe-plus-body relative to the test-particle makes the 
total gravitational field at the particle zero. 

It is convenient to begin by calculating the potential at a test-particle that is 
at rest in a universe containing no irregularities. Since our field equations have 
the same form as Maxwell’s, we can use electrodynamic formulae to calculate 
the potential, and to bring out the analogy with electrodynamics we use a similar 
notation and terminology, but we emphasize that in this paper we shall be 
concerned with purely gravitational phenomena. 

Retardation effects are taken to arise in the same way as in electrodynamics, 
so that the contribution of any region of the universe to the potential at a 
point P at time ¢ is computed by ascribing to that region just the properties that 
are observed at P at time ¢. 

We thus have for the scalar potential (8) 


Say fC 
© = far. (1) 


We use the minus sign in (1) because inertial mass then turns out to be positive, 
but in fact either sign can be used (Section 4(vii)). The vector potential A 
vanishes by symmetry. 

We shall assume that matter receding with velocity greater than that-of light 
makes no contribution to the potential, so that the integral in (I) is taken over 
the spherical volume of radius cr. An assumption of this sort is necessary since 
we have naively extrapolated the Hubble law without considering relativistic 
effects, and should give the correct order of magnitude. A relativistic treatment 
is given in II. 

Since the density is supposed uniform, (1) gives 


= — 2mpc*7?. 


Owing to our assumptions, the numerical factor 27 is only approximate. 

We now calculate the potentials for the simple case when the particle moves 
relative to the smoothed-out universe with the small rectilinear velocity — v(t). 
In the rest-frame of the particle the universe moves rectilinearly with velocity v(t). 
Now at time ¢ there will be observable at the particle, in addition to the Hubble 
effect, a Doppler shift corresponding to v(t) from all parts of the universe. 
Hence, in computing the potential in the rest-frame of the particle at time ¢, 
we must ascribe to every region of the universe the velocity that is observed at 
time ?¢, that is, v(t)+r/r. 

Neglecting terms of order v?/c?, we have 


@ = — 2mpc?7? 


as before. ‘The vector potential no longer vanishes, but has the value 


vp... 
ine | ar. (3) 


Since v is independent of 7, we can take it outside the integral. We then 
obtain 


A= * v(t). 
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Since the change of p with time is very small, the gravelectric part of the 
field is approximately 
1 0A 
E = — grad D— oj 
® ov 
~ at’ 


while the gravomagnetic field is 
H =curlA=o. 

So far we have been concerned with a universe that has no irregularities. 
We now suppose that a body of gravitational mass M is superposed on this 
universe and is at rest relative to it. The field of this body in the rest-frame of 
the test-particle is then 

M,. ¢0v 
“27° ae (4) 
where r is the distance of the body from the test-particle, (= — M/r) is the 
potential of the body at the test-particle, and 
5 dv dv 
‘dt dt® 
The total field at the particle is zero if 


—_ (*3*)5 (5) 


re dt” 

This equation asserts that the system universe-plus-body accelerates relative 
to the test-particle at a rate determined by the mass and distance of the body. 
In accordance with the discussion of Section I, we may re-interpret this result 
by saying that the particle accelerates towards the body, which is at rest relative 
to the universe, and thereby make contact with the Newtonian view-point. 
Indeed we obtain for rectilinear motion a combination of Newton’s laws of motion 
and of gravitation, with the inertial frame determined by Mach’s principle. 

Furthermore, the gravitational constant satisfies the equation 


O+¢ I 


ce G 


c 


or, since 6 <® (cf. Section 4 (iii)), 
GO = —c’*. (6) 
4. Consequences of the theory.—(i) Equation (6) implies that the total energy 
{inertial plus gravitational) of a particle at rest in the universe is zero. It can in 
principle be tested observationally, for, in conjunction with (2), it implies that 
2nGpr?=1 
of 
Gpr? ~1I, (7) 
since the numerical factor is only approximate. 
Equation (6) implies that the gravitational constant at any point is determined 
by the total gravitational potential at that point, and so by the distribution of 
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matter in the universe. This illustrates the fact mentioned above that if local 

phenomena are strongly coupled to the universe as a whole, then local observations 

can give us information about the universe as a whole. With our assumption 

about the structure of the universe, a laboratory determination of G, combined 

with an astronomical determination of 7, enables us to deduce the mean density 

of matter in space (apart from the uncertainty in the numerical constant). 
Taking G=6 x 107-8 c.g.s. units and t=6 x 101® sec we have from (7) 


p~5 x10?" gcm"3, (8) 


A mean density of this amount is much larger than the usual observational 
estimates (~ 10° *° gcm~), but these values refer only to matter condensed into 
nebulae, and the mean density of internebular material may well exceed the 
overall density of material condensed into nebulae (5, p. 45). Indeed, if the 
nebulae are supposed to have condensed from an internebular medium, there is 
no reason why the medium should have been exhausted in the process. It is 
probable that p has an upper limit of ~10-% gcm~, as this appears to be the 
mean density within some nebulae. Hence (8) is not inconsistent with observation. 

(ii) We see from the argument leading to (5) that “‘inertia-induction”’ arises 
from the term 0A/ét, that is, from the “‘ radiation-field’”’ of the universe. Had 
we used a scalar potential such a term would not have arisen, there would have 
been no inertia-induction, and we would not have been able to obtain Newton’s 
law of motion (5). 

(iii) The contribution of matter to local inertia falls off only inversely as 
the distance, since 0A/dt is proportional to the scalar potential (cf. the radiation 
field of an accelerating charge (8, p. 22)). This means that the main contribution 
comes from distant matter—(1) shows that 99. per cent of local inertia arises 
from matter further away than 10° light-years. The fractional contributions 
of the Earth, Sun and Milky Way are 10-*, 10-8, 10-7 respectively (taking the 
mass of the Milky Way as 10“ g); this justifies the neglect of ¢ in (6). According 
to our theory, then, local phenomena are strongly coupled to the universe as a 
whole, but owing to the small effect of local irregularities this coupling is practically 
constant over the distances and times available to observation. Because of this 
constancy, local phenomena appear to be isolated from the rest of the universe. 

(iv) As has already been pointed out, a theory of inertia that is consistent 
with Mach’s principle must account for the fact that Newton’s methods were 
so successful despite their complete lack of reference to the properties of the 
universe. The work of Section 3 shows that our theory satisfies this requirement. 
The universe affects local phenomena at just the two points where Newton’s 
work contains arbitrary elements, namely (a) in the choice of inertial frames, 
and (6) in the value of the gravitational constant. 

(v) Relation (7) is a consequence of many general relativistic models (5, 
cf. also g) as well as of our theory. However, our theory is disproved if (6) 
disagrees with observation, whereas general relativity is not, since it leads to 
many models in which (6) is false. Hence in general relativity the observed 
value of the gravitational constant gives no information about the amount of 
matter in the universe. For instance, if the Milky Way was thought to be the sole 
occupant of the universe, general relativity would still be consistent with the 
observed value of G, whereas in our theory there would be a discrepancy of 
a factor ~107. Indeed, if we accepted Mach’s principle for the reasons given 
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in Section 1, we would, on the basis of our theory, predict that the universe 
contained vastly more matter than had yet been observed. 

(vi) The principle of equivalence asserts that the phenomena occurring in 
a gravitational field are the same as those occurring in the absence of the field 
if observations are made from a suitably chosen accelerating frame. For 
instance, consider an observer confined to a closed laboratory. According to the 
principle of equivalence the observer cannot distinguish between the following 
two events: 

(a) a gravitating mass is suddenly placed near the laboratory ; 

(6) the laboratory is suddenly accelerated, e.g. by being pulled by a rope. 
In our theory this is explained by the fact that in case (b) the motion of the universe 
relative to the observer produces the same gravitational field at the laboratory as 
that of the mass in case (a) (cf. the argument leading to (5)). 

General relativity has difficulty in explaining the principle of equivalence 
because it predicts that one gravitating mass in an otherwise empty universe 
produces effects similar to those calculated with a full universe in Section 3. 
Since there is no universe to act as the source of the field in case (4), it is difficult 
to see why the principle of equivalence should be true. 

(vii) It follows from our equations that gravitation is attractive, whereas 
in general relativity the sign of the field is not determined. ‘This result arises 
because our equations imply that the field at a point P due to a mass is decreased 
in absolute magnitude if the mass has a component of acceleration towards P, 
and increased if it has a component away from it, whichever the sign of the field 
(cf. (4)). Thus the field of the universe can only cancel the field of a local mass 
if the universe moves so that the mass, which is at rest relative to it, accelerates 
towards the particle. ‘This is why the sign in equation (I) is arbitrary: had we 
chosen a plus sign, gravitation would have been repulsive and inertial mass 
negative (corresponding to changing the sign of both sides of (5)), so that the 
particle would still accelerate towards the body. 

5. Uniform rotation.—In Section 3 we assumed that the universe-plus-body 
moved in a straight line relative to the test-particle. We shall now derive another 
possible motion of the system in which the universe and body rotate with constant 
angular velocity about an axis through the centre of the body perpendicular to 
the line joining it to the particle. 

In the rest-frame of the particle we set up a Cartesian system of axes with 
origin at the body and z-axis along the axis of rotation. If the universe were 
non-rotating, its gravitational potential near the origin would be 


A,=A,=A,=0; ®=-1, (9) 
choosing units so that G=c=1. If it rotates, its potential in the x, y plane near 
the origin is 

A,=wy, 
A 


A,=0; 


y= wrx, 


@ = —|1+w?(x?+4")]!? 


= —(1+wr?)t2, 


This follows from the transformation properties of four-vectors (10). It is only 
true for distances from the origin for which the potential of a non-rotating 
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universe differs negligibly from (9), since Rosen’s result was derived for a 
four-velocity which has the values (9) everywhere in the non-rotating frame. 
This restriction is of no importance for phenomena on a galactic scale since the 
potential of the universe at the centre of the galaxy differs from that at the edge 
by less than one part in 10". 
We thus get 
0A 


E = — grad ®— 73 


wr 
~ (t+ 0% 
~w'r, for wr<t. 
The field of the body, neglecting its rotation, is 
M 


3 r. 
The total field is zero if 


> =w*r. (10) 

In the rest-frame of the universe, this is the usual Newtonian equation for 
circular motion. In the rest-frame of the particle, however, Newton’s law of 
motion only holds if we introduce a fictitious centrifugal force-field. Equation (10) 
shows that in our theory this field is not fictitious, but arises from the gravitational 
effect of a rotating universe, in agreement with Mach’s principle. 

In contrast to the case of Section 3, the gravomagnetic field is not zero (cf. 
the magnetic field of a rotating charge-distribution). In fact, we have 


H =curlA=2w. 


Since the test-particle is at rest, this field has no effect; but for another test- 
particle constrained to move relative to the first with uniform velocity v, there 
would be in its rest-frame a gravelectric field 


vaH =2vaw (11) 


acting on it. This can be seen directly by taking this particle to be at rest and 
determining E in the new frame. (Cf. the relativistic deduction of the Lorentz 
force e(E+vaH) from the expression eE.) 

The field (11) just corresponds to the Coriolis field of Newtonian theory, 
and like the centrifugal field it is here ascribed to the rotating universe. Thus, 
in our theory, we can regard the Earth as stationary and a Foucault pendulum 
as pulled around by the gravomagnetic field of the rotating universe. 

Another feature of this interpretation of the Coriolis field is that the bending 
of light in a gravitational field is observable in the laboratory. In fact, the 
experiments of Sagnac (11) and Michelson and Gale (12) showed that in a 
rotating frame of reference light does not travel in straight lines. If we interpret 
this phenomenon according to our theory, we may suppose the frame to be at 
rest and the universe to be rotating around it. The resulting gravitational field 
then bends the path of the light. According to this point of view the Sagnac 
effect is of the same nature as the refraction of light passing through the 
gravitational field of the Sun. A quantitative treatment of the motion of light 
in a gravitational field requires, however, the more elaborate equations of II. 
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The results of this section enable us to reconsider the problem of absolute 
rotation from the standpoint of our tentative theory. In fact, rotation does not 
give rise to any problems not already raised by translation, but the conflict 
between absolute and relative motion has usually been discussed in terms of 
rotation. The main argument for the view that there is an absolute standard 
of rotation has been the fact that local experiments can detect rotation— 
experiments in which the rest of the universe has not appeared to play an 
essential part. The present theory answers this argument by showing how 
the rest of the universe can play an essential part in local phenomena. The 
results of local experiments can be interpreted as giving us information about 
the universe as a whole—an interpretation which is verifiable by direct observation. 

6. Limitations of the theory.—We have so far derived Newton’s law of motion 
for two special cases: rectilinear motion and uniform rotation. A satisfactory 
theory of inertia must, of course, derive Newton’s law for all possible motions. 
Now it follows from the relativistic form of Newton’s law (13) that for general 
motions inertial forces are derived from a tensor potential and not from a vector 
potential. This means that the gravitational field of the moving universe, which 
in our theory gives rise to inertial forces, must also be derived from a tensor 
potential, so that our present theory is incomplete. 

Furthermore, in a theory based on a vector potential it is difficult to give a 
consistent relativistic discussion of the structure of the universe as a whole. It is 
also difficult to describe the motion of light in a gravitational field. In paper II 
we shall describe a theory, using a tensor potential, in which these difficulties do 
not arise. 


I am grateful to my colleagues for commenting on the manuscript, to the 
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Summary 


The motion of a semi-infinite gas cloud bordering upon a vacuum is 
studied by the method of characteristics. Problems previously discussed 
by Burgers and Copson are illustrated by diagrams; they involve gas flows 
which are always continuous. It is shown that unless certain conditions are 
satisfied by the initial state in the cloud the subsequent flow cannot remain 
indefinitely continuous; shock-waves must appear within the gas. An 
example of the breakdown of continuity is given. A singular solution is 
considered in which all the characteristics of one family are concurrent in a 
point of the (x, t) plane. Finally, it is shown that the high eventual velocity 
of these gas clouds, which was considered to be a physically unsatisfactory 
feature of the previous solutions, is not to be avoided by variation of the 
initial inhomogeneities. 





1. Introduction—The equations which govern the flow of an ideal gas, 
neglecting the effects of viscosity and heat-conduction, have recently been used for 
the study of the motion of interstellar clouds. These clouds consist almost 
exclusively of hydrogen, and the value of y, the ratio of specific heats for this gas, 
may be taken as %. For this particular value of y an analytical solution of the 
equations of motion is easily obtained. 

Burgers* considered the case of a semi-infinite homogeneous gas cloud which 
is at rest at some instant, its front adjoining a vacuum; he gives formulae for the 
characteristics of the motion which instantly begins. These are the characteristics 
of a “‘simple wave’’ (see below). The face of the gas moves forward into the vacuum 
with a velocity 3a, where a is the speed of sound in the gas cloud in its initial state. 
At the same time a wave moves back into the cloud with velocity a, this wave 
separating the part of the cloud which has begun to move from the part which 
remains at rest. 

McVittief took up Burgers’ work and tried to extend it to less artificial initial 
circumstances, starting with a stationary non-homogeneous gas cloud. On the 
surface of the vacuum the velocity of sound c is zero, so that there is no initial 
discontinuity of pressure, and the value of c rises continuously towards the interior 
of the cloud. The discontinuity of initial velocity of the gas cloud in Burgers’ 
work is removed by this alteration; unfortunately, however, according to the 
physical model which McVittie adopted, the whole of the inhomogeneous cloud is 
never moving at-once, whereas such a cloud must begin to move immediately 
throughout its whole mass under any likely equation of state. 


* J. M. Burgers, K. Ned. Akad. v. Wet., 49, 588-607, 1946. 
+t G. C. McVittie, M.N., 110, 224-237, 1950. 
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Copson* applied the theory of Riemann invariants to the further example of a 
limited region of inhomogeneity backed by a region of homogeneous gas. In the 
limited region the speed of sound increased continuously from zero at the vacuum 
surface to a, the value in the homogeneous gas, at the inner boundary. 

In the following sections the solutions discussed above are illustrated by the 
method of characteristics, and the same method is used to discuss certain variations 
of the initial conditions which lead to new types of flow within the gas clouds. 
Conditions are given which are necessary if the flow is to remain continuous. 

The method of characteristics is a general one, applicable to any value of the 
ratio of specific heats and adaptable to any proper formulation of the initial value 
problem. With knowledge as to whether a solution will or will not remain 
continuous, it can yield, in diagram form and without detailed mathematical 
analysis, the essential features of the flow; otherwise, or to obtain quantitative 
information, a computation is necessary. In all cases the characteristics diagram 
enables the nature of a solution to be seen at a glance. 

It is shown that the eventual high velocity of advance of gas clouds bordering 
upon a vacuum, which was considered an unsatisfactory feature of the previous 
solutions, is not to be avoided simply by varying the initial inhomogeneities. 

2. The characteristics.—Let u, p, p be respectively the velocity, pressure and 
density in a one-dimensional unsteady motion, parallel to the x-axis, of ideal gas 
expanding according to the adiabatic law p=Kp’, K being a constant. ‘The 
equations of motion and continuity are respectively 

Ou Ou 10p 
nt ~ a = =- ? ps > (1) 
Idp udp du 
pot + bax * dx" x8) 
The velocity of sound in the gas is c where c?=dp/dp=yp/p. Writing (1) and (2) 
in terms of u and c, the equations become 
Ou du 2c dc 
th +5 
2 0c 2u dc du 2 
y—iat ' y—tdx ‘dx (2a) 
The equations (1 a) and (2a) may be replaced, after simple manipulation, by the 


] p 


i +(u- 2 (2 (4) 


(i) along a curve in the (x, t) plane for which = =u+c, the expression 


2c ; 
: +u is constant, and (54) 


(1a) 


which shows that 


dx ° 
— =u—c, the expression 


dt 


(ii) along a curve in the (x, ¢) plane for which 


2c ’ 
2s = ——— —u is constant. (55) 
Y a 





* E. T. Copson, M.N., 110, 238-246, 1950. 
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These curves are the characteristics of the differential equations governing the 
motion, and the conditions r =constant, s=constant, holding respectively along 
the two families of curves are called the conditions of compatibility; r ands are also 
known as Riemann invariants. For hydrogen, with y = 2, the conditions are :— 


x 
along a ~Bth 36+ = 27 =constant, (6a) 


x 
alon = =u—Cc, 3c —u=2s=constant. 6b 
6 dt 


The two families of curves will be referred to below as the r-characteristics and 
s-characteristics respectively. Through any point P in the (x, ¢) plane at which 
there is gas there will be two characteristics, one of each family. The characteristics 
represent the limits of influence of an infinitesimal disturbance originating at P. 
Given values of u, c at a point (x, 0), the values of r, s corresponding to the character- 
istics through this point on the line t = 0 are known. 

At the intersection of an r-characteristic with an s-characteristic the state of the 
gas is given by the simultaneous solution of the conditions of compatibility in 
(5a, b), viz. :— 


u=r—s, c=(y—I)(r+s)/2. (7) 

If there is a range of x for which the initial conditions are homogeneous, all the 
r-characteristics stemming from this region have the same value of r, and similarly 
all the s-characteristics have the same value of s. The equations (5) therefore 
show that if an r-characteristic (say) from a non-homogeneous region meets 
successively the s-characteristics from a homogeneous region, the values of u, c 
along it will be unvarying and hence such a characteristic will be straight. 


A group of such characteristics constitutes what is termed a simple wave. 
In general, both wu and c are varying along a particular characteristic. 

If it should happen that two characteristics of the same family intersect one 
another, the continuity of the motion breaks down. The values of u and ¢ are no 
longer uniquely determined by the conditions of compatibility. A shock-wave 
begins at the point and time when such an intersection first occurs, and moves in a 
determinable way through the gas.* 

3. Burgers’ solution.—Burgers considers a semi-infinite cloud bordering upon a 
vacuum. ‘The gas is supposed to be initially at rest and to be homogeneous. 
Thus, at ¢=0, u=0 andc =a for x <0, if x =0 is the initial front of the cloud and the 
positive x-axis is directed into the vacuum. Although there is here an initial 
discontinuity in density at x =0, the discontinuity exists only instantaneously, all 
values of pressure and density down to the vacuum state being immediately 
propagated with the appropriate velocity in a simple wave. The solution is 
illustrated in a diagram (Fig. 1). 

A section by t= constant contains homogeneous gas at rest, then a section of a 
simple wave leading to the vacuum. The wave of expansion is represented by a 
fan of characteristics centred on O. The s-characteristics through the points 
t=0, x <0 are straight with gradient dx/dt=—a. The r-characteristics through 
these points are straight until they meet the s-characteristics of the expansion fan. 
They then curve as shown. The characteristics of the fan itself meet only 
r-characteristics from an initially homogeneous region; they are therefore 
straight and constitute a simple wave. At the front of the gas c=o; hence 

* A. F. Pillow, Proc. Camb. Phil. Soc., 45, 558-586, 1949. 
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u=2r=2a/(y—1). Thus, the particle velocity in the surface of the cloud is 
constant and equal to 2a/(y—1) (i.e., 3a for y=%). This is also the gradient of the 
leading characteristic. The rarefaction spreads backwards as well as forwards 
and a glance at Fig. 1 shows that the boundary of moving gas itself moves into the 
undisturbed gas with a velocity 

dx - I 

= =u—c= ttre Wet s=(1-y)r9= —a 
(where ry) =5 )=the initial value of r, s on this boundary), i.e., with the speed of 
sound in the homogeneous cloud. All the r-characteristics which initially set out 
from the region x <0 cut a particular s-characteristic of the fan at the same angle ; 
they eventually touch the leading characteristic at infinity. Thus, small signals 
from the interior of the cloud, which would be propagated along the characteristics 
in the (x, t) plane, never quite reach the front of the cloud. 


’ 


4 





REGION OF HOMOGENEOUS \ if 
GAS AT REST, BOUNDARY BETWEEN GAS 
EXTENDING TO INFINITY AT REST AND IN MOTION 


Fic, 1. 


The analytical expression for these r-characteristics is given by Courant and 
Friedrichs* : 


2 +I 
x= —  at— Y 
vs i 
where t = fy is the instant at which such a characteristic first meets the fan. Fora 
real gas, since y >1, this shows that the whole gas finally has the velocity given by 


lim (dx/dt),_,,=2a/(y—1). Thistakes the value 3a for y = §. 


t> @ 
4. Inhomogeneous initial conditions.—In order to remove the initial discontinuity 
of velocity, McVittief suggested the introduction of an initial inhomogeneous 


* R. Courant and K. O. Friedrichs, Supersonic Flow and Shock Waves, New York, 1948. 
+ Loc. cit. 


Aly (t/ty) 2+) 
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distribution of density in the gas cloud at rest at time t=0, but without any 
discontinuities. It is clear from a characteristics diagram that the inhomogeneous 
gas must instantly begin to move; for the intersection of r- and s-characteristics 
of differing constants immediately gives a non-zero value for the particle velocity w. 

Copson has shown, in his paper, that the variation of x and ¢ along characteristics 
in the region in which both r and s vary may be expressed in terms of an auxiliary 
function ¢ defined by 


&(r)=[ arla(r)le-odr (720), (=) =), (=) 


where [x(r)],_9 is the value of x at t=0 expressed as a function of the Riemann 
invariant r. Then 

_ 9 fbr) , H(s)—H(7)| _ 9 

= 5 * (ese fas (say) 


o eu oot __o% 


me et = — 5 (say), 


r+s (r+s)* Or 
and x =£+(4r —2s)t/3 =n +(2r—4s)t/3. 
On the face of the cloud r+s=o, since the velocity of sound vanishes. 
A characteristic r =7, meets the face at a time t =¢, given by 


$'" (1) = — 4h. (8) 
This formula was given by Copson and follows, by taking a simple limit*, from 
the expression for ¢. 
Now, if (8) is to give ¢, >0, clearly 4’ (r,) <0. At t=0, x(r)=¢'(r)/2r and 
r=c=1/w(say). After a little calculation the condition is found to be 


0?x 
Lae ..<° - 


If it is supposed that 7 is a monotonic decreasing function of x then a further 
necessary condition for preservation of the continuity of the motion is that the 
time t=f, increases with 7,; physically this requires that wavelets of the same 
family do not cross before reaching the surface of the cloud, but arrive at the 
surface later the farther they were from it att=o. Thus dt,/dr,>0, implying 
¢"")(r,) <0, or the equivalent inequality 

0*x 08x 
358 +58 >0. (10) 


The conditions (9), (10) are both necessary for continuity of the flow; they are 
not likely to be sufficient, except perhaps for simple distributions of [x(7)],_». 

5. Copson’s solution.—Copson discusses the use of Riemann invariants in the 
solution of the problem for y=§%, and ends with an example of an initial finite 
inhomogeneity between the vacuum and the semi-infinite extent of homogeneous 
gas. Hesupposes that [7(x)],_9 increases parabolically from zero at the front (x =h) 
to the value b=3a/2 at x=0, and remains constant for x<o. In 0<x<h, 
r=b,/(1—x/h). This distribution satisfies the requirements of the previous 
section and the corresponding solution is represented graphically by a diagram 
(Fig. 2). The initial inhomogeneity is confined to the strip OA. 

* The appropriate ¢ for negative values of s is ¢(—s). This satisfies the requirements of the 


continuation of the solution across the characteristic s=o and ensures a finite value of t for the 
arrival of the characteristic r=r, at the surface. 
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On the face of the cloud at any time t, c=0. Anr-characteristic which meets 
the face (a possibility in this case) reflects at this free boundary into an s-character- 
istic; since c =0, it follows that s= —r and dx/dt =u for both r- and s-characteristics. 
Hence, (i) the r-characteristics envelop the vacuum face in the (x, t) plane, and (ii) 
the velocity of advance of the cloud is equal to the local particle velocity (27). The 
maximum velocity of advance of the front corresponds to the maximum value of r 
(i.e. 3a/2); the maximum velocity is therefore 3a as before. In the characteristics 
diagram a section by t =0 yields a uniform region, a “‘ mixed’ region and a vacuum. 
In the uniform region both r and s are constant; hence, the boundary OD is 
straight and has dx/dt= —a. Let B be the point at which the r-characteristic from 
O touches the wave-front. From the discussion of the previous section, a simple 





VACUUM 





Fic. 2. 


wave (but no longer centred on a point) must fill the whole plane to the right of 
OBC, first of all (in time and space) because of the intersections of r-characteristics 
from a uniform region with the s-characteristics which come from between O and 
A; and then because of their intersections with the s-characteristics which are the 
continuations of r-characteristics after touching the front of the cloud (e.g. KLM). 
‘To prove that the characteristic OBC does, in fact, constitute the boundary from 
the point B onwards, consider an r-characteristic QM arbitrarily chosen from those 
to the right of OB. OB and such r-characteristics are all equally inclined to a 
given s-characteristic; hence QM meets the reflection BC at infinity. QM was, 
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however, arbitrarily chosen. Hence, OBC continues, after touching the front, 
as the front itself and BC is straight. From the time at B, t,, onwards, the mixed 
region is no longer present; a section at time ¢ >t, consists only of a simple wave 
lying between a uniform region and a vacuum. 

An example of an inhomogeneity extending to x = — 00 was given in a paper by 
Martin* which appears to have been overlooked by earlier writers. He gives the 
solution for y = $ and r= \/(—.x/2) at t=0 (x <0) in an ideal gas bordering upon 
avacuum. ‘This initial condition satisfies the necessary conditions for continuity ; 
the flow is in fact continuous and the lines of constant density (which, of course, 
include the surface of the vacuum) are parabolas. 

6. Alternative initial conditions—Two further examples of flow with y=3 
are illustrated in Figs. 3 and 4. Fig. 3 corresponds to the case where the speed 
of sound increases linearly between x =A andx=0. Equation (8) indicates that ¢, 
is constant for this case, and in fact the r-characteristics of the mixed region now 
converge toa point B. A point-centred simple wave at B again reduces the pressure 
at the wave-front to zero. ‘This is the limiting case of Copson’s solution when the 
envelope shrinks to a point. (Burgers’ solution is the limiting case when the 
initial spread of c is confined to a point.) 
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Fig. 4 illustrates a further possibility. Here [x(r)],.=(I1-7)?, I>x2A 
{h>0).t This satisfies the condition (9) only for certain values of r and never 
satisfies (10); in this case, therefore, the flow cannot remain continuous. The 
r-characteristics now possess an envelope inside the gas cloud. The inter- 
section of characteristics of the same family gives multiple values of u and c at the 
points of intersection; this indicates the occurrence of a shock-wave within the 
flow, beginning at the point and instant when the one-valued character of the flow 

* M. H. Martin, Am. }. Maths., 65, 391-407, 1943. 


+ The initial point of the shock-wave in Fig. 4 occurs at time t() where t(h) moves continuously 
towards zero as h tends to zero. When h=o the shock-wave begins instantaneously at t=o, 


4 





50 D. C. Pack Vol. 113 


breaks down. ‘To a first approximation the development of the shock-wave is as 
shown by the envelope of the characteristics in the diagram.* In an exact 
treatment the variation of the strength of the shock-wave with time has to be taken 
into account, the r-characteristics no longer being tangential to the shock-wave 
after the initial point. A particular case has not been calculated in detail on 
account of the labour of computation. Certain deductions may, however, be 
made by a brief examination. 
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As gas crosses the shock-wave it gains in entropy. Ajjdotted line in Fig. 4 
represents (diagrammatically) the path of the particle which was at the point of 
origin of the shock-wave at the critical instant when it began; this line separates 
those particles which have passed through the shock-wave and experienced an 
increase of entropy from those which have not. When the shock-wave reaches the 
front of the cloud a pressure discontinuity exists there; this must, as in Burgers’ 
problem, resolve itself instantly into a wave of expansion. The effect on the 
velocity of advance of the wave-front may be seen from the condition of compat- 
ibility for a characteristic which passes to the right of the shock-wave. Along 
such a characteristic the relations of Section 2 require modification, and the 
condition of compatibility is in factft 


s 
3e+u— | R?cdS =3a, 
So 


where S denotes entropy, Sp its value in the gas in its original state, and ? is a 
positive constant. Since, by the second law of thermodynamics, S can only 
increase, and since ¢ is an essentially positive quantity, the value of u when c 
vanishes is greater than 3a. 


*In the complete mathematical solution the r-characteristics are tangential to the envelope 
and the s-characteristics are cusped where they meet the envelope. In the diagram the 
characteristics are only taken up to the envelope; their continuations beyond these points are of no 
physical significance. 

+ See, for example, Pillow, loc. cit. 
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7. Conclusion.—The motion of a semi-infinite gas cloud bordering upon a 
vacuum has been examined by the method of characteristics. Characteristics 
diagrams make physical details of the flow easy to see, and the problems of Burgers 
and Copson have been examined by means of them. ‘These problems concern a 
gas with ratio of specific heats y = , corresponding to gas clouds of monatomic 
hydrogen. Burgers’ problem involved an initial discontinuity of density at the 
front of a homogeneous cloud; Copson’s stipulated an initial region of inhomo- 
geneity between the vacuum surface and the main homogeneous cloud, without 
any actual discontinuities in density. Burgers’ solution is continuous for t>0; 
the particular example given by Copson is also continuous. It has been shown 
above that unless certain conditions are satisfied by the initial state in the cloud the 
subsequent flow cannot remain indefinitely continuous. An example of the break- 
down of continuity is given. 

An initially linear distribution of values of the speed of sound in the gas at rest 
leads to a singular solution in which all the characteristics of one family intersect at 
a point of the wave-front at the same instant of time t= 7 (say). The solution is 
subsequently continuous. 

It is shown that, when the flow remains continuous, the wave-front ultimately 
moves with velocity 3a. When a shock-wave appears in the flow the wave-front 
attains an even higher velocity. Thus, the high eventual velocity of advance of 
these gas clouds, which was considered a physically unsatisfactory feature of the 
previous solutions, is not to be avoided by variation of the initial inhomogeneities. 


Institute for Fluid Dynamics and Applied Mathematics*, 
University of Maryland, 
U.S.A. : 
1952 October 1. 


* On leave of absence from University College, Dundee, University of St. Andrews. Now at 
the University of Manchester. 
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Summary 


When the scattering of light is coherent, the exact solution of the problem 
of line formation in the Milne—Eddington model, with a linear Planck function 
B,=b, +6, 7,, has been given by Chandrasekhar. In Section 2 this 
solution is found by a new method, and in Sections 3-5 the method is used 
to solve the parallel problem when the scattering of light is completely 
non-coherent. This solution is exact when 6," and «,6," are independent 
of the frequency v and the depth z, «, being the continuous absorption 
coefficient. The solution agrees with that found by Sobolev for the case 
Kk, =k, b,=B, b,™=o. In Sections 6-8 the magnitudes of the terms 
appearing in the solution are discussed. Sobolev’s assertion that the 
resulting central intensities of absorption lines are much lighter than those 
found on the hypothesis of coherent scattering is shown to be true for lines 
at the limb. No definite conclusions can be drawn for other lines until the 
H-functions involved have been computed, but the same result is likely to hold. 





1. Introduction.—In the theory of line formation, it has been assumed, until 
recently, that the scattering of light is coherent, i.e. that the frequency of the 
radiation re-emitted within the line depends only on the radiation absorbed at 
the same frequency. ‘There are, however, several causes leading to a redistribution 
in frequency, and an investigation of the effects of non-coherent scattering has 
become necessary. The reasons for this have been fully discussed by G. Miinch 
(5) and (more recently) by M. P. Savedoff (6). Comprehensive lists of references 
will be found in both of these papers and Savedoff reviews all previous attempts 
at solving the problem of line formation when the scattering is either partially or 
completely non-coherent. Unfortunately he dismisses the important paper by 
V. V. Sobolev (7) as ‘‘extremely complicated” and ‘“‘unclear’’. Sobolev’s 
paper is certainly ‘‘unclear’’, but it gives the first exact solution of the problem. 
This is for the case in which scattering is completely non-coherent and the Planck 
function B, and the continuous absorption coefficient x, are both independent 
of the frequency v and the depth z. 

Savedoff’s own method is complicated. He also deals with the case of 
completely non-coherent scattering and he begins by transforming the equation 
of transfer (see Section 3) into a non-homogeneous integral equation. He then 
uses iteration to obtain an approximate solution when the Planck function is 
linear. ‘The resulting residual intensity in an absorption line gives interesting 
results for the wings of the line, but appears to become indeterminate in the 
centre, the place where the discrepancy between theory and observation has 
always been greatest. This is a weakness in Savedofft’s method, though I do not 
think it should be difficult to overcome the indeterminacy. 

In (7), Sobolev asserts that the central intensities obtained from his exact 
solution are a great improvement on those obtained with coherent scattering. 
For this reason his paper deserves serious consideration. 
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Sobolev’s paper is written in Russian; his assumptions are not clearly 
stated; the paper contains misprints, and his deductions are questionable, because 
he neglects terms which are, as a matter of fact (see Section 8), of the same size 
as those which he retains. 

It was whilst trying to elucidate this paper that I began to investigate 
extensions of Sobolev’s results. The case of partially non-coherent scattering 
can be approximately represented as the sum of coherent and completely 
non-coherent scatterings in a suitable ratio, and I therefore restricted myself 
to the case of completely non-coherent scattering, but I replaced Sobolev’s 
assumption that B, is a constant by the assumption that B, is linear. 

The method used by Sobolev breaks down when B, is not constant. In order 
to solve the problem for a linear B,, I had first to evolve a new method of solution 
for the corresponding case of coherent scattering. This is given in Section 2. 
It uses Chandrasekhar’s ideas, but is much simpler than that given in (2), §84, 
and it can be applied when the scattering is non-coherent. This is done in 
Sections 3-5. 

The solution of the problem (equation (46)) is too complicated to be of much 
value as it stands. In Sections 6 and 7, I investigate the functions involved in 
considerable detail. Finally, in Section 8, I consider the sizes of the various 
terms appearing in (46) and I show that Sobolev’s assertion about the improvement 
in the central intensity is true for lines at the limb. No final conclusion can be 
drawn until the H-functions involved have been computed. 

In an appendix I give an alternative method for finding the scattering function 
of Section 4. The second method is longer but is easier to follow. 

2. Coherent scattering.—With the usual notation, the equation of transfer for 
the Milne-Eddington model can be written 

1 
p= (+ odor [Ide KB (1) 
pdz -1 
where z is the depth below the surface, «, is the continuous absorption coefficient 


and o, is the line-scattering coefficient. We assume that x, and o, are independent 
of depth and we write 


z z 
i= | p(k, +o,) dz, 1, = | pk, dz, (2) 
0 0 
yy Kyy A,= I (1 + Ny) = Ky) (x, +9,). (3) 


Then 7,=A,¢, and 
B, = bo af bz, _ b, + r,b,t,, (4) 


if we take the first two terms in the Taylor expansion of B,. Substituting into 
(1) we get * 


dI,(t,, 1 7m | ; 
poet —ltnp)— HAA) J Taltt! dy A (6604208). (3) 





* This is (2), § 84, equation (24), except that Chandrasekhar drops the subscript v and allows 
for thermal emission. We shall retain v because it becomes a variable in Sections 3-8. The result 
of allowing for thermal emission is to change the definition (3) of Ay into Ay=(1+€yy)/(1+W)- 
Then t= Apty (1 +€yy) and 4, has to be replaced by 4, /(1-+e 7»). On making these changes in 
the final solution (16), the required solution is obtained. 

The same changes apply to the equation of transfer (19) for non-coherent scattering, but 
because this still contains o,’, the changes cannot be immediately applied to the solution. It is not 
difficult, however, to introduce the necessary modifications into the subsequent analysis. We 
shall take €,=o0 throughout. 
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The equation (5) has to be solved subject to the boundary condition 
I(0,-p')=O = (0<p'<1) (6) 
together with the condition (implicit in Chandrasekhar’s solution) that /,(¢,, 1) 


is at most linear in ¢, for large t,. Now a particular solution of (5), which does 
not satisfy (6), is 


L{ty 1) =b,9 +A bt, +4); (7) 
as is readily verified by substitution. We therefore write 
L(t, 2) = 6,9 +A,b, Mt, +2) +1,%(t, 1). (8) 
Then /,*(t,, ) satisfies the integro-differential equation 
dl,*(t,, ; a, 
Sete) 1,644) Mra) {Tetltne!) de (9) 
together with the boundary condition 
1,90, —w')= Abu’ —6, — (0<p' <1). (10) 


Moreover /,*(t,, 4) must be at most linear in ¢, for large t,. 

We require the emergent intensity /,*(0,u). This is the sum of /,*(0, 1), 
where /,*(t,, 4) is the solution of (g) subject to the boundary condition 

I,*(0, —p’)=0 (o<p' <1), (11) 

and J,*(0,), where this is the diffusely reflected intensity corresponding to the 
incident intensity (10). 

It can be shown that f, unless A,=0 (which is not so), 

I, he H) 0. 


Hence I,*(0, w) = 1, *(0, 2) 
5 1 
cenget it S ’ , A,b, ’- 5, d € I2 
=| (1, o')( te )dp (12) 
where by (2), §33, equation (41), 
S(u, pn’) =(1-A,)  A(w) Av’ 13 
(15 e') =( ea (u)H(u’) (13) 
and H(,) is the solution of 
* H(p’) 
H(z) =1+ (1 —A, eH | Ala. (1 
(1) ( )pH() we 4) 


From (12) and (13) we have 


‘H(p') 
i, ’ _ 4 —A, H | Ab, Poe b,© e ; d , 
(0,4) =3(1 )H(u) a! be ) yn (n 


1 2 
= (rAd) f (wn 2) He) dy 
1 
= x=A,)6,H(u) f(x Eo) Au) de 





+ This can be seen as follows : with the notation of (2), § 26, 


n ~ ligt 
1,*(ty, w= H(1— Ay) % lat */(1 +k), 
a= 
where the L, are determined by 


n 
xX Le/(1—pika)=0 G=, «005% 
1 


a= 


n n 
Since II (1—pkg) x EZ Le/(1—pke) is a polynomial in uw of degree n—1 with n distinct zeros, it is 
i= = 


a 
identically zero. Hence every Lg=o, and in the limit, as n>, 1,*(t», u)=0. 
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and by (14) this is 
T,*(0, 2) =(b, +A,b,%u){H(n) — 1 — 3(1-A,)tgH(u)} + 3(1—A,)A,F,2,A(H), (15) 


where a, = [ wrt) du. 
By (2), §38, Theorem 1, . 
es (1 wis A,)% =A,?, 
and so (8) and (15) give 
T,(0, 4) = (6,0 +A,b,%u)A, 7 H(1) + (1 —A,)A,b, M2; H(z). (16) 
This is Chandrasekhar’s solution ((2), §84.4, equation (66)). 
3. The equation of transfer in the case of completely non-coherent scattering.— 


When the scattering is completely non-coherent, the scattering term on the 
right-hand side of (1) is replaced by 


foo) 1 ioe) 
ho, | o, dv’ | ly a'|| o, dv’, 
0 -1 0 


and the equation of transfer is 
dl wo 1 ~ 
m ~ =(x,+0,)1,— 0, a oy dv’ bt ly ay | oy dv’ —«,B,. (17) 


In this we can, without loss of generality, take 
«© 
J gd «%. (18) 
0 
In doing this we are really replacing o,, x, and p by 


@ @ @o 
o,f o, dv’, «| o,dv’, p| oydv’, 
0 0 0 


but /,, 7,, 7, and A, are unchanged. 
On using (2), (3), (4) and (18), (17) becomes 
dI,{t,, #) : ts se @ w 
K ~ dt, = I ft,, H)— 3(1 —A,) I, oy dv [tlt )du —A,(b, +A,p, t,). (19) 
This has to be solved subject to the boundary condition (6) and the condition 
that /,(t,, 2) is at most linear in t, for large t,. 

In order to find J,(0,) by the method of Section 2, we have to find a 
particular solution corresponding to (7) and a scattering function corresponding 
to (13). We shall deal first with the latter. 

4. The scattering function for completely non-coherent scattering.—In this 
section | shall follow the method due to Ambartsumian (1) which was used by 
Sobolev. This method is explained in great detail in (3), §31, and I have 
introduced modifications in Sobolev’s analysis in order to bring the notation 
into line with that of (2) and (3). An alternative method following that used by 
Chandrasekhar in (2) is given in the Appendix. 

In Sobolev’s paper (7), «, is assumed to be independent of v but this is an 
unnecessary assumption. 

Consider a semi-infinite, plane-parallel, scattering atmosphere subject to 
incident radiation of frequency v, all of whose rays make a fixed angle cos pg 
(O0<pyp <1) with the inwards-drawn normal. Let the monochromatic flux per 
unit area perpendicular to the rays be F,. Then the incident intensity for any 
direction —y’ and frequency v’ is represented exactly by 


1,0, —p') = $F, 8(% — v’)8(u9 — #’), (20) 
where 6 is Dirac’s 5-function. 
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Let the scattering function for the atmosphere be 


S(v, v', oye’), 
so that the diffusely reflected intensity J7*(0,) of frequency v and direction p, 
due to an incident intensity /°(0, —,’) of frequency v’ and direction —,’, is 
given by 
(0,0) = dv’ f 19°(0, —p') SD ay (21) 
0 0 2p. 
In particular the diffusely reflected intensity due to the incident intensity (20) is 


F,,,S(», Vo by Ho)/4t 
Consider the passage of the radiation through the layer of thickness dz below 
the surface of the atmosphere. The equation of transfer is (17) without the term 
x,B,. This can be written 


> oo 1 
dl,= Oe Ceo odo’ { Ty dp’, 
Be 2p 0 -1 


from which the absorption and scattering processes can be seen. By the 
principle of invariance, radiation travelling downwards through the layer is 
diffusely reflected in the lower surface of the layer according to equation (21). 
After tracing all processes in the passage of the radiation, we get the following 
equation (to the first order in dz) for S(v, vo, 4, Uo): 


S(v, Vor Hy Ho) _ S(v, Posteo) fy _( er 2] [2-¢. pom we 4 210 n,p dz 
4e 4u! 





: d, d. 
ened (* o, dv’ | Sv, v's HH) v> HH) = > ba hn a o, dv’ 
4 Jo Jo 2p B 44 Jo 
rl , , dz oe) . 1S ’ Z . ” d ” 
. | Ss He Pio a) ays 5 OS sit { (% ¥"s yw”) du” 
0 2p 4 Jo 0 2p p 
cs) 1s Pe , 
x | ovdv’ | (v', ey ; Bo) 

0 0 ou 

In the limit, as dz 0, this gives 


S(v, ae ered el Ce dv’ eee ae 














Oy Oy, B Ho Gey - 
If? S(v’, vo, d 
x (: - | oy” aw’ f (v's Yo» Hs Ho) “) : (22) 
- 0 0 Oy Oy, B 


By the principle of reciprocity, 
S(v, VY; Bo» b) = S(v, V, My fg). 


From this and (22), by a rather lengthy argument which I omit, it can be 
shown that 


S( Vg Vy Mg) = S(v, Vo» Bs Ho), 


v ") du’ 
bormy=r4 3[ oy dv’ [eee (23) 


o,0y a 


and so, if 


then (22) gives 
S(v, Vos ity Mo) m Lio h(v, L)P( Vo, Ho) 
O,Oy, (a, +K,)ptg + (oy, +K,y,) 
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On substituting from (24) into (23), we have 


ov’, w’) ; ‘ 
$v, u) =1 + dudlv, 4) [0 2 dy’ leone rt (2) 


From this it follows that 4(v,) is a function of p/(o,+x«,). 








Let 
k=min(«,+<¢,) (26) 
¥ 
and write 

Kp i 
(= > x= , 27 
K, +o, Ky +oOy ( 7) 
H(x) =4(v, 1). (28) 


Then (25) becomes 
_ & 
H(x) =I+ Pa H(x) f°. 


When the order of integration is changed in (29), the limits for v’ have to be 
changed to ones depending on the form of the curve 


K/(oy+Ky’) 
2 dv! } H(x') 3. (29) 


x+x' 


x =k/(ky+o,) («= min (x, +0,)). 


The line scattering coefficient o,, is symmetrical about the centre v’ =v, of the 
line considered and it decreases rapidly to a very small value on either side. No 
appreciable error is therefore introduced if (29) is inverted into 


1 
H(x )=14+)xH(x yf A(x) Ke) ay (39) 
x+X 
where 
2 ao 
K(x’ ==| oy di (31) 
KJ x2") 
and 
v(x’) =Y, if x’<kK (x, +o,,), 
22) 
N=K (Kye) t+ Onyx) if x’ >K (x, +9,,). (32) 


In terms of the new variables, the scattering function (24) is given by 
o,oy xx’ H(x)H(x') 
ae ae 

where x and x’ are given by (27). 

If x, is a constant for all v, then x,=« since o,—-0 as |v—v,|> 0. The 
formulae (30)-(32) are then identical with those of Sobolev. 

When the function K(x’) has been determined, the solution H(x) of (30) can 
be found by the formulae given in (2), Chap. V. Note that the H(x) of (30) is 
not the same function as the H() of (14). 

5. The solution of equation (19).—If, following the method of Section 2, we 
look for a particular solution of (19) similar to (7) by substituting 


I(t, 4)=A,+ B,t,+C,u, 


where A,, B,, C, are independent of t, and y, we find that there is little hope of 
succeeding, even with more elaborate functions, except in special cases, and 
then (7) is itself a solution. Let us therefore write, as in Section 2, 


T(t, 2) =5,% +A, b(t, +p) + 1,%(to H)- (34) 


: (33) 








58 Ida W. Busbridge Vol. 113 


‘Then on using (18) it is found that 7,*(t,, 4) must satisfy the equation 


* - . 
ad cat Pa ») =1,°(t,, fh) — 3(1 -2) ode’ | Tete) ae +8(ty 1), (35) 


where 


g(t,,u)=(1—-A, oy{(b,% — b) + (Abt, —A ob t,,)} dv’ 
wa 


@ Zz 
=(1=A) {ap{(b,0— 6) +(b,,— by) | pds} de’. (30) 


Now ay is greatest at v’ =v, and decreases rapidly to a very small value on either 
side. If v is near to »,, 6, — 6, and b,%, —b,«, are both small when a, is not 
small, and the integrand in (36) is small over the entire range of integration. 
This suggests that g(t,, ) is negligible for values of v near to v, and those are the 
important values in the theory of spectral lines. 

The function g(t,,u) is identically zero if 6, and 6,'°«, are independent of v, 
and then the following equations are exact. Sobolev’s results (but not his 
method) are obtained by taking 5, = B (a constant), 6,” =o. 

If we neglect g(t,,u), we have 


dl,* t,, fe) : cl ’ , 
# » =1,%(t,,n)—3(1 -A,){ teal ree an 
v 0 . 


together with the boundary condition 
T,*(0, —p’)= Ab Ow’ -—b,O (O<p'’ <I, O<v'< 20) (38) 


and the condition that /,*(t,,.) must be at most linear for large t,. Now, as in 
Section 2, the only solution J, *(t,,) of (37) subject to the boundary condition 


I,*(0, —p’)=0 (o<p'’ <I, O<v'<@) 


is J,*(t,,4)=0. Hence the required solution of (37) is the diffusely reflected 
intensity due to the incident intensity (38), i.e. by (21) 


fe) 1 ‘ , 
T,*(0, Ht) = j, dv’ [. Ayb Pp’ sa b) set du’. 


By (33) and (27), this is 


re cK Ky +Gy’) a! f 
1,*(0,n) =3H(x)——— | oy dv’ | ae 


Oytkr lg 0 x+% 


x [aedeMer toy) = — 4 | egal 


K 


Sei A £2 K/|(K,¢+ Oy") x’ H(x’) x’ 
= }H(x) —— , Jo, dv’ cles all | FT Oa 
2 (x) pe \, (x, +o, Jo, dv l x+x ( i Ky we v ) x 
by (3). 


Since o, is large only near v’=v,, we can again replace b,«, and 6, by 
b,%, and 6, when v is near to v,, and we have approximately 
b,%, ae , wy +0y) x”? H(x’) 

= (xy +oy)o,dv 7 
is 0 0 x+X 
why+0y) x Hx) 
0 +x’ 


1,9(0, 4) = M(1—A,)H(x) dx! 


5, 


— 41-2) H(x)-— | : (cy +oy)orde" | dx’. 
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Let * 
K(x')=2 | kyoy de’, (40) 
K J oz’) 


where v(x’) is given by (32). Then inverting the order of integration in (39), we 
have, by (31) and (40), 


1,4(0,1) =x) =" (a) [ACD cK) + Kye} 


-4(1 anon [ eat {K(x')+Ky(x')} de’ 
o *+x 
This equation is exact when b, and 6,°%«, are independent of v. 


The function K,(x’) is related to K(x’) by the equation 
1 1 
{ K,(x’) dx’ + J K(x’) dx’ =1; 
0 0 


1 If? K/(Ky + oy") rs) 
J Kewyae'=2[o,2ar’ | dv’ = [° ay, 
0 KJo 


0 0 Ky a 
1 «© (Hy? + Oy’) «© 
Py I , id OyKk, ‘ 
J K,(x’)dv' == { —_ f dx’ = | ei, 
0 KJo 0 0 Ky + Oy 


and adding we have (42) by (18). By (3), (43) can be written 
1 oo 
{ K,(x') dx’ = | eee eA 
0 0 
Since, by (2), Chap. V, Theorem 1, 


1 1 1/2 
I. H(x') K(x’) dx’ =1— [: - j, K(x’) de’ | A 


; j. H(x’) K(x’) dx’ =1—- (J K(x’) a’) : : (45) 


the terms in (41) containing K(x’) can be simplified by the process used in 
obtaining (15). — ~~ the resulting expression into (34) we get 


( 
Io, u)= (6. oak x) {(-A)Ha (f K\(«’)dv’) ‘ +} 
+ 


(1 —-A,)H(x {es an Hx’) K(x’) dx Eel ee) ae 


x’ H(x')K,(x’) 
mt et fF - 8 
y ik x+x’ sal (49) 
This is the exact solution when 4, and «,b,” are independent of v. 
Since x =xp/(x,+0,) =Kpd,/x,, therefore 


OQ, 
5, 4+ 2 x =F, 4 ADO, 


or, by (42), 





This shows that there is a sieileciey of form between (46) and (16). 

We must now investigate the behaviour of the tunctions K(x) and K,(x) and 
the sizes of the various terms appearing in (46). Jn all that follows we shall assume 
that x, =k. 

* A very small error is introduced when K;,(x’) is written in this form, just as a small error is 


introduced when K(x’) is written in the form (31). These errors are removed whenever the 
integrals are re-inverted (as in the proof of (42)). 
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6. The functions K(x) and K,(x).—When x,=« we have, from (31), (32) and 
(40), , 


«© 


K(x) == C 02 dv, (47) 


i 3) 
z 


K,(x)=2 J. 2,4 (48) 


where 
v(x) =v, if x<1/(I+%,), (49) 
x=(I+ ny)! if x>1/(I+ Ty,)s nia 
and y, =o, /«. This will be large; Sobolev, in an example, takes , = 10°. 
In the notation of (4), 
a” "in d 
=O, -— ry 50 
_— | > ) (5°) 
where 
w = 2(v—v,)(In2)!?/Avp, a= Avy(In 2)!?/Avp. 
Since (18) must be true, 
@ @ A oo ‘ 2) dw 
Sf re ee 
and this at once gives 
2(In 2)!?/Avp = we... 
From (51) we now have 
= mg, (v— ve)» (52) 
The classical values for a in the visible part of the spectrum are small; values 
given in (4) range from 0-001 2 to 0-035. However, P. ten Bruggencate and 
J. Houtgast have shown (8) that the observed value of the natural damping 
ratio may be as much as thirty times the classical value and so a may be as large 
as 1°05. The following discussion only applies when a is not greater than 0-01, 
since then (see the appendix of (4)) o, is given approximately by 
o,=0,¢ all for |w| <a, 
(53) 


=o,ar'*a* for |w| >wo, 


Co 
¥ % 


where, for continuity, 
e @,* ar ew, 2. (54) 
This equation gives two values for w,* and we must choose the larger. It is 
impossible to find an approximate formula for this root and it has therefore to 
be evaluated for each value of a. 
In the following discussion I shall estimate the magnitudes of terms when 
a=o-0I, Ny, = 108. 
This is the example chosen by Sobolev. In this case 
W,) <= 2°67. 
First consider K(x). Let 
w(x) =o, | v(x) — v,]. 
On changing the variable in (47) from v to w by means of (52), we have 
, 2 St ei 
K(x)=—= - — o,” dw. 


1/2 
7 KO», oz) 





No. I, 1953 Coherent and non-coherent scattering 


Ifo<x<(1+7,) 1, then o(x)=y, and w(x) =0. _ Hence by (53) 
° dw 
K(x)= am ed + eh -{°S} 
= Ny, {27" "PL 2w9| + 57 8° a? we *}, (59) 
where P(x) is the probability integral defined by 


af «+ 
P(x) = ze - dt. (60) 
For the value of a given by (5 5), Pl 2w,] differs from unity by less than 10~* 
and §7~*?@’w,. *=0-63 x 10°. Thus K(x) is approximately 2 120," 
If (I+, ) '<x<Xq, W here %q is défined by 
wy = 71a, [v(xX9) —»%], (61) 
then 0<w(x) <wy and, by (58), (59) and (53), 


& 2 / 9 9 ¢ 2 Az) a 
K(x) = N»(2 1 as V 2wW| + gn-8 a" wo *) _ zie 2», j ee?” dw 


= ny, {2 ¥? Pl V/2w9] — 27-1 * Pl /2u(x)] + 37-97 a?w9 3}. (62) 
By (49), (61) and (53), ' . 
Xo =(I+ M42) '=(1+7,,€ % 7. (63) 
For the values given by (55), 
Xp = 0-001 3. 
Thus K(x) decreases rapidly in the very small interval [(1+7,)~1, x9] from 
about 0-7 x 108 to about 0-63. 
If xy<x <1, then w( - >w, and 
a dw , - 
Kam S| Se den, atlo(s)) (64) 
Now from (53), when x>Xp, 
[w(x)] y= 7V2(an,_) Nu2) 
and hence by (49), since x)» >(1+7,,)"), 
[w(x)}* =2¥#(an, tx (1 — x). (65) 
On substituting for w(x) in (64) we have 
K(x)==§a-*4(a/n, )Y?x-¥2(1— x)? (xg<x <1). (66) 
For the values (55), K(x) decreases from about 0-63 to 0 in the interval (xp, 1). 
In exactly the same way we get the following approximate values for K,(x): 
K,(x) = P(wo) + 27 *aw_* (0<x<(I+7,,)"), (67) 
K,(x) = P(wo) — Plo(x)] +27 awe * = ((I+7,,)*<* <%), (68) 
Kj (x)= 22 *4(a/n, )¥?x-¥?(1 — x)? (%»<x <1). (69) 
For the values (55), 27 1awy }=0-002 4, P(wp)=0-999 8. Thus K,(x) decreases 
from about I to 0-002 4 in the interval ((I+7, )-?, x9) and from 0-002 4 to o in 
the interval (%p, I). 
1 


We next need J K,(x)dx and for this it is easiest to use the formula (43). 
0 


When x, =x, this gives 


1 r© og 
| Kils)dx= J. it 
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By (52) and (53) 


1 2 % dw . dw 
K,(x) dx —; J —=— + | SoS 
J, if ) qrll2 { o er + Ny, ei mil2gq-lyy? + Ny, 


2 bie I I a\le 
= ow we aoe dw +74 { — : 
wiJo (EP +H, wea w'+Ny, Ny, 


The integral is of the order of »,~* and is therefore negligible compared with 
the second term. Hence 





1 
J K,(*) dx =z! ‘(a/n, )r”. 
0 
Following Sobolev, we shall write 
. 1/2 
A= ( J Kus) dx) = m%(a/n,)". . 
0 


For the values (55), 
Axo0-0I. 
7. The function H(x).—This is the solution of (30), or, by (2), p. 107, 
Corollary 2, of 
B5 > I ('x’H(x’)K(x’) _, 
H(x) -A+3). x+x’ o. 
As for all the H-functions, H(o)=1 and H(x) is continuous and steadily 
increasing in the interval (0, 1). From (71) it follows that 
1/H(x)>A 


(71) 


and therefore that 

I<H(x)—<A for O<x<t. (72) 

Sobolev gives the following approximate formula for H(x): 
2x/A* 


H(x)=1+ 1+2x/A° 


(73) 
I have failed to find adequate justification for this approximation, which is 
stated without proof. It makes H(o)=1 and H(x) a continuous and steadily 
increasing function in (0,1), but it makes H(1)>A™~!, which contradicts (72). 
Moreover the derivative of H(x) is infinite at x=0, whereas that given by (73) 
is 2/A*. If (73) does give a sufficiently good approximation for H(x), then 
H(A) >§A>1. (74) 
In discussing the magnitudes of the terms in (46), I shall only use (74) to 
show that, if Sobolev’s approximation (73) is satisfactory, then the term which 
he neglects is not negligible. If we work solely with the wider inequality (72), 
we can only show that the term may not be negligible. 
8. A discussion of the solution (46).—Sobolev’s solution is (46) with x,=«, 
6, = B, 6, =o, and in his deductions he neglects the term 
f x H(x')K,(*') dx’. 
0 | 6&+% 


(75) 


I shall show that if (73) is a satisfactory approximation, then (75) is of the same 
rl ¢ 


> 


1/2 
order of magnitude as A= (| K(x) dx) , the term retained by Sobolev. 
0 


If we can only use the inequality (72), the proof that (75) is not larger than A 
still holds, but the proof that it is not much smaller breaks down. 
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Since all the functions appearing in (75) are positive, the integral is greatest 
when x=o and we therefore want the size of 


J * H(x!)K,(x') dx’. 
0 


7 1 he 
By (72) and (70), J Hw K(x!) dx! <A J Kula’) a =A. (76) 


Hence (75) is not greater than A. 
When a and », have the values given by (55), then x) =0-001 3 and A=o-01, 


and 
oyeen A>Xq. (77) 


From (63) and (70) it is easily seen that (77) will be true for any fixed a (and 
therefore fixed w,) provided y, is large enough. We shall assume that this is 
the case. 

We have 


J * H(x’)K, (x!) de! > J * H(x’)Ky(x’) dx’ > H(A) J " K(x!) de’, 
since H(x’) is steadily resintiiaile (A, 1). By (77) and (60) 
J H(x’)K,(x’) dx’ > H(A) . 20-84(a/n, )} J é x-U2(1 — x)? de 
= 4?H(A) 245 =sin-at?— 4%(1— ay} 


>4A?H(A), 
since A <o-I in all the cases we are considering. If (74) is true, then it follows 
that 


J " H(x')K,y(x’) dx! > 4A. 
.., 


Thus the term is not negligible if Sobolev’s approximation holds. Since H(A) >1 
(by (72)), the term lies in any case between A and }$A*. 

We shall now investigate the other integrals appearing in (46) in a similar 
way. We shall, however, only obtain an upper bound for each integral. 


The integral f x'2H(x') K(x!) = 
0 X+X 
is less than (75) and is therefore less than A. 
The integral : 
J x! H(x!)K(x") dx’ 
0 


requires more detailed consideration. By (72), 
“1 1 

| x’ H(x’) K(x’) dx’ < 4{ x’ K(x’) dx’. 
0 0 


On substituting from (31) and inverting the order of integration, we have 
ri. io ied 1l+n,)* 

| x’ K(x’) dx’ = aa a,” dv’ f x’ dx’ 

0 0 


0 


re 2 
oy" 
= a? | ——— v 
o (I+n) 
a) o,” ‘ 


—x (I+7,)* ”. 


be | 
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By (52) and (53) 





1 ; We dw * dw 
, , P st o~ii2 iene 
J, el ) ” sei "We { J, (e + Ny) ’ We (7?a7tw* + at 





12y {"{ : = : \ de 
ve 0 (~ vs ny) (x? 2a- yy “- 1») 
+42*(a/n,,)!*. 


The first term is of the order of n,~? and is negligible compared with the second 
term, which, by (70), is }4*. Hence 


“1 
| x’ K(x’) dx’ = 4 A* 
0 
and (79) gives 

“1 


| x’ H(x") K(x’) dx’ <4 A. (So) 
0 


From these results we see that no integral in (46) is greater than A, but, on 
the other hand, there is no reason to suppose that any one is negligible compared 
with 4. In view of this we shall now consider whether Sobolev’s assertion about 
the central intensity of an absorption line is justified. 

When v=», then x=p(I+7,)* and this is very small, but it does not 
necessarily follow that H(x)==H(o)=1, because H(x) increases very rapidly 
near x=0. For this reason we shall only consider the case in which x =o. 

Since A, =(1+7,,)", this is negligible compared with 4, and (46) gives 


rl -1 
I, (0, 0) +b, © {4 -—3 | H(x’)K,(x’) ax’| + 1,24 | x’ H(x’) K(x’) dx’ 
0 


0 
-1 
+ | x’ H(x')K,(x’) as'| ‘ (81) 
0 
To the same order of approximation (16) gives (with an obvious notation) 
(coh), ~ _— —1/2 
_ (0, 0) ~~ b, OA, } 2==b, On, , =, 
Now by (76) and (70) 
“1 
A-}| H(x')K,(x') dx’>$A=}o"(a/n,)!*. 
~0 
Hence 
I, (0, 0), Tieomo, 0)>42"8(an, M5 


for the values given by (55). Thus the central intensity of a line at the limb 
is much lighter on the hypothesis of non-coherent scattering than on that of 
coherent scattering. A similar result can be expected to hold when p>o, but 
it is difficult to establish this without computing H(z). 


In conclusion I should like to thank Professor H. H. Plaskett, Dr M. G. Adam 
and Mr D. W. N. Stibbs for their interest and for helpful criticisms. 
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APPENDIX 
The scattering function found by Chandrasekhar’s method 


We consider the same scattering atmosphere as in Section 4 subject to the 
same incident radiation. In this section it is convenient to denote inwardly 
directed intensities by the superscript — and outwardly directed intensities by 
the superscript +. Then the intensity of the incident radiation is, by (20), 


1y-(0, — 1") =4F,,5(%—¥")8(uo— K’) (82) 


the law of diffuse reflection is given by (21), and the intensity of the diffusely 
reflected radiation is 


I,*(0, w) = pF, S(v, Vor Py Mg) / [Me (83) 

Because t, varies with v, it is necessary to work with the depth z and we shall 

therefore write /,(z, ~) where we have had /,(t,, ~) hitherto. The equation of 
transfer for J,(z, «) is (17) without the term «,B,, viz. 


dI1,(z, ‘og Ts be ais 
pH) (6, +0), n)—3,{ ord’ [TAs n')dw’. (84) 
pas 0 -1 
In this we have used (18). 
We now define the intensity of the reduced incident radiation by 
1,¥(3, 1) = 3F,5(v9 — v)8(Hg + we Cre ondiee, (85) 


Then /,%(z,u)=0 except when v=vg, “= —py, and it satisfies the transfer 


equation 
dI,*(z, HK) R 2 
a = (0 + z, pu). 8 
po = («,+0)Me, H) (86) 
Let 
5,(z, uw) =1,(z, w) —1,%(z, p). (87) 
Then §,(, ) is the intensity in the reduced radiation field, and by (84) and (86) 
it satisfies the equation 
d3,(z, ) _ ” " 
Mods =(k, ‘ o,)3,(z, H)—o, Xz), 


where 


Ire F 1 ; 3 
B)=5| ord | Ten!) dw’ 
0 J-1 
Also, by (82), (85) and (87), the boundary condition for 3,(z, ;) is 
3, (0, —p’)=0, 


and by (83), (85) and (87) 
3, *(0, Bh) =1,7(0, p) = }F,,S(v, Voy Hy Mo) /H- 
From (8a, (82), (83), we have 
ro) 1 
Fo)=1F,, | oy3(vo— vd’ J 9¢0— 1!) de 


’ 


os) 1 d 
+ aF,, | o, dv’ S(v’, vey b's nr 
0 J0 pe 


rr 1s i ' a d ’ 
=F yon 4t+ >| ody’ | (V's Yo bo 


0 0 o,0,, KB 





by the properties of the 5-function. 
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If we put z=o in (88) and distinguish between upward and downward 
intensities, we have, by (90) and (91), 


d3, Ae el 


ac Fy(les+0,)S(% Yortsts) 2-H}, (03) 
[Ave ts =H) 2. = Fo, Fo). (94) 


Now by the principle of inv iii the law of diffuse reflection is unchanged 
if we omit a layer of the atmosphere of thickness z, and since no component of 
the field /,(z, 1) is being generated within the atmosphere itself, it follows from 


(21) that 
™ , , , S VY; v’, ’ ; , 
L+(e,u)= [av f e(s, wD ay’ (95) 
0 0 lad 
(87), (85) and the properties of the 5-function, this gives 
« 1 , , 
B,(2u)= [df Br(z, —p) ED gy 
0 0 lad 
I 7, —o(K o 
* ae F,,S(v, Vos ty Hole Art one , (96) 
On differentiating this with respect to z and putting z=0, we have, by (93) 
and (94), 


£L* y(n 0,)SC vanced 


pb 
pI(0) f , 
= yd S yb bar K, +o, F,S Vv, » ’ 
2u * eds’ {Sov sw!) Ee jan (Xn ») (¥, Vos Hs Ho) 


and this can be written 


. K,+0, Ky, 7 oy, 
}F re 9( Vy Vo» by Ho) { i + i } 
0 


S(v, ,B’)d 
=0, Ho) {1+ 5f" # dv’ f, eee. (97) 
0,0, Le 
Finally, on substituting for J(o) from (92), we get (22), after which the analysis 
of Section 4 follows unchanged. 
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Summary 

A general account is given of the development of the quartz clocks of 
the Greenwich Time Service, and the present installation is described. The 
methods by which clock performance may be studied are described and a 
convenient criterion is proposed. Quantitative results are presented in 
tabular form. The degree of stability now attained represents a significant 
advance upon that of some five years ago. Details are given of development 
work now in progress with a view to improving stability still further. 





1. General description.—The first quartz clock of the Observatory was 
installed at Greenwich in 1930, and was of the pattern developed at the 
National Physical Laboratory (1). This employed a ring crystal supported 
on phosphor-bronze pins and mounted in an evacuated glass vessel held at a 
constant temperature by means of a mercury thermometer. The crystal was 
incorporated in a simple Pierce oscillator circuit, and was designed to oscillate 
at 99°3 kc/s. The divider chain performed successive division by 3, 6 and I1, 
giving an output of 500 cycles per sidereal second to drive a phonic motor. This 
clock, Q3, took its place among the standard pendulum clocks of the Time 
Service, which were also supplemented by data received daily from the National 
Physical Laboratory regarding two similar quartz clocks Q2 and Q6. 

Clocks employing GT-cut plate crystals in a drive circuit based upon the 
Meacham bridge (2) were developed independently at the Post Office Radio 
Branch Laboratories. The oven and temperature-control arrangements were 
essentially a simplified version of the thermostat designed by Turner (3). The 
fundamental crystal frequency was 100 kc/s, which was converted by means of 
regenerative dividers (4) into a I kc/s output used to drive a phonic motor. 
A group of three such oscillators, known as ‘“‘Group 4’ (§) attained a high 
standard of performance, and in 1942 arrangements were made for information 
relating to these clocks to be communicated daily to the Royal Observatory, 
where it proved of the greatest value. Within a year it became clear that the 
long-term performance of these clocks was so markedly superior to even the best 
of the pendulum clocks available to the ‘Time Service that the pendulum clock 
was rendered virtually obsolete as a long-term standard (6). 

It was clearly desirable that similar quartz clocks should be installed at the 
Observatory, so that the Service would not be entirely dependent upon primary 
standards at an external establishment. ‘There was also an urgent need for 
secondary standards having good day-to-day stability for the control of the 
signals, and in this respect also the superiority of the quartz clock had been 
established. ‘The first three clocks supplied by the Post Office were of the 
Group 4 pattern, and were brought into use early in 1944 for the control of the 
GBR radio time signals, thus bringing about a considerable improvement in 
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the day-to-day uniformity of these time signals. Unfortunately, the value of 
these clocks as long-term standards was seriously impaired by the fact that they 
were mains operated and that, under war-time conditions, interruptions of the 
mains electricity supply were of frequent occurrence. It was thus still necessary 
to base the service on the Post Office Group 4 standards. 

During 1945 four groups of three oscillators, similar in general arrangements 
to the Post Office Group 4, were installed in new clock cellars at the Observatory 
station at Abinger (7). These oscillators, together with their ancillary equipment, 
were designed for battery operation, and an extensive installation of batteries 
and control equipment, together with a stand-by diesel generator, was provided. 
This overcame the main cause of disturbance to the oscillators, and longer 
uninterrupted clock runs were made possible. The temperature-control circuits, 
however, were still dependent on a 50 c/s a.c. supply, and in the event of mains 
failures this was provided automatically by a battery-driven motor alternator. 
‘Two similar groups of three oscillators each were subsequently installed at 
Greenwich, and incorporated some minor improvements. Special types of 
phonic motors were developed in order that all requirements of the service 
could be provided directly from the quartz clocks (8). Motors incorporating 
sidereal gearing were used to provide sidereal seconds for the astronomical time 
observations; and more elaborate timing devices, capable of producing the 
various forms of time signal distributed by the Department, made the’ full 
accuracy of the quartz clocks available for direct practical use. It had thus 
become possible for the quartz clocks to perform, with greater accuracy, all the 
functions hitherto carried out by the pendulum clocks: in fact, the manifest 
inferiority of the pendulum clocks led to their complete elimination. 

One result of the widespread employment of quartz clocks was that the clock 
comparison equipment which had been adequate for use with pendulum clocks 
was found to be incapable of meeting the new demands for accuracy. In 
particular, it became necessary to replace the electro-mechanical chronographs 
by fully electronic decimal counter chronometers (9). Moreover, quartz clocks 
may be conveniently intercompared with great accuracy by methods which are 
inapplicable to pendulum clocks. Relative rates may be readily obtained by 
means of a precision beat comparator (10) and integrated time differences 
between standards by means of rotary beat counters (11). It thus becomes 
possible to maintain a constant check on the performance of all oscillators and 
to connect the dividers and phonic motors to those oscillators giving the best 
performance at the time. 

As initially installed, the three oscillators of each group were mounted on 
one rack, in accordance with standard Post Office practice. Although oscillators 
employing GT-cut crystals are not unduly sensitive to mechanical disturbance, 
it was decided to re-mount the oscillators singly on individual brick and concrete 
piers, thus making it possible to carry out work on one oscillator without risk 
of disturbance to the others of the group. In the course of the necessary re-wiring, 
the power supplies were re-arranged to permit convenient isolation of each 
oscillator. New crystal ovens were constructed in the Observatory workshop, 
to an improved Post Office specification which is designed to accommodate 
either GT or ring crystals. The oven is located in a hollow in a concrete block 
on top of the pier, where it is protected from accidental disturbance and draughts. 
‘The temperature-control unit and maintaining amplifier are mounted indepen- 
dently above the oven. Some minor changes were made in the circuit 
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arrangements: high-stability components were employed, and the layout was 
modified to facilitate servicing. ‘The step-up ratio of the transformer in the 
temperature-control circuit was increased to improve the sensitivity, and the 
constant-temperature arms of the control bridge were re-wound on a solid brass 
rod. The period of the oven-control cycle was thereby reduced to about five 
seconds. Duplicate thyratrons were fitted, with separate signal lamps, thus giving 
warning of impending thyratron failure, and permitting replacement without 
disturbance of the temperature control. In the new maintaining amplifiers 
toroidal transformers are employed and precautions have been taken to render 
the amplifier gain still less sensitive to variations in supply voltages or in the 
loading of the 100 kc/s output. Some amplifiers have been fitted with prototype 
long-life valves. During the testing and setting-up of these oscillators it became 
necessary to devise new techniques for alignment and adjustment. 

In view of the excellent results achieved with a bridge drive circuit in 
association with GT-cut plate crystals, the Post Office Radio Branch Laboratories 
proceeded to produce a bridge circuit suitable for a ring crystal. Three 
oscillators, forming the Post Office Group 5, employed ring crystals supported 
by nickel pins. Later oscillators, forming Group 9, utilized softer metal for 
the crystal supports, and in the latest type the ring is suspended on silk threads. 
Very high values of Q (i.e. sharpness of resonance) have been obtained, but 
perhaps the most significant characteristic of the new ring oscillators is the 
substantial reduction of the initial frequency drift. An oscillator of this type 
was installed at the Observatory in 1950 May, and was brought into use almost 
immediately as the working standard E5. For subsequent employment as a 
long-term standard, its useful run began only a week or so after installation. 
A remarkable feature of this particular standard is the small value of the frequency 
drift, which amounts to only I part in 10° per year. A second oscillator of the 
same type, E6, was installed a few months later, and has exhibited a similar 
degree of stability, but is characterized by a uniformly changing frequency 
drift. ‘Two further ring-crystal oscillators were installed towards the end of 
1951. One of these, C5, is performing satisfactorily, but the other, C6, failed to 
reach the high standard of the other ring oscillators. It was accordingly removed 
from service, and the maintaining amplifier was modified in the light of more 
recent experience. The new run of this oscillator shows a satisfactory improve- 
ment in short-term stability, and there is good reason to suppose that the long-term 
performance will be correspondingly improved. 

The ancillary equipment associated with the oscillators has also received 
attention. ‘The regenerative dividers providing an output at I kc/s have been 
replaced by improved models developed at the Royal Observatory. The earlier 
pattern employed a large number of metal rectifiers and a total of fourteen 
transformers. In the new pattern, no metal rectifiers are employed, electronic 
mixing is used, and the number of transformers has been reduced to four. The 
valves were previously operated at nearly full rated anode current, and 
performance was adversely atlected by valve ageing. The valves are now 
operated at low anode current, with consequent saving of power and reduced 
sensitivity to changes in valve characteristics. In order to meet the need for 
accurately spaced seconds pulses, new pulse dividers have been designed. A series 
of identical plug-in units, each giving one output pulse for every ten input pulses, 
may be connected in cascade to perform frequency division by any required 
power of ten. The units function with equal reliability at any input frequency 





70 Humphry M. Smith Vol. 113 


from 100000 pulses per second to one pulse per second. A series consisting 
of five such units supplied with pulses derived from the 100 kc/s crystal frequency 
provides seconds pulses of more than adequate accuracy. Full advantage is 
thus taken of the extremely high stability of the latest crystal oscillators. 

2. Assessment of long-term performance.—For the operation of a National 
Time Service of high accuracy the clocks should provide a time scale that is 
uniform 

(a) from month to month over a period of two years or so; 
(6) from day to day over a period of one month or longer. 

The clocks which adequately fulfil the first requirement are employed for 
long-term prediction. ‘Those with a high degree of day-to-day stability are 
used as working standards in the control of time-signal transmissions. It is 
therefore necessary to assess the quality of the available clocks from these two 
different aspects. 

The problem of relating the performance of a modern quartz clock to the 
ultimate standard of time, the rotation of the Earth, has become increasingly 
difficult on account of the latest developments in quartz clocks. For some years 
the uniformity of the time scale provided by the clocks has surpassed that detined 
by the Earth: this disparity has been growing more pronounced since the 
introduction of improved clocks employing ring-crystal oscillators (12). In an 
attempt to remove the effects of one known cause of irregularity in observed 
time, it has been the practice for the last five years at the Royal Greenwich 
Observatory to apply current corrections for the variation in longitude arising 
from polar motion, or more briefly, “ polar variation”. ‘These corrections are 
based on independent data derived from variations of latitude observed at the 
U.S. Naval Observatory. Another known cause of irregularity is the annual 
fluctuation in the rate of rotation of the Earth. The mean amplitude of the 
fluctuation over a number of years was independently determined by several 
investigators. The degree of agreement between these various results suggested 
that the effect was sensibly constant from year to year, and that considerable 
advantage would be gained by the application of a correction based on the average 
fluctuation observed in preceding years. It has since become apparent that the 
annual fluctuation is subject to appreciable variation from year to year, and that 
the use of an average value may be misleading (13). This fluctuation can only 
be determined from the time observations themselves, and current estimations 
of its magnitude are subject to considerable uncertainty. It has in fact become 
impracticable to obtain a sufficiently uniform time scale by applying to the 
astronomical time observations provisional corrections for polar variation and 
annual fluctuation. A time scale of adequate uniformity may however be 
established by taking the mean of a number of specially selected quartz clocks of 
the highest precision. But such a scale is not absolute unless it can be related 
to the rotation of the Earth. It is therefore necessary to calibrate the scale with 
reference to the general trend of the astronomical observations over a period of 
two or three years. The procedure now adopted at the Royal Greenwich 
Observatory is to use a mean of the best clocks at the various cooperating 
establishments to define a time scale that is, as nearly as possible, uniform, and to 
determine the rate of the mean clock by consideration of the astronomical time 
observations made over an extended period at Greenwich and Abinger. 

The application in practice of these principles may be illustrated by reference 
to the analysis performed for the period 1950 June to 1952 June. Information 
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relating to all available clocks was readily accessible from the current graphs of 
clock intercomparisons. An examination of the performance of the clocks during 
the period under review permitted the selection of the five best clocks. These 
clocks were then referred to the astronomically observed time corrected for 
polar variation and the provisionally adopted annual fluctuation. An adequate 
basis was thus provided for the adoption, for each clock, of a preliminary parabolic 
ephemeris extending over the two-year interval. Where the intercomparisons gave 
clear evidence of a change of rate or drift in any particular clock, due allowance 
was made. In the course of the two years, two of the clocks, E5 and 9A, required 
no such adjustment; the remaining three each had one apparent change of 
drift: gC in 1951 April, Q13 in 1951 May, and EA in 1951 November. Each 
clock was then referred to the monthly means of the observations, corrected for 
polar variation but not for annual fluctuation. Corresponding monthly points 
in successive years were differenced to give a series of annual rates. If the 
fluctuation in the rate of rotation of the Earth were strictly repetitive from one 
year to the next, it would have no effect upon clock rates determined over these 
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twelve-monthly periods. In practice it is found that the fluctuation varies 
somewhat in amplitude and phase, but that this procedure is effective in removing 
the major part of the fluctuation in determined rate. Moreover, the use of a 
long interval considerably reduces the residual effect, and largely eliminates 
errors arising from observational scatter. The efficacy of the process is 
demonstrated by the fact that the plotted points adhere closely to straight lines 
(see Fig. 1), which were used to determine corrections to the parameters of the 
preliminary parabolic ephemerides. 

The five clocks were then combined into a mean clock, which established a 
provisional uniform time system of which the uniformity and scale value 
depended upon the drift and rate adopted for the mean clock. Since the drift 
and rate of each individual clock had been ultimately determined by consideration 
of annual rates, given by the observations, the time system defined by the mean 
clock was fundamentally related to the astronomical standard of time, from which 
the annual fluctuation had been removed. The system was however largely 
independent of the particular values adopted for the amplitude and phase of 
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the annual fluctuation. This time system is therefore suitable as a basis for the 
evaluation of the 12-month and 6-month periodicities of the annual fluctuation 
by a series of simple harmonic analyses, each covering twelve successive monthly 
means of observations corrected for polar variation. It has been found 
convenient to start each harmonic analysis at the beginning of a quarter-year 
so that any given analysis has an overlap of nine months with the previous and 
following analysis. There were therefore five harmonic analyses in the two-year 
period under consideration, the first and fifth being independent of each other, 
and the others overlapping. The five results for the annual fluctuation appeared 
to be in reasonable agreement, and a convenient set of values was adopted for 
the whole period. Corrections for annual fluctuation computed from these 
were then applied to the monthly means of observations, bringing the plotted 
points to the neighbourhood of a horizontal straight line. ‘The mean ordinate 
was computed and applied, fixing the zero-point of the uniform time system. 
The uniformity and scale value having been fixed as previously described, the time 
system was thus completely determined. The remaining scatter was attributed 
to observational error, and amounted to a mean deviation of + 6-2 milliseconds 
for a single monthly point. 
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In Fig. 2 are shown 

(A) the uncorrected observations (monthly means) referred to the mean 
clock ; 

(B) the correction for polar variation based on latitude observations made 
in the U.S.A.; 

(C) the observations corrected for polar variation, but still subject to 
annual fluctuation, together with the curve based upon the mean 
annual fluctuation adopted for the period; and 

(D) the observations corrected for polar variation and annual fluctuation, 
with the mean ordinate removed. 

It will be noticed that the polar variation correction, Fig. 2(B), shows a 
preponderance of positive values, arising from an apparent discrepancy in the 
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adopted mean latitude of the U.S. observing station. This is not of serious 
importance, as it affects only the zero-point of the uniform time system, and thus 
has no effect whatever upon the determined annual fluctuation, the determined 
astronomical time (G.M.T.), or the assessment of clock performance. An 
adjustment of the adopted mean latitude would give rise to a discontinuity in 
the uniform time system, and is thus better delayed until a revised mean can be 
adopted from, say, ten years of results. For the same reason, changes in the 
adopted position of the Earth’s mean pole, which are equivalent to adjustments 
of adopted mean latitude for many stations, are best made as infrequently as 
possible. It is of interest to note that during the period considered the effect 
of the polar variation on time observations at Greenwich was nearly equal in 
amplitude and opposite in phase to the effect of the annual fluctuation. Further- 
more, any errors in the adopted corrections for polar variation will be reflected 
in the determined annual fluctuation. It has been found that the corrections 
derived from the U.S. latitude observations are substantially confirmed by those 
based upon the results published subsequently by the International Latitude 
Bureau. 
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The establishment of a uniform time system has been described in some 
detail, since this is of fundamental importance in the assessment of long-term 
clock performance (Fig. 3). The mid-monthly values of the daily rates of 
typical clocks, in terms of this uniform time system, are listed in ‘Table 1. The 
stability in rate from month to month of the two best clocks, E5 and gA, is 
indicated in Fig. 4. The clock 9A has been in continuous operation since 1948 
January, and has evidently settled to a practically uniform frequency drift. The 
clock E5, on the other hand, was not installed until 1950 May, and appears to 
have settled immediately to a frequency drift which was not only practically 
uniform but also very small. 
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The most important requirement in a long-period standard clock is that 
it should run uniformly with a constant frequency drift for as long a period 
as possible. ‘Twelve months of uniform run are needed before even one 
determination of annual rate can be made. A uniform run of two years provides 
twelve monthly values of annual rate, and represents the shortest period that 
can give a reliable indication of frequency drift. For any period of approximately 
uniform run, the first differences of the monthly values of rate will represent 
the monthly drift in the rate of the clock. Fora perfect clock the first differences 
will be constant, and the second differences zero. The magnitude of the second 
difference will indicate the departure of the clock from uniformity of run. It 
was therefore decided to evaluate the mean absolute monthly second difference 
for each clock as a numerical criterion of clock performance. The appropriate 
criterion values, in milliseconds per day per month per month, are quoted at 
the foot of Table I. The requirements of a modern time service are such that a 
clock should maintain a criterion value not exceeding 0-10 for two years or so. 
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It may be noted that the criterion so determined is largely independent of 
errors in the adopted rate and drift of the reference time system. In fact, over 
an isolated period of two years, it may be regarded as independent of the 
astronomical determinations of time, since, for such periods, the mean of several 
modern quartz clocks gives a more uniform time scale than that defined by the 
rotation of the Earth. The astronomical observations relating the clocks to 
the rotation of the Earth are however essential both in current work and in 
retrospective assessment of performance. Even the best modern clocks are 
liable to slight changes of rate and drift which cannot be detected by clock 
intercomparisons until some time has elapsed and an appreciable error has 
accumulated. The astronomical observations, although subject to considerable 
scatter, serve to protect the currently adopted time from gross errors which 
might otherwise arise on this account. In retrospective assessment, the practice 
of judging the stability of quartz clocks in terms of their relative performance 
cannot be extended over periods greatly in excess of two years. In order to 
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control long-period changes in rate and drift, it is necessary to refer to an 
independent fundamental standard. The only practicable fundamental standard 
at present available for this purpose is the rate of rotation of the Earth, as 
indicated by observations of stellar transits. A more uniform fundamental 
standard of time would be provided by the period of revolution of the Earth 
around the Sun. 


TABLE I 


Mid-monthly values of daily rates of selected clocks 
(milliseconds per day) 


Abinger Greenwich Dollis Hill 
Standard E5 E6 gC EA 
Ring or Plate R R R R 
Silk, Pin or 
Clamp Sup- 5 5 P 
port 

1950 

July +47°83 —26:61 —0-32 
August +47°69 —27:03 —0-°76 
September +47°61 —27°47 —1°22 
October +47°58 —2785 —1°83 
November +47°44 —40°39 - — 28-26 -2°35 
December +47°26 —41°98 - 28-67 ‘69 
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1951 
January +47°24 —43°32 29°10 . —58-89 
February +47°25 —44°64 — 29°63 65 ~59°46 
March +47°Il —45°92 30°09 . — 60-09 
April +46°:95 —47°20 30°50 : — 60°69 
May +46-go -~ 48°50 -30°85 5° —61°14 
June +46°85 —50°04 -31°14 = — 61°66 


"54 : — 62°26 
‘92 — 62°77 
‘21 — 63°17 
66 — 63°72 
- 33°01 =F -64°15 
40 7°55 64-60 


July +46°84 
August +46°76 
September +46°51 
October + 46°34 
November +46°30 
December + 46°27 
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1952 
January + 46°23 
February + 46°16 
March + 46°05 
April +46°00 
May +45°92 62 
June +45°50 —76: +72°31 = 


Sr —789 —65-01 
28 -8:22 —65-40 
‘78 —860 —65°79 
‘23 —807 —66°15 
‘67 —9°28 —66°61 
‘21 —9°73 —67°18 
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Criterion 0-08 ; 0°53 0:06 0°07 0-08 : O15 


Observations made over many years have revealed that the Sun, Moon 
and planets are departing from their predicted positions by amounts that are 
proportional to their mean motions (14). This is attributed to an error in the 
time scale against which the observations have been recorded: the observed 
secular accelerations of the Sun, Mercury and Venus are thus associated with 
a progressive retardation in the rate of rotation of the Earth, amounting to an 
increase of one millisecond per century in the length of the day. It now appears 
that this may be satisfactorily accounted for by tidal friction. In addition there 
are fluctuations in the observed motions of the Sun, Moon, Mercury and Venus, 
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arising from variations in the rotation of the Earth. These variations have 
been interpreted as occasional sudden changes of several milliseconds in the 
length of the day, but may alternatively be regarded as the cumulative effects. 
of continual random changes. 

It is clearly desirable that the motions of the celestial bodies should be referred 
to a time system unaffected by irregularities in the rate of rotation of the Earth. 
Such a system is defined by the period of the revolution of the Earth in its 
orbit round the Sun, making the sidereal year the fundamental unit. At the 
Conference on the Fundamental Constants of Astronomy (Paris, 1950) it was 
recommended that the unit adopted should be the sidereal year at 1900-0 (15, 
16). The time reckoned in terms of this unit is known as Ephemeris Time and 
the difference between mean solar time and Ephemeris Time is given by the 
following expression :— 


+ 24%-349 + 728-3165 T + 298-949 T* + 15-821B, 


where 7 is measured in Julian centuries from 1900 January 0, Greenwich Mean 
Noon, and B is the fluctuation in the Moon’s longitude, that is, the difference 
between the observed longitude and the amended tabular longitude (14). The 
previously accepted unit of time, the mean solar second, which is now known 
to be variable, may be converted to an invariable unit, the second of Ephemeris 
Time, by use of the same formula. 

The second is one of the three fundamental physical units, and thus enters. 
into many scientific measurements besides those of astronomy. In the majority 
of cases, the errors arising from the use of the variable mean solar second are 
negligible, but clearly an invariable unit of time should be employed. ‘The 
second of Ephemeris Time is specifically intended to provide a unit of this kind, 
but its widespread use is dependent upon the development of satisfactory methods 
for its current determination. ‘This requires regular accurate observations of 
the position of the Moon for comparison with the tabular places. Work now in 
progress at the U.S. Naval Observatory suggests that observations of sufficient 
accuracy may be secured by the employment of a new photographic technique. 
These developments may make it possible for a time-keeping observatory to 
establish a uniform time system, and thus provide a satisfactory basis for the 
assessment of clock performance and for frequency standardization. 

3. Assessment of short-term performance.—It will be evident that the day-to-day 
stability of quartz clocks over periods of a month or so may be satisfactorily 
assessed from the clock intercomparisons alone without reference to the 
astronomical observations. As may be expected, the clocks which have been 
selected in virtue of their high long-term stability also exhibit a high standard 
of uniformity from day to day. An example is afforded by the relative daily 
rate of the clocks 9A and 9C, which is plotted in Fig. 5. ‘Throughout the year 1951 
the variation in rate was smooth, and the range did not exceed 0-5 millisecond 
per day. It may be remarked that the drifts of the clocks are sufficiently similar 
for the relative rate to be plotted direct, without removal of an ephemeris. 

The most satisfactory indication of day-to-day performance is provided by 
an examination of the second differences of the daily values of clock rate. ‘The 
mean absolute value of the second difference affords a quantitative criterion of 
the suitability of the clock for day-to-day extrapolation and interpolation. It has 
recently been pointed out (17) that there is some correlation between the second 
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differences on successive days. It should not be assumed, however, that this 
necessarily impairs their utility as a measure of short-term performance, since it 
has been shown (6) that they may be grouped into three sub-sequences, each of 
which is entirely free from internal correlation. The re-combination of the 
sub-sequences does not invalidate the result. In practice, it is unnecessary to 
separate the sub-sequences and re-combine when forming the mean absolute 
second difference. 
195! 
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Clock intercomparisons of high accuracy are obtainable from the daily 
readings of the rotary beat counters, which give a continuous indication of the 
time difference between any two clocks in units of one microsecond. The 
calculation of the second differences of daily rates is readily performed by means 
of a National Accounting Machine. These daily values are plotted to give a 
convenient picture of the day-to-day stability of the clocks. The mean absolute 
value of the second difference of daily rate provides a useful criterion of short-term 
performance. The intercomparisons between the two clocks E5 and E6 
throughout the year 1951 gave a mean absolute daily second difference of relative 
rate amounting to 16 microseconds per day. Since the two clocks were of 
comparable short-term stability, the criterion value was approximately 11 (in the 
same units) for each clock. The intercomparisons between C5 and C6 for a 
period of some four or five months after the modification of the C6 maintaining 
amplifier gave a mean absolute daily second difference of 19 microseconds per 
day, corresponding to a criterion value of 13 for each clock, if the stabilities were 
equal. It is clear that a short-period criterion value of 10 to 15 is attained by 
modern quartz clocks. While this accuracy cannot be fully utilized in Time 
Service operational work, the precision of the daily intercomparisons provides 
a valuable indication of any slight deterioration in the stability of the clocks 
which might, by its cumulative effect, cause a significant deterioration in the 
long-term performance. It may sometimes be possible to associate an 
enhancement of these small variations with particular physical causes, such as 
alterations in the values of certain components, thus providing information of 
great value in the modification and design of quartz-crystal clocks. 

A graph of the direct intercomparisons of E5/E6 and C6/E6 with a parabolic 
ephemeris removed in each case is given in Fig. 6. It will be noted that the 
deviations do not exceed 0-1 millisecond in the course of a month. 

4. Present developments.—The high stability of the performance attained by 
modern clocks leads naturally to the consideration of the possibility of further 
progress. It has been suggested that the ultimate limit of accuracy is imposed 
by the crystal itself, rather than by the ancillary circuits. If this were indeed the 
case, consideration would have to be given to the development of standards 
utilizing some control element other than a quartz crystal. Tests on the ancillary 
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circuits have shown, however, that in many cases the crystal is not operating 
under the best conditions. It has in fact been possible, by improved design of 
the associated circuits, to obtain a higher standard of stability from crystals of 
conventional types. 

In particular, attention was directed to the reduction of the amplitude of 
oscillation of the crystal. This cannot be achieved in the conventional bridge-drive 
circuit, where the amplitude of oscillation is controlled directly by a tungsten 
filament lamp. As a first step, the tungsten filament lamp was replaced by an 
indirectly heated thermistor, for which the heating current was derived from a 
secondary amplifier. This gave some reason to believe that higher short-term 
stability would result from a lower amplitude of oscillation. In the most recent 
drive-circuit, the bridge has been replaced by a circuit which is electrically 
equivalent, and in which the crystal current is less than 30 micro-amperes. The 
maintaining amplifier uses no transformers, as these were thought to be a source 
of instability in the circuits previously employed. ‘The achievement of higher 
precision in the crystal oscillator makes it advantageous to control the temperature 
of the crystal within closer limits and it is estimated that a new control circuit 
and oven divide ambient changes by a factor of 500 
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Reliability of operation over long periods is of paramount importance, since 
the clock may often suffer greater disturbance from a failure of temperature 
control than from a short interruption of oscillation. The previous designs used 
the 50c's mains supply both for the control bridge and for the heater windings. 
In the new experimental oscillators the current is supplied by a battery-operated 
325¢,s oscillator which is amplitude-stabilized. The oscillator valves are in 
parallel, and removal of either valve has negligible effect on the frequency and 
amplitude. The controlling amplifier, which employs duplicated hard valves in 
place of thyratrons, has been designed in such a way that the valves may be 
replaced without altering the control-temperature by more than half a millidegree. 
The first oscillator incorporating these refinements employed a GT-cut plate 
crystal, and was put into operation in 1951 November. A second oscillator, 
incorporating minor improvements, was installed in 1952 July. Results obtained 
so far have been encouraging. 

5. Conclusions.—It is now nearly ten years since the quartz clock, by virtue 
of its superiority as a long-term standard, displaced the pendulum clock in the 
fundamental work of the Greenwich Time Service. Almost at once it became 
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appatent that among the advantages of the new type of clock were a considerably 
higher standard of short-period stability and far better facilities for accurate 
intercomparisons. The purpose of this account is to show the further advances 
which have subsequently been made in the quality of the clocks themselves and 
in the technique of using them in the operation of a national time service. 

The past years have afforded an unrivalled opportunity for the gaining of 
experience and the investigation of the problems concerning quartz clocks. 
In those observatories which have not yet acquired quartz clocks comparable in 
stability with those available to the Greenwich Time Service, there is perhaps 
inadequate realization of the standards now attainable, combined with an 
instinctive reluctance to forsake the relatively simple pendulum clock in favour 
of the admittedly more complex quartz-controlled oscillator. The pendulum 
clock has the valuable property of being capable of prolonged uninterrupted runs, 
and it has been claimed that this feature is sufficient to justify the retention of 
pendulum clocks as fundamental standards, and that the long-term performance 
of the two types of clock is otherwise comparable (18). 

Considering first the question of length of uninterrupted run, it has been 
found that modern quartz clocks have attained a high standard of reliability. 
Given continuity of power supplies, the period of uninterrupted run may easily 
extend over several years. Among the clocks mentioned in this account, Fr has 
run continuously since its installation in 1947 July, 9C since 1947 October, and 
gA since 1948 January. The clocks at Abinger have not run for such long 
periods, as the installation has been modified in recent years. It is worthy of 
note, however, that in the six years since the Abinger clocks were transferred to 
battery operation there has been only one failure of a crystal oscillator. 

As regards performance, the superiority of the quartz clock is firmly 
established. It is now some years since pendulum clocks were employed as 
standards in the Greenwich Time Service; but up-to-date information on the 
performance of high-precision pendulum clocks at Paris is published regularly 
in the Bulletin Horaire. Mean daily rates are quoted for each month, and these 
may be compared with the corresponding figures for quartz clocks given in 
‘lable I. It is also possible to calculate the mean absolute second differences of 
the monthly rates, thus giving a numerical criterion of performance which may be 
compared with the values deduced for quartz clocks. During the two years 
1950-1951 one of the Paris pendulum clocks, Leroy 1372, attained a standard 
which can have been only rarely equalled by any other pendulum clock. The 
mean absolute monthly second difference was found to be 3-4 milliseconds per 
day. An examination of the other pendulum clocks shows that a value of 10-0 
would represent a good average performance, while a criterion of 50-0 or over 
would be classed as poor. Results are also quoted for the quartz clocks at Paris. 
The best of these has a criterion of 8-5, whereas the other two have values of 
about 20-0. It will be evident that, despite recent claims to the contrary (19), 
even the best pendulum clock does not represent a serious challenge to the typical 
modern quartz clock with a criterion of 0-1 millisecond per day per month per 
month. 

The high short-term stability of a quartz clock has permitted a new standard 
of accuracy in the control of radio time signals, and the application of frequency 
measuring technique has made it possible to carry out clock comparisons with a 
speed and precision hitherto unattainable. But it is perhaps in the realm of 
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long-term performance that the greatest advances have been made. Qtartz 
clocks were first introduced in the Greenwich Time Service as long-term standards, 
and it is due to the high stability of the present clocks over periods of two or three 
years that it has become possible to determine the nature of the annual fluctuation 
in the rotation of the Earth. In current practice the annual fluctuation is 
determined by referring the astronomical observations to each of a number of 
selected clocks. ‘The various solutions show a high degree of agreement, and the 
resulting mean value can be adopted with confidence. 
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value in the operation of the National Time Service. The author wishes to take 
this opportunity of expressing his appreciation of the continued cooperation of 
the Post Office Radio Experimental and Development Branch Laboratories and 
of the National Physical Laboratory. 


Royal Greenwich Observatory, 
Abinger Common, 
Dorking, Surrey : 
1952 October 9. 
References 


(1) L. Essen, Proc. Phys. Soc., 50, 413, 1938. 
(2) L. A. Meacham, Bell System Technical Journal, 17, 574, 1938. 
(3) L. B. Turner, Journal I.E.E., 81, 399, 1937. 
(4) R. L. Miller, Proc. I.R.E., 27, 446, 1939. 
(5) C. F. Booth and F. J. M. Laver, Journal I.E.E., 93, Part III, 223, 1946. 
(6) W. M. H. Greaves and L. S. T. Symms, V.N., 103, 196, 1943. 
(7) H. M. Smith, Annales Frangaises de Chronométrie, 2° Série, 3, 227, 1949. 
(8) H. M. Smith, British Science News, 2, 169, 1949. 
(9) H. M. Smith, British Science News, 2, 207, 1949. 
{10) H. B. Law, Journal I.E.E., 94, Part III, 38, 1947. 
{11) W. A. Marrison, Bell System Technical Journal, 27, 510, 1948. 
(12) H. M. Smith, Proceedings I.E.E., 98, Part II, 143, 1951. 
(13) H. M. Smith, 1.E.E. Monograph, No. 39. 
(14) H. Spencer Jones, .W.N., 99, 541, 1939. 
{15) Bulletin Astronomique, tome XV, fascicule 4, 291, 1950. 
(16) D. H. Sadler, W.N., 111, 624, 1951. 
(17) N. Stoyko, Annales Frangaises de Chronométrie, 2™° Série, 3, 139, 1949. 
(18) R. Forga, l’Astronomie, 65° année, 112, 1951. 
(19) A. Danjon, Annales Frangaises de Chronométrie, 2° Série, 3, 365, 1949. 





ON THE MASSES OF SATURN’S SATELLITES 
Sir Harold Jeffreys, F.R.S. 
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Summary 


The principal determinations of the masses of Saturn’s satellites are 
rediscussed, with special attention to the standard errors. A new 
determination of the mass of Titan is made from the motion of the orbital 
plane of Iapetus, and an attempt is made to estimate the errors remaining 
in Woltjer’s solution from Hyperion. The values resulting are as follows, 
expressed as fractions of the mass of Saturn: 

Mimas (6:69+ 0:20)10-*; Enceladus (1:27+0°53)10~7; 
Tethys (1-141 + 0°030)10°°; Dione (1-825 + 0°061)107-°; 
Rhea (0°4+ 3'8)10~*; Titan (2-412+0'018)107?; 
and Struve’s constants for Saturn are 
k ay? =0°024 305 + 0000 048, 
l/aogt=(11 0+ 2°0)10~. 





1. The mass of ‘Titan was determined almost simultaneously by S. Newcomb 
from the motion of the apse of Hyperion and by H. Struve from a comparison 
of the motions of the orbital plane of Iapetus and the apse of Titan. The results 
differed by a factor of nearly 3, but successive corrections to the theory of 
Hyperion have produced a rough agreement. There seems, however, to have 
been no new solution from Iapetus. Struve’s discussion* is statistically very 
unsatisfactory. It is repeated with little alteration by Tisserand.+ Using the 
same data, I applied the method of least squares, supposing that standard 
errors in 7 and @ (longitude of the node) were in the ratio sinz to I, and giving 
equal weight to all data. The result was (the mass of Saturn being taken as 1) 


My; = (2°324 + 0°075)1074, 
I/my, =4 308 + 138 (s.e.). 


Struve’s value was given as 4678 + 350. The standard error suggests that with 
modern data the method might yield an uncertainty under I per cent. 

J. Woltjer{ gives seven determinations from Hyperion, ranging from 
3767 to 4172, and adopts the value 4033. 

H. Struve§ discussed the observations of all the satellites and made important 
contributions to the theory. G. Struve| rediscussed the observations, with 
others up to 1928, and derived new masses. Estimates of uncertainty are given 


*H. Struve, Beobachtungen d. Saturnstrabanten, Supp. 1 aux Observations de Poulkova, 1888. 
+ F. Tisserand, Théorie de Mécanique Céleste, 4, 98-102, 1896. 

t J. Woltjer, Ann. d. Sternwarte Leiden, 16, Part 3, 1928. 

§ H. Struve, Publ. Obs. Centr. Nicholas, 11, 1898. 

|| G. Struve, Veréff. Berlin-Babelsberg, 6, Hefte 1-5, 1933. 
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for the equations of condition, but the final masses and their uncertainties are 
based on judgment and not on least squares solutions. Accordingly it cannot 
be said that their uncertainties are known, which is unfortunate in view of their 
importance for an understanding of the constitutions of these bodies. In the 
course of the present work it was found that G. Struve’s uncertainties ranged 
from half to twice the s andard errors. 

2. The following analysis is mainly based on G. Struve’s data. Some of the 
equations of condition were rediscussed but no important discrepancy was 
found. Uncertainties have been converted to standard errors. 

The motions of the nodes and apses derived from the secular terms in the 
disturbing functions due to the ellipticity of Saturn and the attractions of the 
satellites are denoted by +m; they are theoretically equal and opposite, to the 
accuracy available. Approximate commensurability of the mean motions 
produces systematic changes of the secular motions for the pairs Enceladus 
Dione and Titan—Hyperion. ‘The solar effect has been subtracted from all 
observed values. The theoretical values as functions of the masses are given by 
G. Struve (Heft 4, p. 53) for the first six satellites. A discussion of the 
observations of Titan and lapetus is given in the posthumous Heft 5. 

For the pairs Mimas—Tethys and Enceladus—Dione, though the masses are 
small, good determinations are possible because fairly large perturbations are 
produced by approximate commensurability of the mean motions. 

The mean motions in degrees per day given by G. Struve will be adopted 
(Table I). The value for Hyperion is Woltjer’s. 


TABLE I 


Mimas 381°994 442 +0°000013 
Enceladus 262-731 940 5 +0:000 002 
Tethys 190°697 950 +0°000 005 
Dione 131°534.972 9 +0°'000 OO! 2 
Rhea 79690 088 1 +0°000 000 7 
Titan 22°577 O15 08 + 0000 000 40 
Hyperion 16-919 989 6 +0:000 002 6 
Iapetus 4°537 995 36+.0-°000 000 40 


The observed values of Bm, in degrees per year, are in Table II. 


TABLe II 

wa 6 
Mimas 36560 +0°10 36523 to-1o (Heft 4, pp. 18, 47) 
Enceladus orn Pee 
Tethys nas 72°227+0°'031 (p. 42) 
Dione . . aes (p. 33) 
Rhea Lz 10°20 +004 (p. 40) 
Titan . + O° 0°492+ 0°029 


The mean motion of the apse of Titan is given (Heft 5, pp. 12-13) as 
31-41 +021 per year; from this must be subtracted the solar part 1-34 
(H. Struve, p. 172), leaving the value given. The value for the node is a rough 
determination from the data in G. Struve’s table (Heft 5, p. 11). The values 
for Iapetus need further discussion. 

I consider first the periodic disturbances. 
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2.1. Mimas and Tethys.—The theory is outlined by Tisserand (pp. 129-138). 
G. Struve (Heft 4, p. 10) writes the perturbations of mean longitude of Mimas 
in the form 
Al=A sin p(t —7)) + B sin 3u(7 — 79) 

and determines 

A = — 44°:390 + 0°122, 

p= (5°:086 4 + 0°-008 5)/year, 

T, = 1 866-27 + 0-06. 
For Tethys A is replaced by A, and the same values of yu, tg are adopted; then 

A, =2°:065 +0°-028. 


The ratio my,/my, then follows from 


at __P 4 =0°093 04. 


ame 
The value of my, depends on the use of elliptic functions, but substantially 
= 
Heyyy 


where y, y, are the inclinations of the two orbital planes to Saturn’s equator. 
Allowing for the uncertainties of these and of yu gives a proportional standard 
error of + 2 x 0-001 7 + 0:025 +.0:009 = +0°026. Then 


My.X ’ 


Mr. = (I +.0°026) = 1-141 x 10-*(1 + 0-026). 


z 
576 400 


The relative uncertainty of my; is 0-014 from that of A, ; that of A is unimportant, 
but those of yu, y, y, enter again; then 


My = (I + 0°026 + 0-014) 


14 960 000 
=6-69 x 10-8(1 + 0°030). 


G. Struve gives 837000 and 915000 as extreme values of I my,, corresponding 
to about 1-5 times the standard error. On the other hand, he gives the uncertainty 
of the mass of Mimas as 5—10 per cent, the latter being too high on a similar scale. 

2.2. Enceladus and Dione.—The largeness of the perturbations in this case 
also depends on the fact that the mean motions are nearly as 2 to I, but the 
theory is much more difficult. In the case of Mimas and Tethys the argument 
of the libration is 4/,—2/—@-—6,, and thus depends on a product of the 
inclinations and not on the eccentricities. For Enceladus and Dione it does 
depend on the eccentricities. H. Struve used a linearized theory, which is 
used in part by G. Struve. Woltjer* pointed out that higher terms need 
attention, and G. Struve quotes his solution also. It allows, however, only for 
perturbations of Enceladus by Dione, and a somewhat fuller theory is needed. 
The following method is that of H. Struve extended to allow for certain 
second-order terms. Letters with suffix I refer to Dione, without them to 
Enceladus. 

* J. Woltjer, Bull. Astr. Inst. Neth., 1, 23-24 and 36, 1922. 
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The rates of change of the relevant orbital elements are given by (H. Struve, 
1898, p. 178) 
dn 


* = 3m,n"{ Ae sin (2/, —1—w) — Be, sin (2/, —-/—w,)}, (1) 


le ‘ 
= =mnA sin (2l,—l—w), 
cd 


[= = Pn—m,n = cos (2/,-—l/-m), 

dn, ° : . 

> he —6mn,*{A,e sin (21, —1—a) — Bye, sin (2/, -—1—a,)}, 
de, 
at 

dw, 


dt 
A=0°7532, B=0:2714, A,=1°1943, B,=0°4348, AB,—A,B=0-003 4. 
Put (7) 
V5 =(2n,—n)t, V =2l,-l=V,+V’, (8) 
h=esin(w—V,), k=ecos(w—V,), h=e,sin(w,—V>), 
k=e,cos(w,—Vo), (9) 
v=Pn—(2n,—n), v,=2n,—n— fn. (10) 
Then 
dh * 
_ =vk—m,nA cos V’, (11) 
. = —vh+mnA sin V’, 
at 
th 
= = —v,k,+mn,B, cos V’, (13) 
dk, 
dt 
d*l_ dn , ? lie on ° 
> =3m,n>(— Ahcos V’ + Aksin V’ + Bh, cos V’ — Bk, sin V’), (15) 
dl, dn, ° F ’ - a 7) ° 7) - 
> dale tes —6mn,?(— A,hcos V’ + A,ksin V’ + B,h, cos V’ — Bk, sin V’), (16) 
d*V’ 
dt* 


= —mn,B, sin (2l,—1—m,), 


B ’ 
= Bn, +mn, @, 008 (2h — 1-2), (6) 


=v,h,—mn,B, sin V’, (14) 


= P(h cos V’—k sin V’)— Q(A, cos V’ —k, sin V’), (17) 


where 
P =12mn,?A,+3m,n°A, (18) 
O =12mn,*B, +3m,n*B=BP/A. (19) 
For any constant «, y 
V'=a+yt 


leads to an exact solution of (11) to (14) and (17). This corresponds to the 
addition of constants to n, m, and the epochs. But it also indicates that the period 
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equation contains y* as a factor, and that we can choose the constants so that V’ 
has no constant or secular part. This identifies , m, with the observed mean 
motions. 
For V’ small, if we neglect V*, there is a constant solution 
h=h,=0, V’=o0, 


m nA 


k=hy= =, ky =k = - 


72 


—. (20) 


V1 


Accents denoting departures from this constant solution, we put 


h=csinyt, k'= Y ¢-cos yt, hy=c,sinyt, k,’=—- Las cos yt, V’=Asin yt 
Vv Vy 


(21) 
and retain the first powers of c, c,, A. Then, cancelling a factor y"A, we have the 
period equation 

Pky __QRw 
fF P-y 
. Ava Byv,A 
with c= — ¢,= =e 
; alee Y 

To the second order in ko, yg the roots are given by 

k 

y?—v? = Ph, (: ~ ow.) (23) 


v2 —v2= —Ohyo (: - aa :) ‘ (24) 


> 
iliac 


=I (22) 





2 
2 


> @ 
_ Vy 


and 


({ Ak Bk Ak Bk ‘ 
Al= — 3m,n? {( — —- —— :)a sin yt + ( i 10 ja sin tt, 2« 
z : | Y —yp" a — v\" Y vi" oy p2 yy? ae v" 1 V1 ( 5) 


Al, =6mn,* J (* on ae sin yt + ( Aiko — ByRio ) A, sin nth, 


2 2 7 >. an = 
<2 rT a Y1 t ye? 


_  2mn,* Ay.) mM 
~  myn* A ad al. (26) 














k, is small and ky, still smaller, and if we neglect m the solution reduces to that of 
Woltjer. Note that A/, A/, are not small compared with A, A,, since y*—v* and 
y," — v,° contain the small factors m, m,. 

G. Struve gives /, e and w for mean dates, each summarizing observations 
over about a year. With e anda he gives also 2/, —/, which never differs from w 
by more than a few degrees. Thus k from (9) is practically the recorded e, and 
the mean value of the latter is an estimate of ky. Thus hp is the forced eccentricity 
of Enceladus. The variations of e and of esin(2/,—/—m) are however too 
small to give useful determinations of the periodic parts of h and k. The values 
of e give 

k, =0°004 48 + 0-000 21. (27) 
G. Struve’s value (p. 33) is 0-004 44 + 0-000 09 (p.e.) but the uncertainty seems 
to be a mistake. 

For Dione the mean eccentricity given by Struve’s table (p. 34) is 


&, =0°002 II + 000009. (28) 
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Struve gives 0:0022I without an uncertainty. (Some of these small changes 
are due to differences in the allocation of weights where two sets of observations 
correspond to one year. I have checked most of the summary values by means 
of x”, and took the larger uncertainty in cases of disagreement, which seldom 
occurred.) As for Enceladus the periodic variations h,’, k,’ are too small to be 
detected from the variations of e and w. 

The periodic parts of / are found to be (G. Struve, p. 23) 


Al = (14':39 + 1':20) sin {63°-75 + 5°:28 + (32°-51 +0°-30)t} 
+(14'-06 + 1'-19) sin {117°-28 + 5°47 + (93°*14 + 0°40)t}, (29) 
and direct comparison of this with /, gives 


Al 
— —1 =x= +0-065 + 0-025. (30) 
Al : 

The speeds of the harmonic terms are estimates of y, y,, but for their 
interpretation we need also preliminary estimates of fn, f,n,. The free 
eccentricity of Dione (é,) leads observationally to 

do, /dt = 30°75 + 0°21; (31) 
there is no direct evidence at all for Enceladus. H. and G. Struve proceeded 
by interpolation from the values for Mimas and Tethys. I extend their method 
by fitting a form Cn’? + Dn- to all the observed values out to Rhea, giving equal 
weights to each. This amounts roughly to regarding the secular motions as 
mainly due to the ellipticity of Saturn and perturbations by Titan, with a small 
part due to roughness of the Titan term and to the effects of the other satellites, 
which is taken as random. The object is only to get preliminary estimates of 
Bn, B,n,, which can be improved later. The result is, with 8m in degrees per year, 
n in degrees per day, 

C=3°45008 x 10-4, D=+53, 
with the following results. 

O-C 
Mimas 3605-26 + 0°14 = 365°40 — 0°17, +0°20 
Enceladus 152°57 + 0°20 =152°77 
Tethys 72°23+0:28= 72°51 — 0-28 
Dione 30°35 +0°40= 30°75 0:00 
Rhea 9°43 +0°67 = 10°10 +010 


The estimated standard error of one value is 0:23. This is more than the 
apparent uncertainties of the observed values, and it will be worth while to 
attempt greater detail. It is clear, however, that the two observed values for 
Mimas differ by 2-5 times the standard error of their difference, and that the 
interpolation gives Sn for Enceladus with as much accuracy as the speeds of 
the two harmonic terms already possess. For Dione the interpolate happens 
to agree exactly with the observed value. I take 


Bn =152°-77/1* + 8(Bn), (32) 
Byn,= 30°°75/1°+8(8,n,), (33) 


where the corrections are unlikely to exceed 0°:3/1°. Then 


2n,—nN=123°°44/1°, v=29°-33—8(Bn), v,=92°-69+8(A,m,). (34) 
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I take also 
m=1:5x10-(I+p), m,=1°8 x 10-%(1+4)). (35) 
Then 

ky = 0-004 46 (I + 4,) (I + 355(8n)), (36) 
hy = — 0-000 0260(1 +p) (1-5 8(B, m,)), (37) 

Pk a 
4 = (2:056 + 0-238 + 3°870 p,) 10-8 (I + 355(Bn)), (38) 

Ohy a < a) 8 —11 is 
<* = — (4°32 + 4°82 4+ 3°81 p,) 10-" (1 — 5 8(Bn,)), (39) 


n 
whence, with sufficient accuracy, 


¥ = (1-105 7 + 0-0106 p+ 0-172 34) (I — 35 5(Bn)), (40) 
7 = 32°43 + 0°34 + 5°59 H1 — 8 (Bn), (41) 
and 
0°34 #+5°59 Hy = +0°08 + 0:30 +8 (Bn). (42) 
From ky 
0004 46 (4, + 355(Bm)) = + 0-000 02 + 0:000 21. (43) 
Also 
(y1° — 4")/v,7==10-% (44) 
and is negligible. ‘The ratio A/,/A/ leads to 
0-066 I 4 — 0-066 I uw, = — 0-001 I + 0°025. (45) 
The mean eccentricity of Dione is c, = é, ; the amplitude of the mode with speed y, 
for Enceladus is Struve’s q, 


n\? ’ , , 
iit di, (*) =5 459’ ¢, (1 +444) =14'06 + 1-20. (46) 


In c, uz, we take c, from (28); then 
C, + 0°002 I pp, = 0-002 58 + 0-000 22. (47) 
5(8n) and 8(8, ,) have been neglected where their effects are negligible in compari- 
son with the uncertainties already present. On account of (44) the “ periodic’”’ 
term in the disturbing function due to Enceladus produces a negligible effect in the 
motion of the apse of Dione. 
Then the equations of condition for j, 4, ¢, are (42), (43), (28), (45), (47). 
The solution is 
= —0°027 + 0°37 8 (Bn) + 0°35, (48) 
= +0°030 + 0047 8(Bn) + 0-034, (49) 
€, =0°002 19 — 0-000 02 8(fn) + 0-000 08. (50) 
Comparing with the equations of condition, and ignoring 5(f), we get x* =3°2 on 
2 degrees of freedom. Then 
m = 1°46 x 1077 {1 +. 0-388 (Bn) + 0°36}, (51) 
m, = 1°854 x 10° {1 + 0:0468 (Bm) + 0-033}, (52) 
1/m=6 850000, 1/m,==539 400. (53) 
G. Struve gives 6 622 000 and 541 300, with uncertainties of 20 and 5 per cent. 
That for Enceladus is too low; that for Dione is about 1-5 times the standard error. 
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The solution is not quite so satisfactory as it appears, for a reason given by 
G. Struve (p. 37). The argument of the y, term is essentially 2/, -/—m,, and the 
constant term in this is (with probable errors) 
— 173°°4 + 3°°4+ (507° 44’-0 + 0'8) — (199° 25'-8 + 0'-8) = 134°-90 + 3°4. 
The observed value is 117°-28 + 3°-69, and the difference is 17°62 + 5°:15 (p.e.) 
or +7°:62(s.e.). This would have contributed 5-4 to x” had it been included in the 
comparison, making x” = 8-6 on 3 degrees of freedom. It did not arise in the actual 
solution because the phase affects none of the equations of condition, but there are 
also certain persistences of sign in the residuals that suggest unconsidered long- 
period perturbations. If these have been inadequately eliminated in the method 
of solution the uncertainties may need to be increased by a factor of 1-5 or so. 
2.3. Rhea.—This has a forced eccentricity due to Titan, and the pericentre 
oscillates about that of Titan. G. Struve (Heft 1, p. 11) adopts my,= 74; and 
derives the libration constants by least squares. The residuals of e after allowance 
for the libration give a mean (—2+6-7)10-°; the adopted mean was 0-000938. 
‘Taking the mean eccentricity as proportional to my, we get 
g6 + 6-7 P 
gd 4033 


I 
— =4117+ 288. 
mr, 4117+ 





M7, = =0°000 242 9(I +0°070), 


2.4. Titan.—For the mean motion of the node on the proper plane G. Struve 
adopted H. Struve’s value of 0°-506/year, and he apparently did not try to improve 
it. Forthe apse and node I take the values in Table II. 

An estimate of the mass of Iapetus is derived from the inclination of the orbit 
to Saturn’s equator, namely 

Ms, = 2°46 (1 +.0°78) x 10-8, 
The uncertainty has as usual been converted into a standard error. ‘The result is 
only an indication that the mass is comparable with those of Tethys and Dione. 

Examination of possible effects of Hyperion shows that its mass is less than 
smo Of that of Titan. 

2.5. Hyperion.—The theory of this satellite is far the most complicated of all. 
Woltjer’s treatment* is the fullest. ‘The essential points may be indicated by 
comparison with the perturbations of Enceladus by Dione. The chief differences 
are that the mean motions are nearly in the ratio 4 to 3 instead of 2to1; the reaction 
of Hyperion on Titan can be neglected; and the mass of Titan is about 100 times 
that of Dione, so that approximations have to be carried to a much higher order. 
As the counterpart of kj, is negligible, that of y, reduces to v,. 

The resemblances are as follows. (1) The ky term in the motion of Enceladus 
has a much larger counterpart. ‘This is a forced eccentricity of about 0-1 combined 
with a motion of pericentre 4n—3n’, accents referring to Titan. (2) The free 
oscillations about the resulting periodic orbit, in the case of Enceladus, are 
represented by the term in mean longitude of speed y. When the departures from 
a circular orbit are described in terms of an eccentricity and a moving apse, and 
the free eccentricity is less than the forced one, the apse can be described as librating 
about a point with constant motion. Hyperion shows a similar feature. (3) The 
y, term in the mean longitude of Enceladus is a forced term due to the eccentricity 


* J. Woltjer, Ann. d. Sternwarte Leiden, 16, Part 3, 1928. See also Bull. Astr. Inst. Neth., 1, 
13-16 and 31-36, 1922. 
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of Dione. The eccentricity of Titan produces a corresponding disturbance of 
Hyperion. 

The ditferences are quantitative matters and are shown chiefly by the slowness 
of convergence of series for Hyperion. The disturbances are not so far from 
linearity that a partly linearized theory similar to that for Enceladus would be 
seriously inerror. Inspection of Woltjer’s series shows that higher terms, though 
not negligible, are much smaller than the main terms. In these conditions it 
appears that further corrections can be treated by adapting the theory of Enceladus, 
the chief features having already been treated fully by Woltjer. 

Woltjer (p. 66) quotes the following estimates of 1/my,;=1/m' found from 
Hyperion. (a), (b), (c) ate his own. 

(a) Mean motion of libration argument 3.986 

(b) Mean motion of pericentre 4080 

(c) Mean motion of node 3 767 +140 
(d) Hill’s method (computed by D. Brouwer) 4143 

(e) Eichelberger (special perturbations) 4172+58 
(f) Longitude of pericentre (Samter) 4125 

(g) Libration (Samter) 3910 


The uncertainty of (c) is large essentially because the inclination of the orbital 
planes is small. Eichelberger’s result* is from the direct computation of pertur- 
bations of the elements, but rests on only about six months’ observations. This is 
only about } of the libration period, but as it stands the uncertainty should be valid. 

Samter’s discussiont is interesting historically because it shares with 
Eichelberger’s the result that a mass m’= 1, is definitely too low, and also 
on account of incidental remarks. 

Brouwer’s method] is an up-to-date version of G. W. Hill’s. A periodic 
solution is found by numerical integration over half a synodic period; thus the 
eccentricity of Titan and the libration are neglected. The motion of the apse 
and the libration period do not contain the first power of e’, and the second 
power is neglected by Woltjer. On the other hand, the libration, which is large, 
affects the motion of the apse to an important extent, as pointed out by Samter. 

Woltjer assigns weights }, } to (c) and (e). For the others he says that 
incompleteness of the theoretical computations is the chief source of uncertainty. 
He assigns weights } to (a) and (g), 4 to (b), (d), (f). Adopting these weights 
gives standard error per unit weight 118 on 6 degrees of freedom, and 

1/m' = 4033 +71. 
(The standard error is my computation.) ‘The obvious criticism of this analysis 
is that the errors of the various estimates are largely due to omissions. in 
computation, and they are unlikely to be independent. It seems to me that 
Eichelberger’s method is the only one not open to thi; objection, and that its 
uncertainty cannot be raised to 190 or so without other evidence. The data 
for (d), (e), (f), (g) are included in those for (a), (b), and the methods of (d), (f), (g) 
take less full account of higher terms than (a), (b) do. I think, therefore, that all 
the useful information contained in these solutions is already contained in 
(a), (b), (c), (e). If only we could be sure that the computations for (a), (b) were 
complete, they would also include (e). 
*W.S. Eichelberger, Publ. U.S. Naval Obs. (2), 6, Appendix 1, B1-17, 1911. 


+ H. Samter, Sitz. Preuss. Akad. Wiss., 1051-9, 1912. 
t D. Brouwer, B.A.N., 2, 119-120, 1924. 
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Woltjer does not take m and n’ as data. Instead, he takes nm’ and the mean 
eccentricity of Hyperion. With these and an adopted value of m’ the mean 
motion of pericentre can be calculated, and n is fixed by the consideration that 
the mean motion of 4/—3l’—m must be zero. His expansions are in power of 
p=m't?, and he starts with the numerical values m’!=4033, e=0-104 3, 
amplitude of libration 36°. m’* and e” are neglected. On these grounds, since 
e” is about ;, we might expect that this neglect would not produce errors 
of more than 6 or so in I/m’; but this seems to be far from the case. 

Adapting the results for Enceladus, with necessary changes in notation, we 
should have 

v=Pn—(3n'—4n), vy, =y, = —B'n +(3n'—4n), (1) 
P=3m'n* A, O=3m'n*B, (2) 
e=ky= =, kyo =9, (3) 


y2—v* = Phy. (4) 


Woltjer’s equation for the mean motion of the libration argument 
corresponding to (4) is 
y=1°571 64 n'y/m’, (5) 
with no term in m’. The secular terms fn have not been introduced at the start ; 
a correction for them is introduced only at the very end in dealing with the apse 
motion. In the calculation y has been supposed expanded in a series Lv,py’, 
and the terms of order »” in the period equation are of the form 
2, 24°% Ry OR, 6 
BY de ~~ Op 06’ (6) 
where R,=M/2a, @ is the argument 4/—3/'—m and R, is the part of the 
disturbing function not containing e’. Hence 0R,/0@ contains m’e, and the right 
side of (5) needs a factor (e/0-104 3)!” and is analogous to the right of (4). The 
term v* has been omitted. Now v is the difference between the actual mean 
motion of the apse and the part due to the secular terms in the disturbing function, 
including those due to Titan, and is therefore of order m’. Its neglect is therefore 
consistent with the neglect of m’, but it is not negligible. It can be assessed as 
follows. ‘Tisserand (vol. 4, p. 108) gives for the effect of Titan 
229 
= =m'n (453+ 2 = cos v) (7) 
with V =6, cos V near —1. The ratio of the terms is —6-g1. The actual mean 
value is — 18°-656 2 in a year, to which the Sun, the ellipticity of Saturn, and the 
attractions of the other satellites contribute +0°-278. We have therefore to 
share the difference — 18°-934 in the ratio 1 to — 6-91, and the separate contribu- 
tions are 3°:20 and —22°14. Thus v= —22°-14/1°. From observation 
y =0°562 039 + 0°-000 023/1%, 
whence 
I — v*/y? =0-988 36. 


Then Woitjer’s estimate m’ ==. needs to be revised to 


3986 


ai at note: Se 0-988 36(1 —e) = 2-479 6 x 10-“(I —e), 


3950 
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where €=0°104 3(I +€). (11) 
The uncertainty of the observational values is negligible. That of the 
computations is not, but this correction has probably included the most 
important effect depending on m”?. 

To adapt the equation for ky it might appear that the definition of v should 
again take into account the secular term due to Titan, in which case m’ should 
be multiplied by 

22°*14 
15°-93 
I think however that this would be erroneous. If other disturbances were absent 
Woltjer’s calculation of the apse motion would be correct; hence the part of Bn 
due to Titan has already been taken into account, and his equation should be 
right except possibly that the term, already small, due to the other bodies may 
need a slight change. I therefore write his equation from the apse motion as 
ee 
m= 4080 
The eccentricity is given as 0-104 19 + 0:000 25 (pp. 59, 64), whence 
€ = —0°00I I +0°002 4. (13) 
Direct comparison of the estimates of m’ gives «= +0-0058, which contradicts 
the observed value. If we substitute (13) into (10), (12) we get respectively 
m’ = 2-4823 x 10-4, 2-4483 x 10-4. The discordance is about halved, and the 
mean would be 


(I +e). 


(I +¢€)=2°4510 x 10-41 e). (12) 


m' =2°4053+0°0170, 1/m’=40560+ 28, 
on I degree of freedom. 

The ditfe:ence is still larger than we should normally expect from the terms 
neglected. It should be noticed, however, that as the most important part of 
the perturbations arises near conjunction, the effect of e’* is essentially compared 
with (a—a’)/a, and may easily be of order 1 per cent. It is likely to be more 
important in the libration speed than in the apse motion, since the libration 
requires a further differentiation of the disturbing function. On an average e”* 
will tend to increase negative powers of the distance between the satellites and 
therefore, probably, to increase the perturbations for a given mass of Titan. 
On this ground I think that it is incorrect to take the mean of the two estimates, 
and that allowance for e* might easily change both to about ;,. The 
uncertainty of this extrapolated value can hardly be put under I per cent. 

3. Secular motions.—Calculated expressions for these are given (logarithmic- 
ally) by G. Struve (p. 53) in degrees per year. k/a,” and //a)* are constants 
depending on the figure of Saturn and the ring. In ordinary form the coefficients 
of the constants in fx for the various satellites are given in Table III. 

Tas_e III 
Mimas Enceladus Tethys Dione Rhea Titan 
1°48440 X10" 6:20243 x 10% 2°93777X 10% =1'23530 108 3°837X 10% 2°023 X10 
3°949 X 10° 1002 x 103 3°097 X 10? 7°94 X10 1°265 X10 1°241 X 107! 
2601 « 10° 5°384 x 104 1°467 x 104 oi 


2°947 X 10° es 2592 < 10° 3°744 > 6°842 x 10° bas 
6-792 x 10* 2°885 x 105 ss 1°333 > 1°392 X10 4°059 X 10? 
2°093 x 104 4°717 X 108 1°509 X 10° ay 4°385 x 104 =7°153 X 10? 
5849 x 10° 1-019 x 104 1°862 * 104 5°183 X oll 1°705 X 108 
3°831 x 5°732 X 10° 8-268 « 10? 1:287 > 2°596 x 103 


3°162 X 10? 
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The observed values are given in Table II. For the constants I take the 

following values. 

— =2°5 x 10-7(1+10-* pp); — =7°0 x 10-4, ; 

a a 

my; = 6°69 x 10-8 (I + 0-030); My = 1°46 x 1077 (I +.0°36); 

My, = 141 x 10-8 (I + 0-026); My, = 1°854 x 10-8 (I + 0-033); 

May = 2°5 X 10-8 ug; My, = 2°5 X 10-4(1 + 107 pry) ; 

my, = 2°5 X 10-8 wy. 
In addition we have the equation from the eccentricity of Rhea 

My, = 2°5 x 10-4(0'971 640-068); ws = —2°8 +68. 

In the first place I made a least squares solution from the above data alone. ‘The 
effects of the first four satellites were calculated from the start ; the uncertainties of 
their masses contribute +0-01 to the calculated Bn for Mimas and Enceladus, 
+ 0-017 for Tethys, and + 0-oofor Dione. Thusthese uncertainties can be ignored 
except for Tethys, and for Tethys they can be treated by raising the standard error 
of Bn from 0-031 to 0-036. The contribution of Iapetus for Titan is 8 x 1074 py. 
According to the preliminary estimate 1, is of order I, and it would have to reach 4 
for its contribution to be as large as the standard error for Titan’s apse. At present, 
therefore, we ignore p14. 

The first solution gave a rather large x* and made me suspect that some of 
the smaller uncertainties, especially the very small one for Tethys, might be 
underestimates. This was checked as follows. The maximum displacements 
due to equal changes of w and @ would be nearly proportional to ae and asiny, 
where a is the radius of the orbit and y the inclination of the orbit to Saturn’s 
equator; hence aeo,, and asiny.o, should have comparable values in all cases 
for equally good observing over equal periods. The same should apply to 
aco, and ao,, with extensions to the various librations. With a in seconds of arc, 
the results were mostly between 0”-02 and 0”-03, except for the two for Mimas, 
which were 0”-05 and 0”-07, and for the node of Titan, which gave 0-055. The 
two estimates from Mimas already differ by 2-6 times the standard error of their 
difference, strongly suggesting that these uncertainties need to be doubled; 
but none of the standard errors, when so reduced to acommon basis of comparison, 
is ou standingly smaller than the rest. The high accuracy for Tethys is due 
mainly to the rather large inclination. The large y? was therefore treated by 
doubling all the smaller observational uncertainties. It was not thought necessary 
for ‘Tethys to double the contribution of the uncertainties of the masses of other 
satellites. ‘The equations of condition, so modified, are as follows. I take 
Mo = — 28+ poy My =I + py’. 


Mimas apse 0°37I py + 2°76,’ + O-OIpy + O°00I py = + 1°87 + 0°20 

node O°37I py + 2°76,’ + 0-OT py + 0°00 ps3 = + 1°50 + 0°20 
Tethys node 0°073 49 + 0°22,’ + 0°05p19 + 0°002p3 = + 0°087 + 0-065 
Dione apse 0°03 I py + 0°06," + 0°13 py + 0°003 p43 = + 0°20 + 0°42 
Rhea node O-010py + GOT py’ + 0:006p3 = + 0°12 + 0°08 
Titan apse 0°000 5pl9’ + 0-000 Ip,’ +0°004 3, = +0°0078 +0-0070 

node 0°000 5p + 0-000 Ip,’ +0°004 32 = —0°001 4 +0029 

Rhea eccentricity Hs = —2:8+6°8 
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The normal equations, reduced to unit weight, are as follows. 
872g + 55°20py'+1°32H2 +0°070~3= + 33°5, 
55°20p9 +392°4py" +4:05p2 +0°25¢3 = +236°5, 
I-32Mg + 4°05p,'+1°0874.+0°026p3= + 2°67, 
0°070p9 + 0°25, + 0°026u,+0°029Qu,= + O-180. 
The solution is 


Mo = +0°0+1'9; py’ =+0°60+0°28; po=+0°241°5; pg=+O0-Q+6°1. 


x” on 4 degrees of freedom is satisfactory, and the data for these satellites are 
consistent among themselves. ‘The uncertainty of 6 per cent in the mass of 
Titan, however, calls for improvement, which may be sought in the data for 
Hyperion and Iapetus. 

3.1. Secular movements of Iapetus.—The relevant formulae are in H. Struve’s 
paper of 1888. G. Struve gave data up to 1927, but his paper ends abruptly 
after linear forms have been fitted; presumably he would have brought 
H. Struve’s work up to date had he lived longer. The equations can be 
written as follows. 6 is again written for longitudes on the ecliptic. 

sini 5 = —<XK,{ sinicosi,—cosisinz,cos (6—6,)} 
x {cosicos?,+sinisini,cos(6—6,)}, 
di 
dt 
‘These were integrated numerically ; the values taken for the Sun and for Saturn’s 
equator, referred to the ecliptic and equinox of 1890-0, were 


= XK, sini, sin (6—6,){cosicosi,+sinzsini,cos(@—6,)}. 


09 =112° 42’, ig =2°203’, 
1 = 167° 58-79 + 0'-46, 7, =28°04'-44 +025. 
The last is based on G. Struve’s determinations from the orbital planes of the 
five inner satellites (Heft 4, p. 48) for 1889-25. These means, weighted according 
to the stated uncertainties, gave x7 =3-2 on 8 degrees of freedom; there was 
therefore no occasion for modifying the weights, as Struve did. 6, was adapted 
to 1890-0 by allowing for precession. These values could be taken as constant. 
For Titan the motion of the orbital plane is appreciable. G. Struve’s solutions 
(Heft 5, p. 11) were used. Precession had already been removed from his values 
of 6y; The satellites out to Rhea could be combined with the ellipticity of 
Saturn, but Titan needed separate treatment. (H. Struve had taken its orbital 
plane as identical with Saturn’s equator, but as these planes can be inclined at 
about 38’ and the orbit of Iapetus is inclined at about 10° to both, this might have 
produced an error of order 5 per cent in the mass of Titan.) We have also 
Ko =4':073/1°, Ky, =0"694/1", Ky =2"779(1 +0-0Ips)/1". 

K, has a standard error of about I per cent, mainly arising from the mass of Rhea, 
but this can lead to one of only 0-3 per cent in the mass of Titan. ‘The initial 
values adopted for 1890-0 were @ = 142° 05’,2=18° 24’, and numerical integration 
was carried out to cover the interval 1770 to 1930. The secular motion is slow 
and it was enough to regard corrections to the initial values as carried through the 
entire period. ‘The result was compared with the observed values (adapted to 
the fixed equinox used). These were much less accurate up to 1880-2 than later. 
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I took the standard error of 6 to be coseci times that of 7 in each interval, with the 
following results. 


Crude residuals Revised residuals 


for) 0; Cg OG; 
, 


Up to 1880-2 rr’r 8 3'°6 IrI’-0 35 
1885-1927 20 062 1-7 0'*54 


The unknowns were estimated with these uncertainties and the standard errors 
redetermined from the new residuals. No substantial further correction was 
needed, and the results were 


86(18g90-0) = + 0'-56 +0'-89, 
81(1890-0) = + 1-59 + 0-35, 
fe™ ~FISFt. 
The separate estimates of uw; from @ and 7 were —4-8+4:0 and —5:7+2-4. 
Effectively the uncertainties are on 16 degrees of freedom in comparison with 


usual standards, since two standard errors have been estimated from 20 equations. 
The residuals are as follows. 


Observer O-C 
6 z 


, , 


Herschel and Bernard —20°7 —2°7 
Bessel —5: —4°7 
Jacob +9": —1I-0 
A. Hall +13 +0°2 
A. Hall +7°9 + 3:0 
H. Struve +0°6 +0*2 
H. Struve and Bernewitz +17 —0'7 
H. Struve and Bernewitz +1-0 +04 
Alden — 3:0 —0-6 
Alden +0Or1 —0'4 


The residuals for the second period are substantially greater than the apparent 
observational uncertainties would suggest. Since a similar result was found 
for the inner satellites it looks as if some additional source of disturbance affects 
secular motions. 

It has been indicated above that uncertainty of the mass of Rhea contributes 
one of about 0-3 per cent to the mass of Titan. The uncertainty of the orbital 
plane of Titan is about 0’-8, compared with a mutual inclination of 10°, and 
therefore contributes about 0-1 per cent. The effects of the uncertainty of 
Saturn’s equatorial plane and of changes in the plane of its orbit would be still 
smaller. ‘Thus the above uncertainty is probably valid and we have from lapetus 


I/my = 4243 + 93- 
This value is consistent with Eichelberger’s. If, as suggested above, Woltjer’s 
value from Hyperion should be revised to about ;4,, it also would be consistent, 
but a value near 4; would be hard to reconcile with either the present value or 
Eichelberger’s. 
G. Struve gives also data for the eccentricity and pericentre of Iapetus. The 
pericentre moved 28° in the interval, with reference to a fixed equinox, but 
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inspection of the residuals left by Struve’s linear solution (and indeed the 
uncertainty that he himself gives) shows that the standard error of dw/dt is not 
less than I part in 30. Since only about a third of the movement is due to Titan 
it follows that these data cannot give my, with an uncertainty under Io per cent, 
and detailed computation would be useless. The reason for the lower accuracy 
is tha: e for Iapetus is much less than the sines of the inclinations of the orbital 
planes. 
4. Combination of data.—The solutions for my; are equivalent to the following 
values of ps3. 
x’ 
Inner satellites +0°9+6:1 0-6 
Iapetus —55+21 408 
Hyperion 
Eichelberger —4F1I+13 0-2 
Woltjer (adapted) —29+1'0 O04 


2:0 on 3 df. 


The weighted mean is —3'54+0-74. The uncertainty given to the adapted 
value from Woltjer’s results is, of course, wholly provisional; but it would be 
impossible to reject these results and no better course suggests itself. The best 
method of all, of course, would be to extend the theory of Hyperion to take account 
of e’* and possibly higher terms in my,, but it is already clear that this would be 
an undertaking comparable with the lunar theory. 

The above value of 4, was substituted in the normal equations derived from 
the secular movements of the inner satellites. ‘The solution was 


Mo = +O°2+1°9; py’ = +0°58+0°28; po=+O°I5+1°5. 
Comparing with the data that contain these three parameters gives yx” =4°9 on 
4 degrees of freedom. Then 
k/ ag? =0°024 305 + 0-000 048, 
l/ag =(11-04 + 1-96)10°4, 
My, =(0"4 + 3°8)10-%, 
My, =(2°411 5 + 0°018 5)10-4,7 
I/my,=4147 + 32. 


The theoretical values of the secular motions of the inner satellites were finally 
computed directly from the original formulae with these values. The details 
are as follows; the separate contributions are given to indicate their relative 
importance. The third decimal was kept in the calculation in case small positive 
values might add up. 


Mimas Enceladus Tethys Dione Rhea Titan 
k 360-78 150°75 71°40 30°02 9°33 0°4917 
l 4°36 I'tl 0°34 0°09 O'ol 0'000 I 
Mi — 0°02 0°00 000 
En O04 8 0°04 "00 
Te 0°08 0°33 see O15 0°000 5 
Di 0°04 o'09 0°28 eee 008 0001 3 
Rh 0°00 0°00 O01 0°02 =— 0°000 7 
Ti o'09 O14 0°20 » 0°63 — 
Total 365°40 152°43 72°272 30°60 10°06 0°4943 
Obs f 365-60 0°20 72'227+0°065 30°7510°42 10°20+0°08  0°5012+0°0070 
* | 365:23 +0°20 0°492 +0029 
O-C +020, —0'17 —0'045 +015 +o°14 ++0069, —0:002 
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There are some end-figure changes, but recomputation was not thought 
necessary. G. Struve’s values of k and / (p. 54) correspond to 


k/a,? =0°024 398; l/agt =7-013 x 10-4. 
These do not appear to have been derived from a least squares solution; the 
difference in k is significant and that in/is great. He gets a residual of —0°-26/1% 
for Tethys, in comparison with an apparent “‘ probable error” of 0°-021/1°, but 
passes it with no comment except that increasing the mass of Rhea would make 
it worse. In the present solution this is avoided by having & slightly smaller and 
/ larger, thus fitting both Mimas and Tethys. 

In comparison with the discussion of 2.2, the calculated $a for Enceladus 
gives 5(8m)= —0-34+0-°05. Correcting the masses of Enceladus and Dione 
accordingly, we have 

My, = (1°27 + 0°53)10~*, 
My, = (1°825 + 0-061)10~-8. 


With regard to the still unsatisfactory determination of the mass of Rhea, 
this may be interpreted as saying that there is a probability of about § that it is 
less than sqqq of that of Saturn, one of about } that it exceeds this value and 
3, that it exceeds jig: _As Rhea is brighter than Tethys and Dione a value 
near sg) Would be plausible. 


160 Huntingdon Road, 
Cambridge : 
1953 January 9. 
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Summary 

The theory of the figure of the Earth on the hydrostatic hypothesis, 
developed by Darwin and de Sitter, is applied to models for Jupiter and 
Saturn proposed by W. H. Ramsey and B. Miles. It is found that the 
integral equation for the coefficient of the fourth harmonic is more easily solved 
directly than by conversion into a differential equation. The results show that 
the ratio D/J?, which is equal to 35//12k? in the treatment of H. and G. Struve, 
is substantially larger than for a homogeneous body in a hydrostatic state, and 
in satisfactory agreement with observation for Saturn. Comparison with a 
model given by Darwin indicates that an inequality stated by de Sitter for the 
coefficient of the fourth harmonic is not general. 





1. Darwin, in a classical paper, developed the hydrostatic theory of a 
rotating planet to the second order in the ellipticity.* His definitions of the 
quantities concerned led to considerable algebraic complexity; part of this 
was removed by de Sitter}, but de Sitter published only a few of the formulae 
and made several errors, some of which have been corrected previously, 
especially by Bullard.{ ‘The theory has been applied hitherto only to the Earth, 
but can now be usefully applied also to Jupiter and Saturn, since the fourth 
harmonic in their gravitational poténtials has an observable effect on the motions 
of the nearest satellites. Further, the work of Ramsey and Miles on these planets 
makes it possible to give a theoretical determination of these terms and compare 
with observation. Application of the iirst-order theory to rotating stars is 
perhaps not altogether out of the question. 

The first-order terms in the equation satisfied by the ellipticity lead to 
Clairaut’s differential equation, the numerical treatment of which for the Earth 
is greatly simplified by Radau’s approximation. The second-order terms were 
included by Darwin and de Sitter, who found that they could be treated by a 
simple modification of the equation. Lor Jupiter and Saturn, however, Radau’s 
approximation is seriously wrong and a differential equation has still to be solved 
numerically; and the equation satisfied by the fourth harmonic has so far been 
solved only by numerical integration, except in a few special cases. 

The equations arise in the first place as integral equations, which are 
converted into differential equations and a pair of boundary conditions. It will 
appear in what follows that this is unnecessary, and that for the fourth harmonic 
at least it is easier to apply numerical methods directly to the integral equation. 

In the theory of Darwin and de Sitter it was found convenient to define a 
modified ellipticity, the relation of which to the original ellipticity appears to 

* Sir G. H. Darwin, W..N., 60, 82-124, 1899; Sci. Papers, 3, 78-118. 
t+ W. de Sitter, B.A.N., 2, 97-108, 1924. 
t In addition, in (1), 128/105 should be 32/105; in (11), 32/3 should be 8/3. In the set of 


equations at the foot of the page, S,=3J should read S,=3J(1+4e). (11) is stated to agree with 
Darwin, but Darwin’s form is correct. 


/ 
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contain the coefficient of the fourth harmonic. This surprising result is due to 
the definition of the ellipticity in terms of the ratio of the equatorial and polar 
axes. I have proposed for geophysical purposes that the standard ellipsoid 
should be so defined that, with the outer surface a level surface, it should lead to 
the correct value of the second harmonic in the external field. Then allowance 
for an additional fourth harmonic will not affect this relation. ‘The reasons for 
this course are stronger when the data are purely astronomical. 

The following treatment is much shorter than Darwin’s; on the other hand, 
de Sitter’s is condensed to an extent that makes the intermediate steps very 
difficult to reconstruct, and I think a rather fuller account than his is needed. 

The axes are taken to be a, a, a(I—e); the sphere of equal volume has 
radius 4, satisfying 

6? = a®(1—e). (1) 
The parameter m (de Sitter’s p,) is defined by 
m = w*b3/f M. (2) 
I use the spherical harmonics in the form 
L,=}4-sin* ¢’ = —§P,, 
L,=sin‘ ¢’ —$sin? ¢’ += P,. 
We have the following identities : 
sin? ¢’ =4—L,; sin’ ¢’=L,—§L, +}; 
sin? ¢’ cos? ¢’ = —L,—4L,+2; L2=L,-AL,+¢. 
The equation of the ellipsoid, to order e’, is 
r/b=1— Ze? + (e+ Be*)L, + 3e* Ly. 
A term —x« sin* 2¢ would give one in 4«L,; taking also, with de Sitter, 
e’ =e-Ze, 
we have 
r/b=1— Ae? +(e + e*)L. + (Re? + 4x) Ly, 
differing from (6) only by a fourth harmonic. ‘The quantity e’ is identical with 
de Sitter’s, but, with this definition of e, « does not appear in the definition of e’. 

The internal and external potentials and their first derivatives are continuous 
at the surface. ‘These relations have to be converted into a pair valid at r=d. 
Suppose that /, 0f/Or vanish atr=b+h. Then to the order retained 


fir) =Mr—b—hyf"(b+h) + A(r—b—hy'f"(b+ hy, (9) 
f(n)=(r-—b—-h)f"(b+h) + Y(r-b- vl b+h) 
=(r—b—h)f"(b) + h(r—b—h)f"'(b) + (rv —b-hy f(b +h), (10) 
and to the second order in h, with r=), 
f(b) = 3h f"(b); f(b) = —hf"(6) — Lhe f’"(6). (11) 


‘The external and internal potentials of a homogeneous body are taken as 


ry Z b* b! 
U,=fM {2 +(J +J5) a L,+ Do 3 L,} and 


, gee 
U,= Fab) +f + J Jebel i. 
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Their difference and its first derivative vanish when r satisfies (8) ; adapting these 
conditions to r=b by (11), and treating J as of the first order in e, we find 


a | — > - 16 3,2 4. _— 
J =<e > Jo = xe ° J,= Dy =5e + 3K 5 D, =*. (14) 


162 . 
~ 85 3 
Darwin gets equivalent formulae by using the exact potential of an ellipsoid of 
revolution: The adaptation to a heterogeneous body in a hydrostatic state with 
external mean radius 4 is as usual; if 7, is the mean radius of the surface of equal 
density — the point considered, 
5 a . > Pr,’ a 
=i), pd(3a’*)+ ier i a a 4 +P) 14+ 3Sr 4420 ) bs (15) 


5 r°? 3 r," 


I [* , I tie 75 , 9 
Po= | pda®;  S= 3 pd{a’*(e’ +3e*)} ; 
a’=0 a’ =0 


3 F 
’; r,° 


“Ty 


I ° 9 
pd(e’ + 38e*); P=—| pd{a’7(e + Sx); 


) 


, Md , 
/a@’=r 1 “~a=0 


-b K 
=r, d(=.) > | 
G m7 puns a | 


where in the integrals p, e’ and « are regarded as functions of the variable 
of integration a’. ‘The difference between e’ and e does not affect the terms of 
order e*, in which the accent is therefore omitted. ‘To the same accuracy the 
rotation term is 

kw*r* cos* ¢' = §af pmr’(§ + L,). (17) 


The sum of (16) and (17) must be constant over a surface of constant r’. Then 
using (8), with 6 replaced by 7,, and expressing L,” in terms of L, and L, by (5), 
we find 

(e’ + Se*)py + (2+ Be’)S+ 3 — Ze’) T+ 4 pm+ > pme =o, (18) 


— (he? + 4x)py — Se’ S+ Se’ T+3P+40+pme =o. (19) 
Multiplying (18) by (1 —#e’) gives 
(e' + 3e*)py—3(S+ T) — dpm = fe(pm— 37), (20) 
and multiplying (18) by 2e’ and subtracting from (19) gives 
(36 — 8x)py —6eS + 3P+8O =o. (21) 


(20) and (21) are de Sitter’s equations, except for the correction of an algebraical 
error. ‘The latter can be checked by comparison with the first-order theory of 
the fourth harmonic. 

It is of some interest to show how some of the main results derived from 
the differential equations can be found more easily directly from the integral 
equations. I take the first-order equation, for a harmonic of any order, in the form 


a { I \ , 3 =) - (3 :)} -_ 
Pa =, ce bee d(a ie ae 4 d *n—2 + 0, 3 m). I 
Po 2n+1\r"*3 Pe ( os ¥ (0, Spm). (1) 


This is essentially an integral equation (written by means of Stieltjes integrals) 
with fixed termini. For the general case (zero on the right, 41) we can without 
loss of generality take ¢,, positive at its extreme value; let this be e, taken at r=r’ 
(e,, is continuous by the physical conditions). In (1) replace r by r’ and integrate 
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by parts, which is permissible because ¢, is continuous and p non-increasing. 
Then the right side is 


si “b 
3 {rs > _-——; a’*+*e_dp—r'** wid} , (2) 


2n+I1 yer a’=0 Ja=r a’" 


which is not decreased if ¢,, is replaced by « everywhere (since p is non-increasing). 
Make this substitution and integrate by parts again; then (2) is 


3€ I ‘all » I 
Be —— ae > 
< ye rk | pda’" Sipy'n a | pd i : 
QRN+I FT J nd a’=r’ ais 


The second integral is zero for n=2, or r’=6, otherwise negative for n>2. 
po(r’) is not a function of a’ ; hence the right side is 


3e OT 


~ an+ 17/8 [por’" +? + {pp —po(r’)} da’ *3). (4) 


Unless p is constant there is an r” such that p(a’) <p,(r’) for a’ >r”’, >p,(r’) for 
a'<r"; then the integral is 

n + 3 r\ , ” ’e 
= f{p —po(r’)(a’" —"") da? <0 


and finally the right side of (1) is 
3€ 
— Se 0. (5) 
2n+1°° ” 
But this is less than the left side unless m=1 or e=0. Hence all harmonics other 
than those with n=1 and the main ellipticity term have zero coefficients. 

An argument on these lines is well known to give a proof that e,(b) =0, after 
it has been shown from the differential equation that ¢,, if not zero, must be 
monotonic. ‘The above generalization of the argument does not need the 
preliminary use of the differential equation. 

If n=1, take e,=1/r. Then the right side of (1) is 

b 


ay pda + | pdx) = 


4 
r J a’'=0 


’, (6) 


agreeing with the left side. Hence 1/7 is a solution. 

‘The above argument cannot be used directly to prove that it is the only solution, 
since e, may be unbounded. But put re; =. Then 7 is continuous even at r =0. 
Let its maximum be ¢ at 7’; then the equation (1) is 


: >? F “b 
Pon = p(b)n(b) -— = | a’*n dp — | n dp 
a a’=r 


Ya’=0 


“b ) 
a dp— | dp C 


~a=r 


=<5 | pda™ > ¢ =p. 


(rc 
LT / a’=0 J 

Hence 7(r’)=¢ is false unless the equality sign can be taken in the second line, 
that is, unless 7 is a constant. Hence e, < 1/r is the only solution (with the trivial 
exception that in intervals where p is constant 7 might vary without there being 
any dynamical consequences). 
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For the main ellipticity. term and for the second-order term containing a 
fourth harmonic, the terms of lowest order take the form 


3 {yn sebel) 





Poen = qn 


2n+1 . 
aE ff i € 
‘n+3 U »n—2 
es a’"*3¢ dp’ —r" goad’ +f) (8) 
/a’=0 /a=r 
Take as a trial solution 
f(r) 
ree (9) 
and substitute on the right; then the modulus of the first set of terms, by the 
previous argument, is 


—. 
~ 2n+ ea Nn 7. (10) 


hence the correction is everywhere numerically < = (tn mas and the 


™ 
results of successive approximations converge like z(; a) . Since the 
trial solution is continuous the corrections are continuous, and the limit, being 
the sum of a uniformly convergent series of continuous functions, is continuous ; 
hence all the required conditions are satisfied. 

I have not succeeded in finding a direct proof from the integral equation 
that «, is monotonic for the main ellipticity term. Clearly «, is everywhere >o 
since the first approximation and all corrections are positive. 

The series may converge too slowly for direct computation. In fact if the 
upper bound of the first approximation is «, that of the true solution may be as 


’ 


large as e, which is 2e for n=2 and y¢ forn=3. ‘The former is actually 


2n+I 
2(n —I1) 
attained for a homogeneous body (for which the first approximation is }m). 
This difficulty may be circumvented as follows. Write the equation as 
n= fo + Ke,, (11) 
where K is a linear operator. The first approximation is ¢, =u, =f(r)/pp. 
Substitution of u, on the — gives a second approximation u,. The essence of 


the matter is that {ttm fr) — Ku, 3 forms a sequence of functions, decreasing 
Po 


approximately geometrically. With two values, say m=1, 2, we can choose a 
constant » so that the function 


(I+ p )(u- fe) - Ku;) —p (u.-© -Ku,)| (12) 


is everywhere as small as possible, and this will be done by taking » so that the 
greatest positive and negative values of the expression under the modulus sign 
are equal in modulus. But this amounts to saying that, judged by the residuals 
left in the original equation, 


r ? 
(1 + pity — pe, 2 +K{(1+p)uy— pu; } (13) 
0 
with the smallest possible error. Then (1 +)u.—ju, can be used as the first 


of a second pair of approximations, and solution should be rapid. ‘The method 
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is essentially one of relaxation. In contrast to the use of such methods for 
differential equations, the “point adjustment” needed to smooth out rapid 
fluctuations of the errors does not arise, being automatically taken into account 
in the integration process; all the approximations are “block adjustments ”’. 
On the other hand, it may not be so convenient for most integral equations as 
with the present ones, because these have the feature that the kernel can be 
written in the form 
&1(a')hy(r) + go(a’)h,(7), 

thus containing four functions that can be tabulated once for all; a more general 
type of kernel would have to be fully tabulated as a function of two variables. 

3. Possible variation of D/J*.—The form of a homogeneous rotating body 
would be an exact ellipsoid of revolution. It was shown by H. Struve that his 
parameters k and / would lead to 

l=? 

7 

for all ellipticities, corresponding, in the present notation, to D/J?=%. The 
value found in my discussion of Saturn’s satellites was //k® = 1-87 + 0-33, and it 
is desirable to examine whether the doubtfully significant difference would be 
explicable by variation of density. At first sight this seemed unlikely. For if 
we consider an extreme case where the whole mass is concentrated in a 
homogeneous core, the outer part has no dynamical effect at all, the core would 
take a Maclaurin form, and the same relation would hold. Thus rather wide 
departures from homogeneity might not affect D/J?. Examination of Bullard’s 
solution for the Earth, however, did show an appreciable departure, and direct 
examination of the data for the constitutions of Jupiter and Saturn was desirable. 

3.1. Solution for Fupiter.—Possible distributions of density in Jupiter and 
Saturn have been given by Ramsey and Miles.* They worked out the second- 
order theory of the ellipticity, but not the fourth harmonic. The details were 
not published, but Ramsey and Miles kindly gave me a typed copy. Of their 
models, those denoted by J, and S, agreed best with the observational data; 
but those denoted by J,, Sy were only slightly inferior and have the advantage 
for the present purpose that the distributions are smoother. The object of the 
present work is to find out whether the ratio D/J*, for given m, can depart widely 
from the value #? that it takes for a homogeneous body of any ellipticity. 
Accordingly I accepted their values of derived by de Sitter’s method for 
J4, Sy, derived the distribution of e, and then used the integral equation 2 (21) 
to estimate x. I use their values for the mass, radius, and rate of rotation of each 
planet. The results for Jupiter are in Table I, where m =0-083 39. 

« at most is only about 0-1 e* ; this is mainly because the e* terms in the integral 
equation are (in much of the interval) nearly as 1, —2, 1. The corrections in 
the second approximation never exceeded 0-1, so that the method is even more 
successful than there was apparent reason to expect. ‘This is chiefly because « is 
small except where p is small compared with p. 

For the external field 

b 
pb'J = > | p d{r°(e’ +2e%)}, 
0 


“b 
pb’ D=- | pd{r'(e? + =«)}. 
r=0 


? 
~ 





* W. H. Ramsey, W.N., 111, 427-447, 1951. B. Miles and W, H. Ramsey, M.N., 112, 234-243, 
1952. 
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For comparison with a homogeneous body it seemed best to drop the e* term in-/. 
This is because e* and xe are neglected in D, and we know that for a homogeneous 
body D/J* is independent of e. It therefore appeared possible that inclusion 


of a higher term in J and not in D would actually decrease the accuracy. On this 
basis J =0-0260, D =0-00339, of which only + comes from«. ‘Then 

D/J* =5:0. 
Allowance for the e® term would reduce this to 4:8 or 4-9. 


TABLE I 
r(10* cm) e p ins 10% K« (first 

approximation) 
0:0676 
0°0673 
0:0667 
‘0657 
0645 
0632 
0617 
0605 
"0605 
0599 
"0575 
0555 
"0539 
"0524 
"0513 
"0503 
0492 
0484 
0479 
0475 3°56 
"0472 3°66 


II 
14 
‘22 
27 
-30 
*36 
39 
42 
‘89 
‘28 
‘60 
‘92 
23 
Be 
‘80 
05 
27 


6-99 
6-96 
6-9 
6°8 
6-7 
6°5 
63 
6-10 
6-10 
6:0 
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3-2. Solution for Saturn.—This is given, tor Model 54, m =0-1353, in Table II. 
Further corrections were negligible. ‘The results lead to 
J =0°:0438; D=0-0116; D/J* =6-0. 
The conjecture that D/J* for a body in hydrostatic state is nearly independent 
of the distribution of density is therefore seriously wrong for Jupiter and Saturn. 
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4. A model given by Darwin.—'lo investigate the possible range of values 
of D/J* I consider a case treated by Darwin, of a body of unit density with a 
fraction » of its mass concentrated at the centre, the density otherwise being 
uniform. I take the radius also as 1. Darwin determined the value of yu that 
makes « greatest, m being taken as given, but did not explicitly attend to the 
external field. (My pw is Darwin’s »/(1+).) The results above become 


I 
Po = p wr (1 —p)r}, 


S=(1 —p)r(e' + 3c); T=(1 —p)(e,' +316," — e — He’), 
P=(1-plle +i); Q=(1-p)(kp?-K). 
To the first order 
em ' r® 
so a ey 5? 
I+3p Br bi — py —5(L-p)r? 
and to the second order 





; = 


c= ee 
*  8(2+ p)’ 


J =8S, D=8P, 


D 25 2+2p 


J? © (1—p)(2+p)" 
Hence this ratio actually tends to infinity as p tends to 1, and is never less than 
*°, the value for a homogeneous body. The reason is that as p—1, J is of 
order (1 —y)m and D of order (I — 4)m?. 

This result is of some interest in relation to the models J, and S, of Miles 
and Ramsey, which they thought the most probable. In these the heavy elements 
are supposed concentrated near the centre. For J; the nucleus has about x of 
the mass of the whole body, and we might expect D/J* to be increased by about 
7 per cent above the value found for J,.. For Saturn the corresponding quantities 
are about #4 and 10 percent. It is therefore quite likely that the correct values 
are about 5-4 for Jupiter and 6-6 for Saturn. 

5. Bullard’s solution for the Earth.* —This applies de Sitter’s method. Since 
the differential equation for « is derived from Darwin’s solution, it is correct. 
The boundary condition (10) for « should read 

dk = a: ae 

dp —4K—GPi &1 +76 Pr 
the calculations have been carried out with the correct form. ‘The easiest way ot 
finding D/J* is to compare with the International Formula, which gives 
J =0-00164, D=0-0000107, DjJ*=3°98. Bullard finds on the hydrostatic 
theory that «=6-Sx 10%. ‘This represents the distortion of the level surface, 
and since the International Formula is consistent with the outer surface being 
an exact ellipsoid, irrespective of whether hydrostatic conditions subsist in the 
interior, 4« is added directly to D and makes D/J? = 4-98. 

6. An inequality given by de Sitter.—De Sitter states without proof that for 
surface values 

K <jme — fe’. 





* E. C. Bullard, M.N., Geophys. Suppl., §, 186-192, 1946. 
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This relation is satisfied by Bullard’s solution for the Earth and by the above 
solutions for J, and S,. For Darwin’s model, on the other hand, the right side 
is {ue*, which is less than x, except for ~=0 and »=1. It appears, therefore, 
that the above inequality is not general. 

7. Comparison with observation for Saturn.—The values for k and / for 
Saturn* made k//? = 1-87 + 0°33; to reduce to J and D we must multiply by %, 
giving D/J* =5-45+0-96. Thus the theoretical results for Saturn are consistent 
with observation. There is no basis for an attempt to estimate the mass of the 
Ring, such as might have been made had the observational value exceeded the 
theoretical one. 

On the other hand, the results of the present paper suggest that it is reasonable 
to take for Saturn D/J? =6-3 +0-3, which is equivalent to //k® =2-16 +0-Io. 
This leads to an additional equation of condition, in the notation of the previous 
paper, 

0°027 fy — 7°Opy' = —5°7 £0°58. 
Combining this with the previous data I get 
k/ag” =0°024 270 + 0:000 016; L/ag! =(12-42 + 0-52)10-4; 
Mpy =(2°8 + 2°8)10-*; my, =(2°412 + 0-018) 10~4, 


The comparison with observation is as follows. 


Observation x’ 
Mimas apse 36560 +0°20 °6 
node 365°23 +020 I'l 
Enceladus are 
Tethys node 72°227 +0:065 o'2 
Dione apse 30°75 +0°42 foe) 
Rhea node 10'20 +0°08 $-% 
Titan apse 0*5012+0'0070 0°3 
node 0492 +0:'029 roee) 
us from Rhea —28 +68 5 ole) 
1853 from 
Hyperion and Iapetus —66 +14 ; °'o 
1/k? (theory) 2:16 +0°'10 . o'2 


5°9 (6 d.f.) 


The general comparison is quite satisfactory. In so far as the theoretical value 
taken for //k* is the mean of values for the models S, and S,, and the other data 
indicate a most probable value a shade less than this, it may be held that the 
result slightly favours S,, but the evidence is weak. 


St. John’s College, 
Cambridge : 
1953 january 9. 


* H. Jeffreys, M.N., 113, 81-96, 1953. 
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GLACIAL EUSTASY AND THE ROTATION OF THE EARTH* 


Andrew Young 
(Received 1952 December 29) 


Summary 


A rise in sea-level at the rate of 1 cm per century following a reduction 
in the extent of ice in glaciated regions results in an increase of the moment 
of inertia of the Earth about its axis of rotation sufficient to lengthen the 
day by 10~* second per century provided there is no isostatic compensation. 
The consequent apparent secular acceleration of the Moon is 2 seconds of 
arc per century per century. The changes however are negligible if 
immediate and complete isostatic compensation occurs. A _ tentative 
discussion of climatic and glaciological evidence indicates that fluctuations 
of sea-level do occur and that their effect on the Moon’s apparent acceleration 
may be appreciable. This re-opens the question of the cause of the Moon’s 
apparent acceleration, which is usually attributed to the deceleration of the 
Earth by tidal friction. The change in tidal friction due to a rise in sea-level 
is very small and has not been sufficient to produce the observed change in 
the Moon’s acceleration. 


* The full text of this paper is published in M.N., Geophys. Suppl., 6, No. 7, 1953- 
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OBSERVATIONS OF COLOUR TEMPERATURES OF ST ARS* 


VoL. II. REeLatrveE GRADIENTS 


Royal Observatory, Greenwich 
(Communicated by the Astronomer Royal) 


(Received 1953 March 13) 


Summary 


Relative spectrophotometric gradients are given for 250 stars, comprising 
practically all stars of HD types Oes, B and A, brighter than 4™-5 and N. of 
the equator, as well as a selection of F- and G-type stars. The observations, 
which were made photographically with slitless spectrographs attached to 
the 30-in. and 36-in. reflectors at Greenwich, consist essentially in comparing 
(at the same altitude) the blue/red intensity ratio in the spectrum of a pro- 
gramme star with the corresponding ratio in one or more of 25 standard stars 
well distributed round the northern sky. The plates were photometrically 
standardized in an auxiliary spectrograph using multiple slits calibrated by 
the method of full and half apertures. 

The intensity measures in the stellar continuum were made at two wave- 
length regions : one at 6321 A (harmonic mean), the other either at 4285 A 
(for O and B stars) or 4270 A (for A stars) or 4693 A (for later types). 
Corrections were made to the measured gradients to allow for residual 
atmospheric reddening (the comparisons. being made, however, as far as 
possible at equal altitudes) and for inequality in exposure time between 
star spectrum and calibration spectrum. 

The zero of gradient is defined arbitrarily as the mean gradient of nine 
specified stars of type Ao. An estimate of the colour temperature corre- 
sponding to this zero is given from comparisons between the stars and a 
standard lamp operating at known temperature. This enables any listed 
gradient to be converted to a temperature if desired, but temperature values 
are not tabulated because of the uncertainty of the terrestrial comparisons. 

Detailed comparisons are made with other measures of star colour, 
including earlier gradient determinations, colour equivalents, photoelectric 
colours and photographic colour indices. Absorption lines necessarily 
comprised in the broad wave-length ranges selected by colour filters markedly 
affect the correlation with photographic and photoelectric colours. The 
relation between gradient and spectral type is investigated using both HD and 
Mt. Wilson types; and the cause of the abnormally red colour of certain early- 
type stars is discussed. A special series of observations of § Cephei enables 
a comparison to be made between the light curve and a curve showing the 
variation of gradient with phase. 


* The full text of this publication appears as a separate volume of the same title published 
by H.M. Stationery Office, London, 1952 ; price £1 5s. net. 
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